Engineering  Experiment  Station 
University  of  New  Mexico 
Albuquerque,  New  Mexico 

AD  601  / 17 


EE- 107 


THE  EFFECT  OF  THE  ROUGH  AIR- SEA  INTERFACE  ON 
VERY- LOW-FREQUENCY  AND  EXTREMELY- LOW-FREQUENCY  PROPAGATION 


by 

Kenneth  R.  Weiner* 


May,  1964 


*E.  E.  Department 
University  of  New  Mexico 
Albuquerque,  New  Mexico 


This  work  was  performed  under 
Contract  Nonr  2798(01) (FBM) 


Best 

Available 

Copy 


I 


t 

1 

» 

i 

|  ABSTRACT 

i 

l In  the  search  for  extremely  reliable  electromagnetic 
communication  to  submerged  submarines,  the  question  arose, 
"What  is  the  effect  of  the  roughness  and  irregularity  of 
the  sea  surface  on  the  propagation  of  electromagnetic  waves?" 
The  purpose  of  this  investigation  is  to  obtain  an  engineering 
understanding  of  the  effect  of  the  rough  air-sea  interface 
on  electromagnetic  signals  used  in  communication  to  sub¬ 
merged  submarines . 

The  frequency  of  the  electromagnetic  wave  is  restricted 
to  the  ELF.  or  VLF  range.  In  the  initial  part  of  the  investi¬ 
gation,  the  sea  surface  is  assumed  to  be.  a  two-dimensional 
;  (constant  in  one  variable  or  direction)  sinusoidal  surface; 
later  a  doubly  (three-dimensional)  sinusoidal  surface  is 
considered.  The  source  of  electromagnetic  energy  is  assumed 
to  be  a  plane  wave  with  arbitrary  direction  of  propagation 
and  polarization. 

The  fields  on  the  air  side  of  the  sea  surface  are  com¬ 
puted  with  the  aid  of  the  assumptions  that  the  sea  is  a 
perfect  electric  conductor  and  that  the  sea  surface  is  only 
slightly  rough  (i.e,,  the  maximum  slope  of  the  sea  surface 
iS  .ituch  less  than  l)  ,  The  integral  equations  governing  the 
tangential  magnetic  fields  are  formulated  and  solved.  These 
solutions  show  a  variation  of  the  tangential  magnetic  field 
(of  the  order  of  2  db.  frcm  the  flat  surface  case)  depending 


on  polarization  and  direction  of  propagation  of  the  incident 
plane  wave. 

The  fields  in  the  sea  are  computed  by  assuming  the  tan¬ 
gential  magnetic  field  is  continuous  through  the  air-sea 
interface.  The  method  used  in  these  calculations  is  a  numeri¬ 
cal  one  based  on  finite  differences. 

Both  from  the  numerical  solutions  and  a  heuristic  theory 
of  propagation  in  the  sea,  it  is  seen  that  the  perturbation 
of  the  fields  caused  by  the  roughness  of  the  sea  surface 
decays  rapidly  with  depth  if  the  sea  wave  wavelength  is  less 
than  or  the  order  of  magnitude  of  the  skin  depth  of  the  sea 
at  the  frequency  considered;  if  the  sea  wave  wavelength  is 
many  orders  of  magnitude  larger  than  the  skin  depth,  there 
is  little  decay  (at  the  depths  considered)  of  the  perturbation, 
so  that  the  phase  and  amplitude  of  the  fields  in  the  sea  vary 
with  the  height  of  the  sea  vertically  above  them. 
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PREFACE 


This  report  is  concerned  with  a  problem  that  arises  in 
the  theory  of  communication  to  submarines.  In  the  search 
for  extremely  reliable  radio  communication  to  submerged  sub¬ 
marines/  the  question  arose/  "What  is  che  effect  of  the 
roughness  and  irregularity  of  the  sea  surface  on  the  propa¬ 
gation  of  radio  waves?"  Within  the  content  of  the  submarine 
communication  problem  this  report  will  attempt  to  answer 
that  question. 

The  effejt  that  nonuniform  or  rough  surfaces  have  on 
electromagnetic  propagation  is  not  completely  known.  The 
interaction  of  such  surfaces  with  incident  electromagnetic 
energy  is  a  particularly  difficult  problem.  No  attempt  to 
solve  the  above  problem  (where  the  rough  surface  is  taken 
to  be  the  sea  surface)  for  all  frequencies  and  classes  of 
surfaces  appears  feasible  at  present.  However,  for  restricted 
frequency  intervals  and  classes  of  surfaces,  detailed  solu¬ 
tions  may  be  obtained. 

The  radio  wave  propagation  problem  associated  with 
long  range  communication  to  submerged  antennas  is  necessarily 
concerned  with  "low  frequencies."  This  follows  in  part  from 
consideration  of  the  attenuation  of  an  electromagnetic  wave 
as  it  propagates  through  a  conducting  medium.  For  propaga¬ 
tion  through  sea  water,  the  attenuation  is  approximately 
proportional  to  e~  '  ,  where  d  is  the  depth  belc/tf  the  sea 


1 


2 


surface  and  6  is  the  skin  depth  of  the  sea  water  at  a  fre- 

_  wo 

quency  of  the  radio  wave  (6  ~  f~  for  the  frequencies 
considered  here).  The  above  implies  the  need  only  to 
consider  a  restricted  frequency  range  for  the  submarine 
communication  problem. 

With  this  in  mind/  a  study  of  the  electromagnetic 
fields  caused  by  incident  VLF  and  ELF'1'  plane  waves  on  the 
rough  air-sea  interface  has  been  made  and  is  presented  in 
this  report.  Main  attention  is  given  to  the  fields  in  the 
sea,  somewhat  near  the  air-sea  interface  (i.e.,  within 
twenty- five  meters  or  so  of  the  sea  surface),  as  these 
are  the  electromagnetic  fields  presently  in  use  in  communi¬ 
cation  systems.  In  the  VLF  range,  6  for  sea  water  is  a  few 
meters ?  this  implies  the  fields  at  a  depth  of  tens  of  meters 
(a  few  skin  depths)  are  orders  of  magnitude  less  than  the 
fields  at  the  sea  surface.  For  moderate  sea  states,  the 
sea  wave  heights  are  the  same  order  of  magnitude  as  the 
skin  depth  in  the  VLF  range?  the  above  implies  the  radio 
signal  is  greatly  changed  as  it  propagates  downward  through 
the  rough  sea,  as  compared  to  the  signal  under  the  flat 
sea  surface  condition.  For  radio  waves  in  the  ELF  range. 


VLF  (very  low  frequency)  is  usually  taken  to  be  the 
range  from  5  to  J>0  kc/s?  ELF  (extremely  low  frequency)  is 
from  about  1  cps  to  kc/s.  (James  R.  Wait,  Electromagnetic 
Waves  in  Stratified  Media  [Pergamon  Press,  New  York,  1§62 ], 

p.  s.)  . 
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the  skin  depth  is  the  order  of  ten  to  one  hundred  meters 

so  that  the  effect  of  sea  roughness  will  be  somewhat  less 

2 

than  in  the  VLF  range. 

It  is  clear  then  that  an  effect  of  the  rough  air-sea 
interface  is  to  change  or  distort  the  radio  signal  as  it 
propagates  to  a  submerged  antenna.  That  is,  the  rough 
air-sea  interface  has  the  effect  of  introducing  noise 
into  the  radio  signal.  One  of  the  major  problems  in  com¬ 
munication  systems  is  to  preserve  as  good  a  signal  to  noise 
ratio  as  necessary  for  detection  of  the  signal.  In  some 
cases,  as  possibly  (under  certain  conditions)  the  one 
considered  here,  the  noise  properties  of  the  communication 
channel  are  determined  primarily  by  the  propagation  proper¬ 
ties  of  the  time-varying  signal  path  used  in  the  system; 
that  is,  the  noise  created  by  the  time- varying  path  is 
greater  than  the  other  noise  created  in  the  system;  e.g., 
atmospheric  noise,  and  the  noise  created  by  the  time- 
varying  path  is  the  limiting  factor  in  the  communication 
system.  Fart  of  the  propagation  path  to  a  submerged 
antenna  is  through  the  time- varying  rough  air- sea  interface. 
One  of  the  objectives  of  this  study  is  to  gain  an  engineer¬ 
ing  understanding  of  the  distortion  of  the  electromagnetic 
signal  by  the  rough  air- sea  interface,  so  that  ways  of 
alleviating  this  condition  may  be  found. 

- - 

In  sea  water  at  a  depth  of  ten  meters  the  attenuation 
of  an  electromagnetic  wave  at  18.6  kc/s  is  approximately 
47  db.,  while  the  attenuation  at  cps  is  approximately 
.6  db.  A  curve  of  skin  depth  vs.  frequency  is  given  on 
page  55 . 
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There  are  two  major  difficulties  in  the  solution  of 
the  electromagnetic  wave  sea  surface  interaction  problem. 

The  first  is  the  solution  of  the  electromagnetic  boundary 
value  problem  for  a  particular,  completely  specified  rough 
surface.  In  what  follows,  as  is  true  in  most  discussions 
of  boundary  value  problems  involving  rough  surfaces,  only 
an  approximate  solution  is  obtained  for  the  particular 
rough  surface  considered. 

The  second,  and  in  many  respects  a  more  difficu3-t 
problem,  is  the  mathematical  description  of  the  sea  surface.-' 
The  description  of  the  sea  surface  is  statistical,  and  this 
implies  that  the  solution  to  the  "sea  surface- radio  wave" 
problem  would  be  given  in  statistical  terms.  However,  as 
discussed  later  in  this  study,  statistical  results  such  as 
average  field  strengths  are  not  very  meaningful  for  this 
problem. 

In  the  initial  part  of  the  following  investigation, 
the  sea  surface  will  be  assumed  to  be  two-dimensional 
(constant  in  one  variable  or  direction).  The  "basically 
spherical"  earth  is  replaced  by  a  "basically  flat"  earth. 

This  approximation  is  made  often  in  "low  frequency"  propa¬ 
gation  problems,  particularly  when  only  local  fields  are 
4 

considered. 

- 

"The  mathematical  descriptions  of  the  sea  surface 
are  considered  in  Appendix  A. 

4 

Anderson,  W.  L.,  "The  Fields  of  Electric  Dipoles  m 
Sea  Water  —  The  Earth-Air- Ionosphere  Problem,"  Technical 
Report  EE- 88,  Engineering  Experiment  Station,  University  of 
New  Mexico,  Albuquerque,  N.  M. ,  May,  1965/  P*  5* 
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The  theoretical  foundation  for  the  two-dimensional  problem 
is  presented  in  Chapter  1,  starting  directly  from  Maxwell's 
equations  (a  set  of  vector  partial  differential  equations, 
boundary  and/or  interface  conditions,  and  source  conditions). 
It  is  shown  that  for  the  problem  considered.  Maxwell's 
equations  may  be  approximated  by  a  linear  and  time  invari¬ 
ant  operator,  in  which  case  it  is  convenient  to  consider 
only  monochromatic  electromagnetic  fields.  The  concept  of 
vector  potentials  is  then  given,  along  with  a  brief  outline 
of  their  theory  for  the  monochromatic  case.  In  the  two- 
dimensional  problem  the  vector  electromagnetic  boundary 
value  problem  may  be  reduced  to  a  set  of  scalar  boundary 
value  problems  by  use  of  the  vector  potentials.  The  scalars 
used  in  this  reduction  are  the  rectangular  components  of 
the  vector  potentials. 

The  formulation  of  the  scalar  boundary  value  problem 
in  terms  of  integral  equations  is  given.  The  starting  point 
in  this  formulation  is  Green's  theorem,  involving  two  arbi¬ 
trary  functions.  One  of  the  functions  is  restricted  until 
it  is  the  desired  solution  of  the  boundary  value  problem; 
the  second  function  is  chosen  to  facilitate  interpretation 
of  the  mathematical  formulation  in  physical  terms.  An 
unfortunate  result  of  considering  a  plane  wave  as  the  source 
term  and  an  infinite  rough  plane  as  the  scatterer  is  that 
Sommerfeld's  radiation  condition  is  not  sufficient  to  render 
the  solution  unique. 


To  render  the  problem  "well-set, "  a  dstailed  discussion 
of  the  boundary  value  problem  in  terms  of  integral  equations 
is  given.  The  plane  wave  source  condition  is  obtained  as 
a  limit  of  the  usual  source  condition  with  a  finite  source. 
Also,  because  the  virtual  sources  which  are  assumed  to 
exist  on  the  scattering  surface  are  of  unbounded  extent 
as  the  scattering  surface  is  unbounded,  the  radiation  con¬ 
dition  is  imposed  again  by  way  of  a  limiting  process.  With 
the  above  mathematical  formulation,  the  problem  is  then 
"well-set"  and  the  solution  unique. 

To  compute  the  fields  in  the  air  and  on  the  sea  surface, 
the  sea  is  assumed  to  be  a  perfect  electric  conductor. 

Later,  when  considering  the  fields  in  the  sea,  the  electri- 
cal  properties  of  the  sea  are  assumed  to  be:*' 

i)  er  =  8l — relative  permittivity, 

ii)  a  =  4  mhos/m — electric  conductivity  of  the  sea, 

iii)  |i  =  1 — relative  permeability. 

The  above  assumption  of  a  surface  impedance  of  zero  is  a 

usual  approximation  made  in  discussing  the  rough  surface 

£ 

problem;  its  validity  is  discussed  somewhat  later  in  this 


^Stratton,  .T.  A.,  Electromagnetic  Theory,  McGraw-Hill 
Book  Company,  Inc.,  New  York,  1941,  p.  606. 

"Sea  Water"  in  McGraw-Hill  Encyclopedia  of  Science  & 
Technology,  Vol.  12  (Mc-Graw~Hill,  New  York,  i960),  p.  106. 

c 

uMuch  of  the  work  on  "electromagnetic  rough  surface" 
problems  uses  the  assumption  that  the  surface  impedance  is 
zero;  however,  this  assumption  is  not  usually  verified. 

(See  Lerner  and  Max,  "Very  Low  Frequency  and  Low  Frequency 
Fields  Propagating  Near  and  Into  a  Rough  Sea,"  a  paper 
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report.  The  two-dimensional  vector  electromagnetic  problem 
may  then  be  reduced  to  a  set  of  uncoupled  scalar  problems 
for  the  rectangular  components  of  the  vector  potentials. 

The  problem  subdivides  into  two  parts ,  depending  on  the 
polarization  of  the  incident  electromagnetic  wave.  Foi 
different  polarizations ,  different  vector  potentials  are 
used  in  the  formulation  of  the  integral  equations  (i.e., 
scalar  boundary  value  problems)  and  this  is  reflected  in 
the  different  boundary  conditions  applied  to  the  scalar 
field  considered. 

As  stated  above,  the  interaction  of  electromagnetic 

waves  and  rough  surfaces  still  remains  an  unsolved  problem. 

The  theoretical  treatment  of  "rough  surface  problems"  was 

begun  by  Rayleigh  in  his  classic  Theory  of  Sound. ^  A  review 

of  major  theoretical  investigations  is  given  in  Appendix  B. 

An  extensive  bibliography  may  be  found  in  Lysanov's  review 

0 

work  on  Bechmann  and  Spizzichino 1 s  monograph  on  scattering 


presented  to  the  URSI  Spring  19 65  Meeting;  R.  E.  Hiatt, 

T.  B.  A.  Senior,  and  V.  H.  Weston,  "Surface  Roughness  and 
Impedance  Boundary  Conditions,"  in  "Studies  in  Radar  Cross 
Section  XL, "  Ann  Arbor,  Michigan,  The  University  of  Michigan 
Research  Institute,  July  i960,  an  unpublished  report; 

S.  P.  Morgan,  "Effect  of  Surface  Roughness  on  Eddy  Current 
Losses  at  Microwave  Frequencies,"  Journal  of  Applied  Physics, 
Vol.  20,  1949,  p.  552.) 

7 

'Rayleigh,  J.  W.  S.f  The  Theory  of  Sound,  Vol.  II, 

Dover  Publications,  New  York,  1945/  PP*  89“9&* 

0 

°Lysanov,  Y.  P.,  "Theory  of  the  Scattering  of  Waves  at 
Periodically  Uneven  Surfaces,"  Soviet  Physics  Acoustics, 

Vol.  4,  No.  1  (Jan. --March,  195&)/  PP*  1-8. 
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of  waves. ^  It  may  be  stated  that  there  exists  no  theoreti¬ 
cal  solution  for  the  "rough  surface  problem."  However, 
under  certain  assumptions  about  the  rough  surface,  there 
have  been  developed  methods  for  approximate  calculation  of 
the  fields.  In  the  problem  considered,  the  sea  surface  is 
assumed  to  be  sinusoidal  which  is  a  particular  realization 
of  the  sea  surface  (for  brief  mathematical  discussion  of 
the  sea  surface,  see  Appendix  A),  with  wavelength  L  much 
less  than  A,  the  wavelength  of  the  electromagnetic  wave  in 
free  space.  The  technique  used  in  the  calculation  of  the 
vector  potential  on  the  sea  surface  is  an  "integral  equation" 
type  method.  The  physical  parameters  or  constants  of  the 
sea  surface  and  electromagnetic  wave  are  such  as  to  permit 
accurate  calculation  of  the  vector  potential  and  hence  the 
electromagnetic  fields  by  this  method. 

The  solution  of  the  integral  equations  for  the  fields 
or  potential  in  the  air,  but  on  the  sea  surface,  is  considered 
in  Chapter  Two.  In  the  TM  case, ^  the  approximate  method 
of  degenerate  kernels  is  used.  Because  of  the  relative 
magnitude  of  the  physical  parameters  involved,  the  integral 

^Beckmann,  Peter  and  Andre  Spizzichino,  The  Scattering 
of  Electromagnetic  Waves  from  Rough  Surfaces ,  P e rgamon  P r ess , 
New  York,  I963,  PP*  476- 491- 

■^TM — transverse  magnetic  (the  magnetic  field  is  perpen¬ 
dicular  to  a  fixed  direction);  TE — transverse  electric  (the 
electric  field  is  perpendicular  to  a  fixed  direction).  This 
notation  is  explained  in  Harrington,  Time-Harmonic  Electro¬ 
magnetic  Fields,  (McGraw-Hill  Book  Company,  Inc.,  New  York, 
1961),  p.  219/  and  is  used  in  Section  1.2. 
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equations  with  the  degenerate  kernel  accurately  approximate 

the  complete  integral  equation.  (This  result  is  due  to 

Lysanov  and  Meecham.  )  The  approximate  integral  equation 

can  be  solved  by  classical  Fourier  methods.  The  solution 

shows  a  change  or  perturbation  in  the  "worst  case"  for  the 

physical  parameters  considered  of  about  thirty- five  per  cent 

13 

from  the  flat  interface  case.  This  is  in  basic  agreement 

14 

with  Lerner  and  Max,  who  obtain  a  similar  result  by  a 
completely  different  method.  The  above  result  is  also 
shown  to  be  independent  of  incident  angle  (except  glancing 
angle,  which  is  not  directly  considered).  Lerner  and  Max 
considered  only  glancing  angles. 

In  the  TE  case,  the  vector  potential  and  the  magnetic 
field  are  unperturbed,  which  again  is  in  agreement  with 
Lerner  and  Max,  Wait^  and  Morgan.1^  Both  Wait  and  Morgan 

1  Lysanov,  Y.  P.,  "An  Approximate  Solution  of  the 
Problem  of  Scattering  of  Sound  Waves  from  an  Irregular 
Surface,"  Soviet  Physics  Acoustics,  Vol.  2,  1956/  p*  190* 

12 

Meecham,  W.  C.,  "Fourier  Transform  Method  for  the 
Treatment  of  the  Problem  of  Reflections  of  Radiation  from 
Irregular  Surfaces,"  J.  Acoust.  Soc.  Amer.,  Vol.  28  (May, 
1956),  p.  370. 

1^The  tangential  magnetic  field  has  a  variation  of 
approximately  2.8  db.  compared  to  its  constant  value  in 
the  flat  interface  case. 

14 

Lerner,  R.  M.  and  J.  Max,  op.  ext.,  p.  19. 

1-^Wait,  J.  R.,  "The  Calculation  of  the  Field  in  a 
Homogeneous  Conductor  with  a  Wavy  Interface,"  Proc.  IRE, 

Vol.  47,  No.  6  (June,  i960),  p.  1155* 

1 

°Morgan,  S.P.,  op.  cit. ,  p.  353* 
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assumed  this  result  based  on  physical  principles.  Wait 
assumed  the  field  was  unperturbed  in  both  cases,  TE  and  TM; 
Morgan,  however,  solved  the  TM  case,  just  as  Lerner  and 
Max,  by  conformal  mapping  of  the  static  (o>  =  0,  where  cjd  is 
the  radian  frequency  of  the  radio  wave)  problem.  It  is 
interesting  that  Morgan  was  considering  losses  in  "rough 
wave  guides"  in  the  microwave  frequency  range,  which  shows 
that  the  above  results  depend  only  on  the  ratio  of  relative 
physical  parameters,  basically  the  ratio  of  L  to  6,  assum¬ 
ing  ?\  »  L,  and  not  necessarily  on  their  absolute  values. 

To  consider  the  fields  in  the  sea,  the  integral  equa¬ 
tions  for  the  vector  potential  in  the  sea  but  on  the  sea 
surface  is  given.  In  the  VLF  range  for  the  physical 
parameters  considered,  this  integral  equation  yields  an 
"impedance  type"  relation  for  the  vector  potential  much 

like  Leontovich's  impedance  boundary  condition  for  the 
17 

fields.  1  However,  as  this  relationship  applies  to  any 
wave  function,  it  also  may  be  applied  directly  to  the 
electromagnetic  fields.  The  major  phenomenon  leading  to 
this  result  is  the  great  attenuation  of  radio  waves  in  sea 
water.  In  the  ELF  range  the  attenuation  is  not  as  great  in 
terms  of  physical  distances  and  the  "local  impedance  condi¬ 
tions"  need  not  hold,. 


17 

1 Leontovich,  M.  A.,  "Approximate  Boundary  Conditions," 
Investigations  on  Radio  Wave  Propagation,  Part  II,  Moscow: 
Printing  House  of  Academy  of  Sciences,  1948,  PP*  5~^2. 
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By  use  of  the  impedance  boundary  condition  the  effect 
of  the  finite  conductivity  of  the  sea  on  the  fields  may  be 
estimated.  The  solution  originally  obtained  (under  the 
assumption  of  infinite  conductivity)  is  seen  to  approxi¬ 
mately  satisfy  the  "complete"  integral  equation  (where  cp 
and  cp  are  assumed  to  be  related  by  the  impedance  boundary 
condition  —  cp  being  the  solution  to  the  boundary  value 
problem  under  consideration  and  cpn  the  normal  derivative 
of  cp  on  the  surface)  which,  implying  no  basic  change,  is 
necessary  in  the  solutions. 

The  next  section  of  Chapter  Two  is  concerned  with  the 
use  of  quasi- stationary  kernels  and  solutions.  It  is  shown 
that  cp  ,  the  "scattered  field, "  may  be  considered  a  quasi- 
stationary  field,  while  cp^.-  the  incident  field,  and  cp^, 
the  reflected  field  (the  field  reflected  if  the  rough  sea 
surface  were  assumed  to  be  flat),  may  not.  The  use  of 
stationary  kernel  is  valid  in  the  computation  of  cp  ,  but 
not  for  cp.  and  cp 

Now  knowing  the  fields  on  the  surface,  we  wish  to 
obtain  the  fields  in  the  sea.  This  is  then  a  Dirichlet 
type  boundary  value  problem  for  the  rectangular  components 
of  the  fields.  In  some  cases  the  complete  field  may  be 
generated  by  the  use  of  only  one  component  of  the  field  and 
its  derivatives.  In  this  case,  we  may  again  simply  use  a 
scalar  component  of  the  vector  potential.  In  general,  the 
vector  potential  cannot  be  used  to  go  through  the  rough 
interface  correctly  as  there  are  too  many  requirements  on 
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the  potential  at  a  rough  interface.  In  this  case,  each  com¬ 
ponent  of  a  field  vector  may  be  calculated  separately. 

In  Chapter  Three,  the  method  used  to  compute  the  fields 
in  the  sea  is  given.  First,  a  general  discussion  of  classi¬ 
cal  methods,  separation  of  variables,  is  given  along  with 
the  results  to  be  expected  from  such  considerations. 

Basically,  the  above  method  implies  that  if  the  varia¬ 
tion  of  the  electromagnetic  field  with  respect  to  the  x 
variable  is  "too  great, "  the  wave  will  be  attenuated  as  it 
travels  downward  in  the  y-direction.  In  order  t^  see  if  a 
particular  "mode"  will  be  attenuated  significantly  more 
than  the  n  =  0  mode  (plane  wave  propagating  approximately 
in  the  y-direction,  with  propagation  constant  y),  a  crude 
breakpoint  is  chosen.  This  is  if 

L  £  JT~  7rn6 

the  mode  will  be  attenuated  significantly  more  than  the 
n  =  0  mode.  L  is  the  sea  wave  wavelength  for  the  sinusoidal 
sea  wave  considered,  and  n  is  the  index  of  the  mode.  The 
larger  the  n  is,  the  greater  the  variation  of  the  field  in 
the  x-direction  and  the  greater  the  attenuation.  The  above 
implies  that  the  asymptotic  fields  in  the  sea  (i.e.,  the 
far  fields)  tend  to  a  plane  wave  propagating  basically  in 

v 

the  y-direction  with  propagation  constant  y.  Asymptotically 

the  major  perturbation  is  the  n  =  1  mode.  This  result  is 

independent  of  the  shape  of  the  sea  surface;  that  is.  Wait's 

18 

conclusion  with  respect  to  this  result  is  correct,  but  the 

- TR - 

'  °Wait,  loc.  cit. 
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asymptotic  perturbation  does  not  have  the  same  shape  as 
the  sea  surface;  the  latter  is  true  only  for  sinusoidally 
shaped  surfaces. 

For  the  fields  near  the  sea  svrface,  if  the  lower 
order  modes  are  not  attenuated  much  more  than  the  n  =  G 
mode,  there  should  be  little  difference  between  the  fields 
pred  ,ed  by  Wait's  approximations  and  the  actual  fields. 

If  the  lower  order  modes  are  attenuated,  there  will  be  a 
great  difference  and  the  field  will  "rapidly"  become 
approximately  a  plane  wave.  Practically,  the  large  sea 
waves  have  wavelengths  so  long  that  for  much  of  the  ELF 
and  all  of  the  VLF  range,  the  lower  order  modes  are  "un¬ 
attenuated"  ana  Wait's  prediction  is  "relatively  accurate." 
The  smaller  sea  waves  superimposed  on  the  lar^e  ones  may 
have  small  enough  wavelengths  so  that  the  field  perturba¬ 
tions  caused  by  them  will  be  greatly  changed  by  attenuation. 
However,  the  perturbations  caused  by  the  smaller  waves  are 
relatively  small  and  their  effect  on  the  total  field  is 
further  reduced  by  their  additional  attenuation  (over  the 
attenuation  of  the  n  =  0  mode). 

The  above  theory  gives  a  good  basic  understanding  of 
the  effect  of  the  roughness  of  the  sea  surface  on  electro¬ 
magnetic  propagation  used  in  communication  to  submerged 
antennas.  However,  to  place  a  more  precise  meaning  on  the 
term  "little  change"  in  the  fields,  the  numerical  solution 
to  the  propagation  problem  in  the  sea  was  obtained.  When 
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the  numerical  solutions  are  compared  with  theory,  a  good 
understanding  of  the  propagation  problem  results. 

The  method  used  for  numerical  solution  of  the  Dirich let 
boundary  value  problem  is  the  "method  of  lines."  This  is 
a  modification  of  the  classical  separation  of  variables. 

One  of  the  reasons  it  was  chosen  is  its  close  correlation 
with  the  physical  processes  involved  in  the  problem. 

A  finite  difference  approximation  is  used  in  the  x-direction, 
basically  in  the  direction  along  the  surface,  to  determine 
the  "propagation  constant"  in  that  direction.  This  "propa¬ 
gation  constant"  should  be  somewhat  smaller  than  the 
"propagation  constant"  in  the  y- direction.  The  partial 
differential  equations  then  become  a  differential  equation 
in  y  which  is  basically  normal  to  the  surface  with  the 
"propagation  constant"  determined.  The  differential 
equation  is  then  solved.  Here  again  the  rough  surface 
causes  problems.  Below  the  surface,  that  is,  below  a  plane 
tangent  to  the  surface  at  the  lowest  point  (trough)  of  the 
surface,  the  radiation  condition  may  be  applied  and  the 
computation  is  straightforward.  However,  above  this  the 
radiation  condition  does  not  apply  and  the  solution  is  by 
iteration. 

The  results  of  the  calculations  are  given  in  Chapter 
Four,  along  with  interpretation  of  the  calculations.  The 
results  agree  with  what  is  intuitively  expected  by  use  of 
the  classical  arguments  given  above. 
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In  part  two  of  this  study,  the  three-dimensional  elec¬ 
tromagnetic  wave  rough  sea  surface  interaction  problem  is 
considered.  The  method  used  for  the  three-dimensional 
problem  parallels  that  used  in  the  two-dimensional  problem; 
for  this  reason  the  development  given  in  part  two  is  brief 
and  refers  to  the  parallel  development  used  in  part  one. 

The  mathematical  model  of  the  rough  sea  surface  is  a  three- 
dimensional  extension  of  the  model  used  in  the  two-dimensional 
problem. 

The  vector  potentials  are  no  longer  useful  and  the  prob¬ 
lem  is  formulated  directly  in  terms  of  the  electromagnetic 
fields  themselves.  The  Stratton-Chu  equations  are  used  to 
represent  the  solutions  of  Maxwell's  Partial  Differential 
Equations.  The  boundary  conditions  (assuming  the  sea  surface 
is  a  perfect  electric  conductor)  are  placed  in  the  Stratton- 
Chu  equations  yielding  the  integral  equations  (or  vector 
integral  equation)  to  be  solved. 

In  Chapter  J,  the  integral  equations  are  solved.  The 
kernel  of  the  integral  equations  is  assumed  to  be  the  static 
(o,N  =  0)  kernel  and  the  kernel  is  further  approximated  by  a 
simpler  kernel,  permitting  solution  by  the  Fourier  Method, 
as  was  done  in  part  one.  The  solution  yields  the  tangential 
magnetic  field  on  the  sea  surface  which,  as  the  surface 
becomes  constant  in  one  variable,  tends  to  the  solution 
obtained  in  the  two-dimensional  problem.  No  numerical  solu¬ 
tions  are  obtained  in  the  sea  as  the  propagation  properties 
of  the  fields  in  the  sea  are  basically  unchanged  for  the  two 
or  three-dimensional  problem. 


1.0  THE  BASIC  THEORY  FOR  THE  TWO-DIMENSIONAL  PROBLEM 


1.1  Introduction 

The  purpose  of  this  chapter  is  to  present  the  basic  theory, 
assumptions,  restrictions,  and  approximations  involved  in  the 
treatment  to  be  given  below  of  the  question  "What  is  the  effect 
of  the  roughness  and  irregularity  of  the  sea  surface  on  the 
propagation  of  VLF  and  ELF  electromagnetic  waves?" 

The  starting  point  in  the  mathematical  description  of 
macroscopic  electromagnetic  phenomenon  is  Maxwell's  equations.'*' 
Maxwell's  equations  are  a  set  of  partial  differential  equations, 
boundary  and/or  interface  conditions  which  the  electromagnetic 
fields  must  satisfy.  A  general  review  of  Maxwell's  equations 
is  presented.  Then  the  electromagnetic  fields  are  restricted 
to  be  monochromatic  along  with  the  assumption  that  the  sea 
surface  is  stationary.  In  a  later  chapter,  after  the  boundary 
value  problem  has  been  solved  for  a  stationary  sea  surface, 
the  effect  of  the  motion  of  the  sea  surface  is  considered. 

The  vector  potential  method  is  introducted,  and  it  is 
shown  that  some  electromagnetic  boundary  value  problems  may 
be  formulated  in  terms  of  the  vector  potential . 

The  concern  here  is  only  with  "large-scale"  phenomena 
and  in  the  past  it  has  been  confirmed  that  solution  of  Maxwell's 
equations  does  represent  the  actual  measurable  quantities.  (See 
Stratton,  op.  cit , ,  p.viiand  Harrington,  op.  cit . ,  p.  1  ) 
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A  mathematical  model  of  the  physical  situation  involving 
the  rough  sea  surface  is  given,  along  with  a  brief  comment  on 
the  validity  of  the  assumptions  and  approximations  involved  in 
the  construction  of  this  model.  Later,  some  of  these  assumptions 
are  again  reviewed  to  give  a  further  estimation  of  the  accuracy 
of  the  solution. 

Using  as  a  basis  the  mathematical  model  just  described, 
the  electromagnetic  boundary  value  problem  is  formulated  in 
terms  of  integral  equations.  Starting  with  Green's  theorem 
for  two  relatively  arbitrary  functions,  restrictions  are  placed 
on  one  of  the  functions  so  that  it  represents  the  solution  of 
the  boundary  value  problem  considered.  The  other  is  an  aux¬ 
iliary  function  arbitrarily  chosen  as  the  "free-space  Green's 

i 

function"  which  permits  immediate  interpretation  of  the  integrals 
in  terms  of  physical  processes.  The  two-dimensional  electro¬ 
magnetic  boundary  value  problem  is  formulated  in  terms  of 
scalar  components  of  the  vector  potentials,  yielding  the  inte¬ 
gral  equations  to  be  solved  in  the  next  chapter. 

1.2  Review  of  Maxwell's  Equations  and  Monochromatic  Fields 

In  a  source-free  region,  the  set  of  partial  differential 

2 

equations  included  in  Maxwell's  equations  is 


curl  e  =  -  Id 

(a) 

curl  it  =  t  ^ 

(b) 

div  cf  =  p 

(c) 

div  it  =  0 

<a) 

In  rationalized  m.k.s.c.  units  which  are  used  throughout 
the  remainder  of  this  paper  unless  specifically  stated  otherwise. 
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The  constitutive  relations  which  characterize  the  electro¬ 
magnetic  properties  of  the  medium  in  which  the  electromagnetic 
phenomena  occur  are 


$ 

=  \S 

(a) 

=  ee 

(b) 

-?■ 

J 

=  ae 

(c)  . 

For  the  problem  considered  here,  it  is  sufficient  to  assume 
|i,  e,  a  are  scalars  and  are  constant  with  respect  to  variations 
in  position  and  electromagnetic  fields  (i.e.,  the  medium  is 
linear,  homogeneous,  and  isotropic,  which  is  assumed  to  hold 
throughout  the  remainder) .  The  changes  in  media  will  be  re¬ 
flected  in  the  application  of  boundary  or  interface  conditions . 
Also,  for  the  moment,  it  will  be  assumed  |i,  e,  a  may  be  con¬ 
sidered  time  invariant,  and  the  problem  solved  for  this  "static" 
condition  and  later  correction  is  made  for  the  time  variation. 
Under  the  above  assumptions  and  the  assumption  ( see  below)  that 
the  boundary  or  interface  conditions  are  linear  and  time  in¬ 
variant,  Maxwell's  equations  become  a  linear  and  time  invariant 
operator . 

In  the  case  where  Maxwell's  equations  are  linear  and  time 
invariant,  it  is  convenient  to  use  monochromatic  fields.  Let 
all  the  electromagnetic  fields  have  a  time  variation  of  the 
form  cos ( cot  +  9) ,  where  cd  is  the  radian  frequency  for  all  the 
fields .  The  vector  time  functions  can  then  be  obtained  .from 
the  complex  vectors  by  the  relation 

e(r,t)  =  Re{^(  r,  ioo)  eia>t|  (1.2,5) 

where  r  =  (x,y,z)  —  the  ordered  triple  of  the  rectangular 
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y  _  y 

coordinates  of  position  ,  i  .e . ,  e(r,t)  =  e(x,y,z,t).  Similar 
equations  hold  for  the  other  electromagnetic  fields. 

For  monochromatic  electromagnetic  fields ,  in  a  linear, 
homogeneous,  isotropic,  source-free  region,  the  (complex) 
partial  differential  equations  that  the  electromagnetic 
fields  must  satisfy  become 


curl  t  =  -zif 

(a) 

curl  if  =  vf 

(b) 

(1.2.4) 

div  t  =  0 

div  it  =  0 

(c) 

Z  -  icup. 

Y  =  0  +  inie  . 

(1.2.5) 

Another  set  of  requirements  the  electromagnetic  fields 

must  satisfy  are  the  boundary  or  interface  conditions . 

These  reflect  the  electromagnetic  properties  of  the  different 

media.  It  can  be  shown  from  the  set  of  partial  differential 
3 

equations,  that  at  the  interface  between  two  different  media 

i)  tangential  electric  and  magnetic  fields  are  continuous 


nx[^+( r)  -  E~(r)]  =  0 

ax[g+(£)  _  t-(r)]  =0 

where 

if+(f)  =  lim  ll’(f) 
r->S 
reV 

S  —  surface  ( interface)  bounding  volume  V+ 


(1.2.6) 


3 

^This  approach  that  the  boundary  conditions  are  derivable 
from  the  partial  differential  equations  is  taken  by  H.  Bremmer 
in  his  "Propagation  of  Electromagnetic  Waves,"  in  Handbuch  der 
Physik  Band  XVI  Elektrische .  Felder  und  Wellen  ( Berlin: 
Springer,  195a,  pp.  4^4-425) . 
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n  --  outward  normal  to  S  (pointing  out  of  V+)\ 

Similar ily 

^“(r)  =  lim  ^(r) 
r->S_ 
rev 

ii)  normal  electric  and  magnetic  displacements 
n  •  C$+(£)  -  ^"(r)  ]  =  0 

n  •  C^+(r)  -  il“(r)]  =  0  .  (1.2.7) 


Equatj.ons  (1.2.6)  and  (1.2.7)  are  the  interface  conditions. 
As  stated  above,  the  interface  conditions  are  linear  and 
time  invariant. 

There  are  also  sets  of  boundary  conditions;  they  may 
be  applied  only  under  somewhat  idealized  conditions  (for 
Maxwell's  equations  the  interface  conditions  always  apply) . 

A  convenient  though  somewhat  idealized  model  of  some  real 

5 

materials  is  that  of  a  perfect  electric  conductor.  A  dis¬ 
cussion  of  the  boundary  conditions  at  a  perfect  electric  con- 
ductor  is  given  by  Stratton;  they  are  limiting  cases  of  the 
interface  conditions. 

i)  tangential  electric  field  is  zero. 

n  x  =  0  (1.2.6) 


If 

For  the  definition  of  outward  normal,  see  J.  W.  Gibbs 
(The  Scientific  Papers  of  J.  W.  Gibbs,  Vol.  Two  (New  York: 
Dover,  l%lj ,  p.  JZ)  . 

5 

A  perfect  electric  conductor  is  a  material  in  which  the 
electric  field  is  zero. 

^Stratton,  op.  cit.  pp,  485-484.  In  this  case,  it  is  not 
assumed  that  the  surface  current  is  zero;  however,  this  is  not 
a  physically  realizable  problem. 
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.ii)  normal  magnetic  field  is  zero. 

n  •  lt+  =  0  .  (1.2.9) 

The  above  set  of  boundary  conditions  is  considerably  more 
convenient  than  the  general  interface  conditions  of  equations 
(1.2.6)  and  (1.2.7) . 

Unfortunately,  the  above  sets  of  conditions  (partial 
differential  equations  and  interface  conditions)  may  still 
not  be  sufficient  to  completely  determine  the  electromagnetic 
fields.  A  "boundary"  condition  may  be  necessary  in  the  "far 

7 

field"  if  the  medium  is  unbounded.'  That  is  to  say,  even 
though  in  the  "far  field"  there  may  be  no  change  in  the  medium, 
it  may  still  be  necessary  to  restrict  the  electromagnetic 
fields  by  placing  added  requirements  on  them.  Intuitively, 
if  the  sources  and  virtual  sources  are  bounded  in  extent,  the 
fields  far  from  these  sources  (i.e,,  far  fields)  must  be 
"outward  traveling  waves."  In  fact,  it  can  be  shown  in  the 
far  field  that  the  electromagnetic  waves  are  approximately 
outward  traveling  plane  waves . 

o 

Mathematically,  Sommerf eld’s  radiation  condition0  (for 
two-dimensional  space),  is 

lim  JT  {^P  +  ik}  =  0  (1.2.10) 

where  p  =  ,/x^  +  y^ 


Ibid,  pp.  485-486. 


°Sommerfeld,  Arnold,  Partial  Differential  Equations  in 
Physics  (New  York:  Academic  Press,  1964) ,  pp.  189-19^ • 

Jones,  D.  S.,  The  Theory  of  Electromagnetism, (Oxford: 
Pergamon  Press,  1964 ) ,  p.  93 • 
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and  (1.2.10)  holds  uniformly  in  direction.  Equation  (1.2.10) 
holds  for  a  solution  to  the  scalar  wave  equation  and  there¬ 
fore  holds  for  each  rectangular  component  of  the  vector  elec¬ 
tromagnetic  fields.  The  Sommerield  radiation  condition  suf¬ 
fices  when  the  sources  (both  real  and  virtual)  are  of  finite 
extent.  Unfortunately,  in  the  problem  considered  here,  both 
the  sources  and  virtual  sources  are  not  of  finite  extent.  In 
this  somewhat  more  complex  problem,  the  equations  are  obtained 
as  limiting  cases  (the  computations  are  given  below  in  the  sec¬ 
tion  1.4) .  Toward  this  end,  a  second,  somewhat  more  general 
"radiation  condition"  is  used.  It  is  assumed  that  the  medium 
has  some  losses  (a  >  0;  k  =  k.-  ^2?  k^,  k2  ^  0)  ;  later  the 

limit  as  a  tends  to  zero  is  taken  to  obtain  the  solution  for 

a 

a  lossless  medium.  The  last  requirement  on  the  fields  is 
that  they  must  satisfy  the  source  condition.  This  require¬ 
ment  is  considered  in  detail  in  the  section  on  integral  equa¬ 
tions.  The  problem  of  satisfying  the  partial  differential 
equations,  interface  conditions,  source  conditions,  and  pos~. 
sibly  some  form  of  the  radiation  condition  (i.e„,  finding 
fields  that  satisfy  Maxwell's  equations)  is  a  mathematically 
well-set  problem,  the  solution  of  which  is  unique. 

As  the  use  of  monochromatic  fields  simplified  the  elec¬ 
tromagnetic  boundary  value  problem,  so  the  use  of  vector  po¬ 
tentials  in  many  cases  further  simplifies  the  boundary  value 
problem.  The  great  power  of  the  vector  potential  method 

- 5 - - - 

^This  view  is  taken  by  Baker  and  Copson,  (Baker,  B.  B., 
and  E.  T.  Copson,  Mathematical  Theory  of  Huygen's  Principle, 
Second  Edition  (Oxford:  clarendon  Press,  1953),  P*  154. 
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(along  with  dividing  the  electromagnetic  fields  into  TS  and  TM 
modes)  ±s  given  by  Schelkunof f' ^  although  the  method  is  con¬ 
siderably  older. ^  Mathematically,  the  starting  point  for  ob¬ 
taining  the  vector  potentials  is  the  equations  for  linear, 
homogeneous,  isotropic  and  time  invariant  media  (which  is  as¬ 
sumed  below) . 

div  ^  =  0  ->  ^  =  curl  f  div  =  0  -*•  it  =  curl  ^  (1.2.11) 

12 

The  general  realtion,  xn  a  source-free  region,  is: 

=  -  curl  7?  +  i  curl  curl  it 

§  =  curl  ^  ^  curl  curl  ?  (1.2.12) 

However,  an  additional  problem  involved  in  the  use  of  the  vec¬ 
tor  potentials  is  the  interface  condition.  In  some  cases  the 
simplicity  gained  by  having  to  consider  only  a  single  scalar 
component  of  the  vector  potential  may  be  lost  when  the  inter¬ 
face  conditions  are  applied.  In  general,  interface  conditions 
are  difficult  to  apply  to  potential  functions  (this  can  be 
done  only  in  certain  coordinate  systems  and  only  if  the  inter¬ 
face  is  a  surface  generated  by  one  of  the  variables  equal  to 
a  constant) ?  even  in  the  many  problems  where  the  interface  can 

- - 

Schelkunof f,  S.  A.,  Electromagnetic  Waves  (New  York, 

D.  Van  Nostrand,  1943),  pp.  l27~l29. 

■^'Browich,  T.  S.,  "Electromagnetic  Waves,"  Phil.  Mag., 

Vol .  38,  1919,  PP.  144-164. 

12 

Harrington,  Op.  cit . ,  p.  129 .  .  . 


24 


be  satisfied,  ~ew  components  of  the  vector  potential  are 
needed  to  satisfy  these  conditions.1-^ 

In  the  two-dimensional  problem,  the  vector  potentials 
are  used  to  compute  the  fields  on  the  sea  surface  which,  for 
the  frequency  considered,  is  assumed  to  be  a  perfect  elec¬ 
tric  conductor.  There  is  no  cross -coupling  between  modes  (TE 
and  TM)  and  a  single  scalar  f  iel^fr  may  be  used  to  generate  the 
complete  electromagnetic  fields  of  the  mode.  This  can  be  seen 
from  the  equations  given  below  which  may  be  obtained  from 
equation  (1.2.12). 


14 


TM  Case  (TM  to  z) 

F  =  I  S2  CO 
X  Y  r)X-)  Z  ^ 


Hx  - 


H  1  O 

'  Y  5y35f  * 


E  = 


y  (~^Z  +  ’K2)  <P 
Y 


Hy  Sx 


Hz  =0 


(1.2.13) 


where 


(v2  +  k2)  cp  =  0 


k  =  -  ,/YZ 

in  particular  k  =  — 
c  o  c 


^Sommerfeld,  op.  cit . ,  pp.  246-65. 

14 

Harrington,  op.  cit.,  p.  129.  In  the  simplest  terms 

TM  to  z  means  H  =  UT  TE  to  z  means  E  =  0 . 

z  ’  z 


TE  case  (TE  to  z) 


E  =  - 
x 


E  =  - 


E  =  0 


3cp 

dcp 

3x 


where 


H .  =  i 


x 


H  =  £ 


H  = 


1 

Z 


3x6z  ^ 

d2  » 

oyoz  ^ 

{-^2  +  *2)  cp 
BzZ 


{?  .2.14) 


(v2  +  k2)  cp  =  0, 


The  solutions  to  Maxwell's  equation  may  then  be  generated  by 
considering  solutions  of: 

( v2  +  k2)  cp  =  0 


with  the  appropriate  boundary  and  source  conditions. 


1.3  Physical  Model  for  the  Two-Dimensional  Problem 

In  this  section  a  more  detailed  description  of  the  two- 
dimensional  problem  is  given.  The  problem  considered  may  be 
stated  as  determining  the  electromagnetic  fields  in  the  sea 
when  a  VLF  or  ELF  electromagnetic  plane  wave  is  incident  on 
the  sea  surface.  We  shall  restrict  ourselves  to  the  computa¬ 
tion  of  the  electromagnetic  fields  in  the  sea  relatively  near 
the  sea  surface  (within  about  twenty-five  meters  or  seventy- 
five  feet)  along  with  a  discussion  of  the  asymptotic  behavior 
of  the  fields  far  from  the  sea  surface.  The  electromagnetic 
fields  near  the  interface  are  the  fields  of  greatest  interest 
in  communication  systems  involving  submerged  antennas. 
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The  first  consideration  in  setting  up  the  model  for  the 
above  electromagnetic  boundary  value  problem  is  to  define  or 
at  least  characterize  the  term  "sea  surface"  mathematically. 

A  discussion  of  the  sea  surface  is  given  in  Appendix  A. 

As  the  mathematical  model  of  the  electromagnetic  boun¬ 
dary  value  problem  with  the  rough  sea  surface  it  will  be  as¬ 
sumed  here  that: 

1.  The  "basically  spherical"  earth  may  be  replaced  by  a 
"basically  flat"  earth.  This  is  an  often  made  approximation  in 
"low  frequency"  propagation  problems  and  it  reduces  the  compli¬ 
cated  spherical  geometry  to  a  plane  surface  geometry.  In  the 
problem  considered,  only  the  local  fields  are  of  interest, 

and  the  above  approximation  is  quite  good;1^  however,  it  does 
introduce  some  mathematical  complexity  as  now  the  virtual 
sources  which  represent  the  scattering  effect  of  the  sea  sur¬ 
face  are  no  longer  bounded  in  extent. 

2 .  We  shall  take  as  a  deterministic  mathematical  des¬ 
cription  of  the  sea  surface 

y(x,z),  =  A  cos  k  x  =  ?(x)  (l.J.l) 

b 


This  is  a  surface  that  varies  sinusoidally  in  the  x-direction 
and  is  constant  in  the  z-direction  (see  Figure  l) .  As  dis¬ 
cussed  in  Appendix  A,  this  is  a  particular  realization  of  a 
random  process.  Consideration  of  the  statistical  signifi¬ 
cances  of  this  fact  is  postponed  until  Chapter  J.  The 

- p- - 

^Anderson,  loc .  cit . 
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relative  magnitudes  of  the  physical  parameters  is  of  great  im¬ 
portance  in  the  solution  of  the  integral  equations  governing 
the  fields  on  the  air  side  of  the  sea  surface.  The  relation¬ 
ships  are1^ 

k  A  £  1/7  usually  k  A  <<  1 

o  o 

L  <<  1  and  k  <<  1 
s 

A  >>  L  and  kQA  <<<  1. 

Where  L  is  wave  length  of  sea  wrve, 

2-n- 

k  =  — —  —  wave  number  of  the  sea  surface 

s  L 

2-rr 

k  =  — i.  —  wave  number  of  the  radio  wave  in  the  air 
o  A 

17 

A  is  wave  length  of  the  radio  wave  in  air.  '  We  note  that  the 
assumed  shape  of  the  wave  is  time  invariant  (i.e.,  the  sea  sur¬ 
face  is  represented  as  a  standing  wave) .  In  practice,  sea 
waves  are  actually  traveling  waves  that  even  change  their  shape 
with  time.  However,  the  velocity  of  the  electromagnetic  waves 
in  air  is  so  much  greater  than  the  sea  wave  velocity  that  as 
far  as  electromagnetic  fields  are  concerned,  the  sea  waves  can 
be  considered  stationary.  However,  the  velocity  of  electro- 

- - T7T - 

uSee  Appendix  A  for  these  relations. 

"^For  the  problem  considered  here,  typical  values  of  the 
parameters  are:  1Q8  *  ^  a  1QH  meters 

lCr  a  l  a  10  meters 
25  ^  A  ^  0  meters 
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magnetic  radiation  is  considerably  less  in  sea  than  in  air, 
and  as  a  result  the  effect  of  the  sea  wave  velocity  is  more 
pronounced  in  the  sea  than  in  the  air.  In  Section  2.9,  the 
effect  of  the  sea  wave  velocity  on  electromagnetic  fields  in 
the  sea  is  considered. 

3.  The  source  of  the  electromagnetic  fields  is  a  plane 

wave  incident  on  the  sea  surface .  One  of  the  reasons  for  such 

an  assumption  is  that  in  the  theory  of  ELF  and  VLF  radio  wave 

propagation,  plane  waves  play  a  major  role.  In  both  ray  theory 

for  VLF  and  mode  theory  for  ELF,  a  general  feature  is  the  local 

18 

fields  are  represented  approximately  as  a  sum  of  plane  waves . 

The  incident  plane  wave  propagating  in  an  arbitrary  direc¬ 
tion  above  the  rough  sea  surface  may  be  resolved  into  two  com¬ 
ponent  plane  waves;  one  propagating  in  the  direction  in  which 
the  surface  varies  and  the  other  propagating  in  the  direction 
in  which  the  surface  is  constant.  A  second  resolution  of  the 
problem  is  made  on  the  basis  of  the  polarization  of  the  incident 
field  (see  Figure  l) .  To  solve  completely  the  plane  wave  two- 
dimensional  surface  problem  at  a  fixed  frequency  and  for  a 
fixed  or  given  surface,  four  boundary  value  problems  must  be 
solved  with  arbitrary  angles  of  incident  for  the  source  waves. 

TU  . 

uThe  sky  wave  in  ray  theory  is  represented  as  a  sum  of 
rays  which  are  approximately  plane  waves  [see  H.  Bremmer, 
Terrestrial  Radio  Waves  (New  York:  Elsevier,  19^9),  p.  89]. 

In  mode  theory,  the  modes  are  TE  or  TM  plane  waves .  The 
ground  wave  is  also  usually  represented  as  a  plane  wave. 


Vertical 

Polarization 


c  z 

TM  to  z 


Figure  1.1.  Diagram 
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It  is  to  be  noted  that  the  use  of  plane  waves  as  the 
sources  introduces  some  mathematical  complexity  over  the  simple 
line  source.  The  source  condition  for  the  plane  wave  source  is 
computed  as  a  limit  of  the  source  condition  of  a  line  source. 

4.  To  compute  the  fields  on  the  air  side  of  the  sea  sur¬ 
face,  the  sea  is  assumed  to  be  a  perfect  electric  conductor. 

The  sea  is  a  very  good  conductor  at  the  frequencies  considered 
.  19 

here.  The  assumption  that  the  surface  impedance  of  the  sea 
is  negligible  does  not  greatly  effect  the  tangential  magnetic 
fields  which  are  used  to  compute  the  fields  in  the  sea.  With¬ 
out  this  assumption  the  use  of  vector  potentials  would  be 
greatly  limited  and  a  somewhat  more  complex  problem  would  re¬ 
sult.  The  justification  of  letting  the  surface  impedance  be 
negligible  is  considered  later. 

After  computing  the  fields  on  the  air  side  of  the  sea  sur¬ 
face,  the  fields  in  the  sea  are  considered  on  the  basis  that 
the  tangential  magnetic  fields  are  continuous.  The  electro¬ 
magnetic  properties  of  the  sea  are  assumed  to  be: 

IS 

“H  =  577  — *  "impedance"  of  air. 

\  =  X  XjSe - "impedance"  of  sea. 

«  X^F  =  (l+j)  (  .17)  f  =  50  kc/s  a  =  4  h  =  nQ 

as  T)  ~  */T 
c 

hcl  «<  % 


f  <  50  kc/s 
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a  =  4  I  mhos/meterl 
— 


>-7j" 


1 


|i  =  47rlO_ '^henry/meter  j 


[^f  arad/meter  j 


€  = 


_  9-10 


-9 


4ir 


Using  assumptions  one  through  four,  the  physical  problem 
of  the  electromagnetic  wave-rough  sea  surface  interaction  may 
be  given  a  rigorous  mathematical  formulation. 


1,4  Formulation  of  the  Electromagnetic  Boundary  Value 
Problem  in  Terms  of  Integral  Equations 

In  the  past  few  decades  there  have  appeared  numerous  works 

in  which  important  results  relating  to  the  solutions  of  boundary 

value  problems  were  obtained  by  techniques  involving  the  use  of 

20 

integral  equations.  Integral  equations  have  been  used  exten¬ 
sively  in  the  theoretical  analysis  of  boundary  value  problems 

/  .  N  21 
(e.g.,  proofs  of  theorems  of  existence,  uniqueness,  etc,). 

However,  integral  equation  techniques  have  only  more,  recently 

been  applied  with  success  to  the  practical  solution  of 

Mikhlin,  S,  G.,  Linear  Integral  Equations  (Delhi:  Hin¬ 
dustan  Publishing  Corp.,  I960),  pp.  145-213. 

,  Integral  Equations  (New  York:  Pergamon 
Press,  T3577 7  pp.  137-33 3. 

Courant,  R.,  Methods  of  Mathematical  Physics,  Vol.  II, 
Partial  Differential  Equations  (New  York:  Interscience,  1^62), 
pp.  240-320. 

21 

See  Mikhlin  and  Courant  references  directly  above 
(footnote  21) . 
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boundary  value  problems  (e.g.,  the  work  on  static  elasticity). 
These  recent  successes  in  the  use  of  integral  equation  tech¬ 
niques  have  led  to  their  application  to  problems  which  cannot 
be  solved  readily  by  other  methods . 

The  same  story  basically  holds  true  in  electromagnetic 
theory.  For  electromagnetic  fields,  the  usual  integral  formula¬ 
tion  is  in  terms  of  the  well-known  formula  of  Stratton-Chu  (in¬ 
tegral  representation  of  the  solution  to  Maxwell's  equations  in 
terms  of  the  boundary  values) .  ^  However,  only  recently  has 

the  Stratton-Chu  formula  and  the  integral  formulation  of  the 

24 

reciprocity  theorem  been  used  for  the  practical  solution  of 
electromagnetic  boundary  value  problems . 

One  of  the  major  advantages  of  the  integral  equation  ap¬ 
proach  is  that  the  boundary  and/or  interface  conditions  along 
with  the  source  conditions  are  immediately  considered  at  the  be¬ 
ginning  of  the  problem.  Also,  because  the  integral  equations 
represent  actual  processes  Involved  in  the  physical  problem, 

- — 77 - - 

Maskbelishvilll,  N.  I.,  Singular  Integral  Equations 
(Netherlands,  Groningen:  Erven  ‘t,  ^oordhoff ,  1953T7 

,  Some  Basic  Problems  in  the 
Mathematical  Theory  o£  Ellas  tic  ity  (Netherlands  .Groningen : 

Erven’ '’No'orHKoH; 

^Stratton,  op.  cit „ ,  pp.  464-468. 

24  ~ —  - 

Godziwski,  Z.,  "The  Surface  Impedance  Concept  and  the 
Structure  of  Radio  Waves  ever  Real  Earth  ( IEE,  1961) . 

Feinberg,  E.  L.,  "Propagation  of  Radio  Waves  Along  an 
Inhomogeneous  Surface,"  Nuevo  Cinento,  Series  10,  Vol.  11, 
no.  1  Suppl • ,  1959,  P.  6ET  ~ 

Harrington,  op.  cit.,  pp.  116-120,  317-338,  340-365. 
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the  integral  equations  permit  direct  interpretation  in  terms 

2^ 

of  virtual  sources.  The  use  of  the  concept  of  virtual 

sources  is  very  useful  in  the  mathematical  formulation  of  the 

If) 

electromagnetic  boundary  value  problem.  Further ,  by  inter¬ 

preting  the  integral  equations  in  terms  of  actual  physical  pro 
cesses,  the  problem  of  obtaining  valid  approximations  is  con¬ 
siderably  reduced.  Generally  it.  is  easier  to  see  valid  ap¬ 
proximations  in  physical  terms  rather  than  in  mathematical 
terms . 

In  a  very  practical  vein,  the  integral  equation  method  in 
volves  a  reduction  of  "dimensionality,"  that  is,  a  two-dimen¬ 
sional  boundary  value  problem  can  be  represented  as  an  integ¬ 
ral  over  a  one -dimensional  space?  similarly,  a  three-dimen¬ 
sional  boundary  value  problem  leads  to  an  integral  in  two  di¬ 
mensions.  From  the  area  of  pure  mathematical  analysis,  though 
it  has  great  practical  implications;,  integral  operators  are 
considerably  easier  to  handle  than  differential  operators,  par 
ticularly  when  approximations  are  necessary. 

The  above  reasons  of  course  are  not  necessarily  suffi¬ 
cient  to  cause  all  electromagnetic  boundary  value  problems 

- jc - 

JStratton,  op.  cit. . ,  p.  467. 

Baker,  B.  B.  and  E.  T.  Cops on,  op.  cit . ,  p.  114. 

26 

^Waterman,  p.  c.,  "Scattering  of  Electromagnetic  Waves 
by  Conducting  Surfaces,"  Wilmington,  Mass.:  Research  and  Ad¬ 
vanced  Development  Division  Avco  corporation,  Dec.  1962,  an 
unpublished  report. 
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to  be  formulated  and  solved  in  terms  of  integral  equations. 

A  major  disadvantage  of  the  integral  formulation  is  that  often 
to  compute  the  fields  at  any  point  an  integration  is  necessary. 
If  only  the  far  fields  are  ne^led,  asymptotic  expansion  of  the 
integral  may  be  used  to  obtain  the  fields  in  closed  form  which 
does  exhibit  directly  the  variation  of  the  field  with  respect 
to  position.  In  the  near  field,  an  integration  of  an  integral 
with  a  very  complex  kernel  is  needed  for  each  point  at  which 
the  field  is  to  be  computed;  this  is  very  laborious  even  with 
high  speed  computers . 

For  the  rough  surface  problem  considered,  the  integral 
methods  of  formulation  are  relatively  convenient.  The  classi¬ 
cal  methods  applied  to  differential  operators  (such  as  sepa¬ 
ration  of  variables)  is  extremely  inconvenient,  particularly 
in  that  the  unknowns  in  this  case  are  not  susceptible  to  a 
physical  interpretation  that  permits  approximation.  However, 
some  general  qualitative  results  can  be  obtained  by  arguments 
based  c-n  the  separation  of  variables  . 

The  usual  starting  point  for  the  study  of  electromagnetic 

boundary  value  problems  by  use  of  integral  equations  is 

27 

Green's  theorem.  '  In  two-dimensional  space,  Green's  theorem 
takes  the  form 


j  {cp(p)  V2!|((p)  -  t(p)  V2cp(p)}dS  =  {cp(  p)tn(p)  -  ^(pHCp)}^ 

(i.4.i) 


- 77 - 

’Kaplan,  Wilfred,  Advanced  Calculus  (Reading,  Mass.; 
Addison-Wesley,  1953) ,  P.  775. 


35 


where 


p  =  (x,y)  —  ordered  pair  of  rectangular  position 
coordinates  in  two-dimensional  space. 

cpn(  o)  =  n  •  grad  cp(  p)  --  normal  derivative. 

cp,  i|r  are  two  arbitrary  scalar  functions  of  position  with  sane 

requirements  cn  their  derivatives. 

Surface  S  is  bounded  by  a  closed  curve  C. 

The  function  i|r  will  be  used  as  an  auxiliary  function  and 

restrictions  will  be  placed  on  cp  such  that  it  will  be  the 

2g 

solution  of  a  scalar  boundary  value  problem.  As  a  first 


step  in  the  above  process,  let 

(V2  +  k2)  *(p,  51)  =  -  l(t-t') 

(v2  +  k2)  =p(p,  p")  =  -  c^fp-p”)  (1.4.2)29 

where 

p’  =  xx+yy--  the  position  vector?  i.e.,  cp  and  i|r 


satisfy  the  Helmholtz  equation. 

Placing  restrictions  (1.4.2)  into  (1.4.1) 
<p(p',P")  P'eS 


0 


( P" , P  1 ) 
0 


+ 


p '/(sue) 

_t(p,p')  <Pn(P,P")  -  <P(p,PM)  fn(p,p*)]d£ 


c 


o"  es  [ 

P"/(suc)J 

Ci.4.5) 


- 23 - 

The  solution  of  Maxwell’s  eq.  ^tions  (i.e,,  the 
electromagnetic  fields)  is  obtained  from  the  solution  of 
the  scalar  boundary  value  problem, 

A^The  6(£-p')  symbol  represents  a  source  term.  The  use 
of  this  notation  is  given  in  Friedman,  Principles  and  Techniques 
of  Applied  Mathematics  (New  York:  JohnTTiley,  I^)b) ,  pp.  154- 
THfe ,  where  tTie  mathematical  properties  of  the  "  6  -  function" 
are  given. 
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where 


SUC  denotes  the  union  of  the  sets  S  and  C  (i.e.,,  the  points 
belonai nq  to  either  S  or  C  or  both) ; 

Des  denotes  0  is  an  element  of  S  (i.e.,  p  is  a  point  in  S. 


Now  to  form  an  integral  equation,  let  p •  ->  C.  The  is  a 

Jc 

"singular"  or  "discontinuous"  integral.  In  Appendix  C  it 

is  shown  that 

p  _  p 

iim  |  ♦n(p,P,)cp(p,p")d4  =  !  'Jrn(o,p,)cp(p,p")d't  -^p(p',p") 

p'Sc  P'9c  (1.4.4) 

p 1  s  s  p "  s  s 

p"  eS 


(1.4.4) 


p  p 

lim  |  cPn(o}P,')’l'(o,D,)(3t  =  1  <Pn(  o,  r"  )  |(  o,  p*  )dt 

p""Sc  c"ec 


p"  ->c 
p"  es 
P'  eS 


c"  ec 
p  *  £  S 


£y(5-,5") 

(1.4.5) 


Applying  (1.4.4)  to  (1.4. 5) 

cp(o',P")  =  2lLy(5-,p')  +  K(5>5")*(ff>o') 


-  tn(  p,  o ' )  cd(  o,  p" )  jd-c  (1.4.6) 


P  ,  P  1  £  C  p"  £  S 


Equation  (1.4.6)  is  the  basic  integral  equation  used  below 
in  the  mathematical  formulation  of  the  electromagnetic  boundary 
value  problem.  It  holds  only  for  closed  curve  C. 

Equation  (1.4.6)  may  be  rewritten  with  t(p,p')  -  G(p,P') 
(this  is  a  restriction  of  the  general  equation),  where  G(o,p') 


is  the  "free-space  Green's  function." 

=  i  h(2)(!<I  =-?•!  ) 


(1.4.7) 


:'.(o',o“)  =  2  s(fiG(c”,H’)  +  (  {cpn(  o,  p")g(  0,  p  ') 

-  (1.4. 8) 
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cp  given  by  (1.4.8)  satisfied  the  Helmholtz  equation  (1.4.2); 

a  further  requirement  of  the  electromagnetic  fields  is  that 

they  must  satisfy  a  source  condition.  In  (1.4,8)  the  source 

condition  is  represented  by  the  term  with  as  a  factor. 

is  the  strength  of  the  source,  eCpG(p*,p")  represents  the 

incident  field,  i.e.,  the  field  at  p'  due  to  source  at  p" 

3  ’  cp 

if  no  scatterer  were  present  jcp^(  p  • ,  p"  )  =  eCpG(p",p')}  •  The 
field  represented  by  the  integral  then  has  the  physical  in¬ 
terpretation  as  the  effect  of  the  scatterer.  The  scattering 
term  is  immediately  interpreted  in  terms  of  virtual  sources 
on  the  scattering  sur: ace  ( or  curve) . 

For  the  boundary  value  problem  considered  here  the 
scattering  surface  (or  curve)  is  of  unbounded  extent.  The 
closed  curve  C  representing  the  scatterer  is  composed  of  two 
parts; 

i)  C  —  curve  (surface)  which  is  the  mathematical 
representation  of  the  sea  surface, 

ii)  Coo  --  the  infinite  semicircle  (hemisphere)  which 
closes  C . 

For  purposes  of  analysis,  the  field  cp  under  consideration 
is  divided  into  three  parts . 

cp  ==  cp  +  cp  +  cp  •  (1.4.9) 

X  XT  o 

where 

cp .  —  the  incident  field  {a  plane  wave  in  the  problem 
consic" ared  here;  however,  the  source  condition  for  the  plane 
wave  will  be  obtained  as  a  limit  of  source  conditions  for 
finite  sources,  i.e.,  the  incident  field  considered  just 


3>B 

below  is  a  cylindrical  wave  resulting  from  a  finite  line 
source) . 

--  the  reflected  field  (a  plane  wave  in  the  problem 

considered  here.  It  is  the  reflected  wave  if  the  scatterer 

were  a  plane,  i.e.,  the  specularly  reflected  field  basically). 

cpg  --  the  diffused  scattered  field  (which  represents 

the  effect  of  the  roughness  of  the  scattering  surface) .  It 

satisfies  the  Sommerfeld  radiation  condition. 

The  integral  equation  (1.4.8)  will  now  be  applied  to 

each  term  of  the  sum  on  the  right-hand  side  of  (1.4.9) • 

For  the  incident  field  cp. 

1 

cPi(p,,P")  =  2[cpi(p',p")  +  J  {ep±n(  p,p")G(p,p') 

C 

-  Gn(p,p')cpi(p,p")}d^]  (a) 

(1.4.10) 

=  2[cp;.(p,p")  +  j  {^in(  P,P")G(  P.P') 

-  Gn(  P*)^  P,P”)}<3^]  (b) 

as 

f  [tPin(P,P")G(p,P')  -  Gn(p,P')q>  (p,p”)}di]  =  0  (1.4.11) 

dc« 

by  Sommerfeld' s  radiation  condition. 

(1.4.11)  need  only  hold  if  the  sources  are  of  finite 
extent.  To  obtain  the  equivalent  integrals  (i.e.,  the  integral 
representing  the  incident  field)  to  (1.4.10  and  1.4.11)  when 
the  sources  are  unbounded  in  extent  (i.e.,  the  source  or 
incident  field  is  a  plane  wave) ,  a  limiting  process  is  neces¬ 


sary. 
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To  start  the  limiting  process,  let 
cPi(  P  ’ ,  P"  )  =  elk?H^  {k(  Pp  '  -  £.?!)} 

P"  =  1<i? 


(1.4.12) 


where  £  is  to  be  the  unit  vector  in  the  direction  of  propa¬ 
gation  of  the  incident  plane  wave,  and  ?  is  a  parameter  used 
in  the  limiting  process . 

Equation  (1.4.12)  represents  the  field  caused  by  a  line  source 


at  p”  of  strength 


e  b.  As  intuitively  expected,  the 


source  point  distance  from  the  observation  point  must  become 
unbounded  (tend  toward  infinity)  and  the  source  strength  also 
becomes  unbounded  to  obtain  a  plane  wave  at  the  observation 
point .  Both  the  above  occur  as  ?  ->  ® . 

Applying  the  limit  ?  ->  »  to  equation  (1.4.11)  for  the  cp^ 
given  by  equation  (1.4.12) 

0  =  lim  {tp.  (  p,  p")G{k|  p-p'  |)  -  G  (k|p-p' |)cp  (p,c")}at 
?•*»  <-c»  1 


(1.4.13) 


=  ,{-i«<-ne+ik£l-'?G(k|^ll)  -  G  (klp-p-De^V^W 

Ur*  ^  i  n  + 


(1.4.14) 


lim  cp.(pt)  -  e+ik^i,pl 


lim  cp  (  p  » )  =  in 
in 


•  .  fe.ke+J*V*' 


a  plane  wave. 


Equation  (1.4.14)  represents  the  relation 

o  =  f  {^(pjGjp.p-)  -  opin(p)G(p,P')}di 


(1.4.15) 


which  implies 


O=cpi(0’)  +  2 J{cpi(p)Gn(p,p‘)  -  cp.n(p)G(p,p')}^  (1.4.16) 


where  ecu  is  a\  plane  wave 
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Then  (1.4.10)  holds  even  if  the  source  is  unbounded  in  extent. 


i.e.,  cp±  is  a  plane  -wave. 

For  the  reflected  field  cp 

2T 

«Pr(p)  =  2^  {«Prn,'(  P‘)g(p,  P*)  “  G^p^'H^p*)^' 
c 

=  2  J  {ff>rrl.(  P  *)G(  p,  P'}  -  Gn,(p,P,)«Pr(P’)}<«’’ 
c 


,  (  P  *)g(  p,  P  1 )  -  Gn,(p,p')cpr(p')}^’  =0 


(a) 

(b) 

(1.4.17) 

(1.4.18) 


by  second  radiation  condition. 

For  the  scattered  field  cp 

8 

4>S(P)  =  2  J{cpsn,(p')G(p,p')  -  Gnt(P,P')cpg(p')}df  (a) 

0 

=  2 ^  {<Psn,(p')G(p>P’)  -  Gn. ( e. p ’ ) ^>s( p ■  (b) 

c  (1.4.19) 

as 

I„Kn'(P')G(5'5')  ~  '3n .  (  P ,  ?  * )  (  P  * )  }d^.  •  =  0  (1.4.20) 

by  Sosranerf  eld's  radiation  condition. 

Adding  equations  (1.4.16),  (1.4.17),  and  (1*4.19), 

4>r(p)  +<Ps(p)  =9i(p)  +2^  {cpn,(p’)G(p,Pt) 

c 

-  Gn,(p,  p')cp(  p')}dt*  (a) 

cp(p)  =  2{cp±(  p)  +  J  {cpn«(pi)G(p,P')  -  Gn,(f),p,)cp(p,)}dt}. 

C  (b) 

(1.4.21) 

Equation  (1.4.21)  is  the  integral  equation  for  the 
field  cp  such  that 

1.  cp  satisfies  the  Helmholtz  equation  (1.4.2)  with  0. 

2.  cp  satisfies  the  source  condition  (i.e.,  the  incident 
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field  is  a  plane  wave) . 

3.  cp  satisfies  the  necessary  radiations  (as  seen  above). 
To  compute  the  mathematical  formulation  of  the  boundary 
value  problem,  it  remains  only  to  apply  the  necessary  boundary 
or  interface  conditions  required  of  cp  on  the  interface  C. 

While  the  complete  boundary  value  problem  may  be  formulated 
by  setting  up  the  integral  equations  on  each  side  of  the 
interface  and  applying  the  interface  conditions y  it  is  very 
convenient  to  assume  the  lower  half-space  bounded  by  C  is  a 
perfect  conductor.  The  effect  of  this  approximation  is  given 
below.  The  above  assumption  changes  the  interface  condition 
to  a  boundary  condition.  It  is  also  convenient  to  consider 
the  field  divided  into  two  parts: 

1.  TE  to  z. 

2.  TM  to  z. 

The  TM  to  z  part  may  be  generated  by  a  single  component 
of  the  vector  potential  For  the  basic  equations,  see 

•30 

Harrington  or  Section  2.3.  The  boundary  condition  is  cp  =  0. 

This  yields  E,  =0  and  H  _  =  0.  cp  is  related  to  the 
tan  normal  n 

tangential  H  field. 

The  TE  to  z  part  may  be  generated  by  a  single  component 
of  the  vector  potential  The  boundary  condition  is  cp^  =0; 
this  again  implies  Etan  =  HnQmal  =  0. 

The  boundary  value  problem  may  then  be  formulated  in 
terms  of  a  single  scalar. 

The  TM  case,  the  integ  l  equation  is: 

^  Harrington,  op.  cit.,  p.  129. 
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^n.C  P*)G(ko|p'-^' |)dt*  =  cp±(p)  .  (1.4.22) 

C 

The  TE  case,  the  integral  equation  is: 

«P(p)  =  2cp  (p)  +  2  G  .(k  |t-p'  I «P( P • ) a-t • .  (1.4.23) 

1  Xl  o 

(1.4.22)  and  (1.4.23)  are  the  integral  equations  to  be  solved 
in  the  next  chapter. 


2.0  SOLUTION  OF  THE  INTEGRAL  EQUATIONS 


2.1  Introduction. 

In  this  chapter,  integral  equations  (1.4.22)  and 
(1.4.23)  will  be  solved.  From  these  solutions,  the  fields 
on  the  air-side  of  the  sea  surface  can  be  found.  Also 
considered  are  the  integral  equations  in  the  sea.  From 
the  approximate  solution  of  the  integral  equations  in  the 
sea  comes  an  impedance  type  relation  between  cp  and  cp  which 
holds  in  a  restricted  frequency  range  and  for  a  restricted 
class  of  surfaces.  Using  this  impedance  relation,  a  simpli¬ 
fied  integral  representation  in  the  sea  is  given.  Then  a 
brief  comparison  of  the  solutions  obtained  here  with  those 
obtained  by  other  methods  is  presented  (a  more  detailed 
discussion  of  other  methods  used  to  solve  "rough  surface" 
problems  is  given  in  Appendix  B).  The  following  sections 
are  devoted  to  the  question  of  the  validity  of  some  of  the 
assumptions  used  (such  as,  if  the  sea  can  be  considered  a 
perfect  electric  conductor  in  determining  the  fields  in 
the  air).  Lastly  considered  is  the  effect  of  the  motion 
of  tne  sea  surface  on  the  fields  in  the  sea.  Wait's^ 
approximate  solution  is  used  in  this  discussion,  as  the  use 
of  the  more  accurate  'numerical  solution  would  only  compli¬ 
cate  the  discussion  without  adding  any  new  information. 

^"Wait,  loo.  cit. 
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A  recent  review  work  considering  the  practical  solu¬ 
tions  of  integral  equations  gives  as  two  of  the  nine  or  so 
available  methods 

1.  Methods  of  finite  differences  and  sums. 

2.  Method  of  degenerate  kernel. 

Although  numerical  computations  using  finite  differ¬ 
ence  methods  have  recently  been  shown  to  be  a  useful  tool 
in  the  solution  of  integral  equations,  such  computations 
proved  to  be  impractical  for  the  solution  of  integral 
equations  considered  here.  The  main  reasons  for  the 
difficulty  were: 

1.  The  curve  C  is  unbounded.  The  integral  equations 
solved  by  the  finite  difference  method  have  had  a  finite 
interval  for  the  range  of  integration. 

2.  The  kernel  of  the  integral  equations  does  not 
decrease  rapidly  enough  with  distance  from  the  source  point 
to  allow  useful  approximations  (i.e.,  cutting  off  the  range 
of  integration  to  a  small  interval  about  the  source  point). 
The  computations  therefore  could  not  be  made  sufficiently 
detailed  to  yield  a  useful  solution. 


2 

Walther,  A.,  "General  Report  on  the  Numerical  Treat¬ 
ment  of  Integral  and  Integro-Dif ferential  Equations," 
Symposium  on  the  Numerical  Treatment  of  Ordinary  Differen¬ 
tial  Equations,  Integral  and  Integro-Dxf f erential 
Equations  --  Proceedings  of  the  Rome  Symposium,  i960. 
Organized  by  the  P ICC  (Basel,  Berkhauser,  I960),  p.  649 ♦ 
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For  the  solution  of  one  of  the  integral  equations 
(1.4.22),  however,  a  modified  form  of  the  method  of  degen¬ 
erate  kernel  was  found  to  be  the  most  practical.  The 
general  solution,  using  the  degenerate  kernel  method, 

is  to  expand  (and  approximate)  the  kernel  of  the  integral 

•5 

equation  as  follows:"^ 

N 

K(x,§)  M  Y  ui(x)v.(§)  (2.1.1) 

i=l 

v(x)  =  t (x )  +  X  P  K(x.  E  )y  (? )dE  (2.1.2) 

(2.1.2)  is  the  integral  equation  to  be  solved. 

N 

y  (x  )  =  t(x)  +  y  ciui(x)  (2.1.3) 

i=l 

(2.1.3)  is  the  solution  of  integral  equation  (2.1.2)  with 
the  approximation  (2.1.1)  where  c^  is  the  solution  of  a 
system  of  N  linear  equations.  In  the  problem  considered 
below,  this  approximation  is  modified  and  takes  the  form 

K(x,?)  «  K(x-S).  (2.1.4) 

With  the  aid  of  this  approximation,  the  resulting  integral 

equations  may  be  solved  by  classical  methods  (using  Fourier 
.  4 

Series).  The  approximation  is  valid  because  of  the  relative 
values  of  the  physical  parameters  of  the  problem. 

•^Ibid,  p.  654. 

4 

Morse,  Philip  M.  and  Hermen  Feshbach,  Methods  of 
Theoretical  Physics,  McGraw-Hill*  New  York,  1955/  Part  1, 
pp.  950-962. 
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The  solution  of  (1.4.25)  is  approximately  the  unper¬ 
turbed  cp/  as  is  easily  verified. 

The  solution  to  the  integral  equations  in  the  sea  is 
based  on  the  high  rate  of  attenuation  for  electromagnetic 
waves  in  sea  water. 

2.2  TM  Case. 

The  integral  equation  to  be  solved  in  the  TM  case  is 

J_  <pn,(p')  G(k0|?-?'|)<U*  =  cp.(p)  ,  (2.2.1) 

'c 

where:  k^  -  wave  number  of  the  electromagnetic  wave  (in  air), 

G(k|x|)  =  -^  ^  (k  jx  |  )  -  the  free  space  (two- 

dimensional)  Green's  function. 

It  can  be  shown  that  if 

(5  '  )2  «  1  (k  §)  «  1 

'max  v  o  'max 

where:  ?(x)  is  the  equation  of  the  sea  surface, 

?'(x)  is  the  derivative  of  5  (i.e„,  the  slope  of  the 
sea  surface)  with  respect  to  x. 

Then 

^n*  (p  '  )  G(k0l  P~P  '  I  )d^'  w  <Pn>  (p  '  )G£(k0l  P“P  ’  I  )d^'  > 
C  C 

(2.2.2) 

where : 

G (k0 1 o- p  '  I )  =-5Hi2)(kJ(^'-)2  +  (?(*)  -  ?(*'))2) 


Of(H0|p-p'  I)  = -5  )(ko|x-x‘ |  )  . 


112-0 


mi*  6 
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The  argument  leading  to  this  result  is  presented  in  the  dis- 

.  g 

cussion  on  the  three-dimensional  problem. 

The  integral  equation  then  becomes 

9n.(p')  Gf (k0IP“P ' I )d^‘  =  ^(o)  •  (2.2.3) 

*C 

This  equation  may  then  be  solved  by  Fourier  Series. 

To  solve  the  integral  equations  for  cPn<»  let 


<Pn(p)=£ane 


-ikx.x 


03 


n=-a> 


-I 

n=i~« 


ane 


-ink  x 
s 


-ikxix 


(2.2.4) 


’’i(p)  =  I  V 


-ik  x 
m 


m=-« 


-ikxix 


CO 


=  Fbme 


-imk  x 
s 


m— —  ea 


-ik  .x 
xi 


(2.2.5) 


where  k  =  2ir/L 
s 

p  gC  . 

The  validity  of  equations  (2.2.4)  and  (2.2.5)  is  estab¬ 
lished  by  considering  the  motivation  behind  such  expansions. 
“ikxix 

The  factor  e  is  present  because  the  incident  field 

(source)  has  such  a  factor.  The  sum  expresses  the  fact  that 

the  periodic  surface  implies  a  periodic  virtual  source 

-ik  .x 

(besides  the  factor  e  xx  )  which  in  turn  implies  a  periodic 
field.  The  expansions  (2.2.4  and  2.2.5)  are  the  most  general 
mathematical  statement  of  this  periodicity  of  the  fields. 


c 

^See  page  121.  This  argument  is  given  for  the  two- 
dimensional  problem  by  Meecham,  loc.  cit.,  and  Lysanov, 
loc.  cit. 
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It  should  be  noted  that 


J^Gf(k0|p-p' |  )  e 


-i{nlc  +  k  )x ' 
X1  / 


dx  '  = 


-ik  .x 

XI  / 

e  (■ 


V(nv*xi>*-  ^ 
(2.2.6)6 


)■ 


Then  substituting  equations  (2.2.4)  and  (2.2.5)  into  the 
integral  equation  (2.2.5)  and  using  formula  (2.2.6),  the  rela¬ 
tions  between  the  sets  of  a  's  and  b  's  may  be  found  by 

n  n 

equating  the  coefficients  of  like  terms  in  the  series. 


5 


Equation  (2.2.6)  is  obtained  from  the  relations  given 


above . 


CO 


r 

jl 


)e~ikx 


H^^(kx)  cos  kx  dx 
o  '  ' 

1 


00 


J  (k  x)  cos  kx  dx 
ov  o  ' 


0  <  k  <  kQ 
k  =  k 

o 


> 


l  0  kQ  <  k  < 

(Bateman  Manuscript  Project,  Vol«  I,  McGraw-Hill,  New  York, 
1954;  p.  4 3,  equation  No.  1.) 


00 


v 

Jo 


N  (k  x)  cos  kx  dx  = 
o'  o  ' 


0  <  k  <  k 

o 

► 

k  <  k  <  00 
o 

) 


(Bateman  Manuscript  Project,  Vol.  I,  McGraw-Hill,  New  York, 
1954;  p.  47,  equation  No.  28.) 


Therefore,  equation  (2.2.6)  is  obtained 


1 
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The  relation  between  the  a  '  s  and  b*  's  is 

n  n 


an  -  bn  <2/(^s  +  kxi)2-  k2  ) 


(a) 


a  =  2ik  .  b 
o  yi  o 


a  _  =  a  es  +2nk  b 
-n  n  s  n 


cp  =  e 


-ik  ,.x  ik  .  ?(x) 
xi  -  ' 


(e 


yx 


) 


?  (x)  =  A  cos  (k  x) 


bQ  —  1  -  2  (^yi^ )  +  •  •  •  M  1 

bl  =  b-l  =  ’f  *iZ  +  ••• 


A 

-  lkyi2 


b-n  "  bn  -  |tkyiA)n 


ao  -  2ikyi 

a-i  =  ai  ~  iks  kyiA 

an  “  "(V*)"  n  >  1 


0>) 

(=) 


(2.2.7) 


which  may  be  neglected  because 
of  the  factors  (k  JV)  «<  1. 


CO 


Hp^(k  |xj)e 

JLo  °  ° 


-  ikx 


dx 


-2i 


CO 


k  4  k 
o 


k  =  k 
o 


It  is  assumed  that  k.  +  nk  4  k  for  any  integer  n  (note  this 
assumes  k  .  ^  k  for  n  =  0).  This  assumption  is  useful  in 
avoiding  the  question  of  the  convergence  of  the  integrals 
above.  The  physical  phenomena  occurring  when  k  =  k  .  leads 
to  a  "resonance"  in  the  fields.  This  "resonance"  is  a  physical 
occurrence;  however,  in  some  mathematical  formulations  the 
fields  become  infinite,  which  is  not  a  physical  possibility. 
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i 


-  ik  .x 

*n  ~  21kyi(1  +  ksA  cos  V>e  X1  =  2^ni^x^f 

{l  +  k.?(x)}  (2.2.8) 

where  ©  .  (x),.  is  the  normal  derivative  of  qp.  on  the  flat 
ni'  'f  Ti 

surface.  The  perturbation  compared  to  the  flat  interface 
is  ?■  wgA  whir’.:  ;  •  than  about  thirty  per  cent.^ 


2.3  TE  Case. 

In  this  case,  the  integral  equation  becomes 

<p(p)  -  2^(5)  +  2J  <p  (p  '  )Gn(ki P~P  '  |  )dl‘ 


(2.5.1) 


In  the  flat  interface  case,  the  solution  is 

<P(p)  =  2<pt(p)  .  (2.5.2) 

Using  the  static  kernel  (Green’s  function)  and  placing 

✓ 

the  above  approximation  (equation  2. 2*2)  into  t&e  integral, 
the  integral  becomes  aero.  Therefore,  no  correction  is 
necessary.  The  justification  of  the  use  of  the  static  kerne’ 
is  considered  in  Section  (2.8).  In  this  case,  the  solution 
is 

•P(p )  =  29i(p)  .  (2-5-5) 

This  is  equivalent  to  the  condition  Htan  =  2(Htan)i. 


Htan 


7 

‘As  Htan  -•  <pn,  the  perturbation  (maximum  derivation)  in 

is  2k  A  or  about  28$  for  k  A  s  1/7  • 
s  s 
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2.4  Integral  Equations  in  the  Sea. 

For  the  fields  in  the  sea,  the  integral  equation  (1.4.21) 
still  holds  only  with  =  0. 

<p(p)  =  2  '  (  P  '  )  G(kclp"P’i)  **  «P(p  '  )Gn- (kciP"P  '  I  )}d^' 

C 

(2.4.1) 

where  k  =  -  the  complex  wave  numbe.  in  the  sea  for 

co 

electromagnetic  v;aves . 

In  the  VLF  range,  for  usual  sea  water  (o  =  4,  \i  =  ^0), 

|k  |  £  „5.  This  implies  G  and  G  tend  rapidly  to  zero  for 
|p_p'  I  >0.  In  the.  order  of  10  meters  G  and  Gn  become 

p 

negligible.  If  <p  and  q>n  are  relatively  constant  on  the 
sea  surface  (i.e.,  <y  and  <pn  are  relatively  constant  for 
distances  of  the  order  of  ten  meters  or  so  (J>b)  along  the 
sea  surface)  and  the  sea  surface  is  relatively  constant  over 
these  distances;  then: 


^To  estimate  the  accuracy  of  neglecting  G  and  Gn  after 

a  distance  of  about  *  /’2  from  the  "source  point"  the  asymp- 

c 

totic  approximations  to  G  and  Gn  will  be  used  (the  errors 
in  the  asymptotic  approximations  are  sufficiently  small  to 
be  neglected  here). 


~ikcp 

e 


G  ~ 
n 


- ..  ik  n 
i  m  (-ik  )  e  c 
c 


n 


00 
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Jjpn. (p ' ) G  (kc i p-p 1 1 )d^'  »  ^n(p ) 

c  c 


J^cp (p  • )  Gn,  (kc|?-p'|)di*  «<p(p)*5^“r 

c  c 


(2  o4»2 ) 


(2-40) 


Ac/2 


(k  x)  dx  I  ^  - - —  P  s(x)dx  =  — ~ — 

c  1  2|*d  ’U  2|kc! 


(.0062) 


2  /  \  1  -"XT  /2  , 

§  (x )  = - e  -  normal 


distribution 


CO 

j  G(k  x)dx 


1  V2 


£  c. 


o. 


G(k^x)dx 


c1  =  .003 


00 


Uo  GJ\x)ay.  „  ^ 


•Jo 


an(v)dx 


«  1 


r 

V2 


Gn(kcx)dx  *  J 


JEr 


?i(x)dx  £  .006  «  1 


G  (k  x)dx 
nx  c  ' 


£  C, 


(kcx  )dx 


•Ol|kc|R. 


For  many  cases  C2  may  be  the  order  of  1  or  greater;  this 
occurs  as  R  -*  «.  While  the  important  contribution  to  the 
integral  no  longer  occurs  near  the  "source  point"  the  result 
is  the  same;  that  is 


cp(p')  Gn- (kc|p-p  '  I  )  <U' 

—  CO 


«  1 


which  is  all  that  is  required  for  the  impedance  boundary  con¬ 
dition. 
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where  R  -  radius  of  the  curve.'' 

For  Ik  |R  »  1  which  holds  in  the  VLF  ranqe  and  for  the 
i  c  1 

sea  states  considered: 


*  "  TT  ^n 
c 


(2.4.4) 


This  is  an  impedance  type  relation  between  <p  and  <pn  on  the 
sea  surface. 


9  r>  ® 

To  estimate  <p  ( p  ' )  G  ,  (k  !p-p'|)  d£',  again  <p  ( p  '  ) 
tJ-oo  n  c 

is  assumed  constant  over  the  meaningful  range  of  integration. 

For  small  |p-p'|  — L  «  l)  the  static  Green's  function 

c 

may  be  used. 


Gn.(\p)  «  n'  .grad  *n  p}  *  ^  j  pn 


Pn  = 


f-S’x  + 


l  XX.  + 


x*-  + 


then  for  p  »  0 


Gn '  (k*P  )  w  77  - *“5  -x  ±.l. - 

S&VTP  (x2  *  §2) 

The  Taylor  series  expansions  of  §(x)  may  be  used  for  x  small 


Gn'(kc0>  "Tir 


rg  "X2 


X2(l+(g')2) 


1  1 
77  7r 


(±+(e  '  )2  )3/2 

R  =  ^ — L — L -  -  radius  of  curvature  of  the  curve 

5" 

(surface). 


'«V  - 


\/ 


* 
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As  the  frequency  decreases  into  the  lever  ELF  range, 
|kc|  decreases  and  for  some  sea  surfaces  the  relation  |k  | 

R  »  1  no  longer  holds.  The  impedance  relationship  between 
cp  and  (p^,  however,  may  still  be  an  excellent  approximation. 

This  relationship  (2.4.4)  is  a  Leontovich  type  of 
boundary  condition, ^  and  may  be  expressed  in  terms  of  the 


Using  the  above  result  and  the  assumption  ?"  is  relatively 
constant  over  the  meaningful  range  of  integration  (i.e., 

L  »  A  )  for  the  problem  considered  here: 

f  <P(p  ■ )  Gn,  (kc|p-p|  )<U  =  <p(o>g  2j  H > (kop )kcpn  d l 

«  'pCp^X  J  x  H^OyO  ^ 

«9(p)^^77 

Kc 

w  <p<5)  si 

C 

So  that  if  5"  and  cp  are  relative  constants  over  a  distance 
of  kc/2  or  so  along  the  surface,  then: 


->  <» 

<P(p')  Gn, 

—  CO 


(kc|?-p'|)dr  <p(p)  -2TR 

c 


10 


Leontovich, 


loc.  cit. 


Brekhovskikh,  Leonid  M.,  Waves  in  Layered  Media,  trans¬ 
lated  by  David  Liberman,  Academic  Press,  New  York,  i960,  pp. 
14-15. 


Skin  Depth  (meters) 


Figure  2.1.  Skin  Depth  of  the  Sea 
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fields  as: 

Etan  "  1c  “tan  •  <2'4'5) 

This  is  true  for  the  fields  in  the  sea  *  however,  as  the 
tangential  fields  are  continuous  at  an  .erf ace  (2.4.5)/ 

it  also  holds  for  the  fields  in  the  air. 

2.5  Simplified  Integral  Representation  in  the  Sea. 

The  usual  integral  representation  of  a  wave  function 
(with  qK  =  0)  is 

9(p‘)  “  J  {G(^c  i  P“P  '  !  )  <Pn(p)  “  4  (p  )  G(kc  j  p-p  '  i  )jd6 

0 

(2.5-1) 

which  includes  terms  representing  virtual  sources  of  two 
types : 

<PX(5')  -  J  {o(*0Ip-p'I)  <Pn(p)  }  d-?  represents  the  fields 


of  a  "single  layer"  source. 


(2.5.2; 


5>,(c')  =  '  G  (k  |p-p'|)  ®(p)}d£  represents  the  fields 

4  #  n  o 


of  a  "double  layer"  source  (dipole). 


(2.50) 


In  electromagnetic  problems,  these  types  of  sources  are 
the  electric  current  and  the  magnetic  current  sources.11 
It  would  be  of  interest  if  a  simplified  integral  representa¬ 
tion  of  the  wave  function  could  be  found,  particularly  if  the 
fields  were  to  be  evaluated  by  numerical  calculation  of  the 


11 


Stratton,  op.  cit.,  p.  467- 
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integral.  The  fields  of  a  "double  layer"  are  particularly 
troublesome,  because  the  calculation  of  the  normal  deriva¬ 
tives  involves  use  of  the  derivative  of  the  equation  for  the 
curve  C.  If  there  exists  some  reason  to  believe  the  fields 
are  relatively  insensitive  to  changes  in  the  surface,  a 
simplified  representation  would  be  possible;  however,  if 
the  field  depended  greatly  on  the  shape  of  the  surface,  no 
appreciable  simplification  would  be  possible.  In  the  problem 
considered,  the  fact  that  the  electromagnetic  wave  is  attenu¬ 
ated  rapidly  in  sea  leads*  to  the  belief  that  a  simplified 
representation  should  be  possible. 

Consider  representing  cp(p')  using  sources  of  a  single 
layer  only. 


<p(p')  .= 


1  £ (p  )  G(kcj?-p'|)<U  • 
C 


(2.5-4) 


If,  in  the  sea,  f  can  be  considered  reasonably  constant  the 


is  given  approximately  by 


<p(p) 


2k 


f(p) 


(2.5.5) 


f(p)  - ~r  <p(p)  -  (2.5.6) 

Equation  (2.5-6)  yields  the  function  f  needed  such  that  the 
representation  of  equation  (2. 5-^0  tends  to  correct  boundary 
value  on  C.  However,  this  does  not  mean  (2.5-^-)  is  a  correct 
representation  of  the  solution  to  the  boundary  value  problem, 
cp  (p  1  )  has  yet  to  be  considered. 
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<pn(5')  = 


HI)  Gn(kc!o-?'l)dt 


1  1  2k  k 

«-  Ho')  =  -  (-T5)  <P  (o  ‘  )  =  — f  <p(p*  ) 

2  2  1  1 


(2.5-7) 


Equation  (2.5*7)  represents  the  impedance  boundary  condition  of 

Section  2.4, 

then 


P  2k 

9 ( p ' )  =  j_  -r  *(S)  g(*cI?-p,I)<u. 


(2.5*8) 


Equation  (2.5*8)  is  an  accurate  solution  to  the  boundary  value 
problem  in  the  sea  only  under  restricted  conditions?  the  condi¬ 
tions  are  basically  those  assumed  in  Section  2.4. 

2.6  Comparison  with  Other  Theories. 

The  basic  results  from  the  solution  of  the  integral  equa¬ 
tions  ±n  the  air  are: 

TM  Case 

®  (x)  2ik  .  (l  +  k  A  cos  k  x)e’llcxix  =  2cp  . (x)(l+k§(x)). 
Ynv  1  yiv  s  s  '  ymv  7fv  s*v  11 


TE  Case 


<p(x)  m  2e_i^kxiX  +  kyi^  »  2cpi(x)f. 


(2.6.1) 


(2.6.2) 


In  the  TM  Case,  Lerner  and  Max  and  Morgan  use  conformal 
mapping  of  the  static  (co  =  0)  problem  to  obtain  the  tangential 
magnetic  field?  the  justification  of  such  a  procedure  is  con¬ 
sidered  in  Section  2.8.  The  results  of  Lerner  and  Max  basical¬ 
ly  agree  with  the  result  given  above.  The  agreement  is  not  com¬ 
plete,  however,  as  Lerner  and  Max  consider  only  grazing  inci- 
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donee,  and  the  above  solution  includes  all  except  grazing  inci¬ 
dence.  The  tangential  magnetic  field  is  obtained  from  equation 
(2.6.1)  by  simply  setting  Htan  =  cpn«  Morgan's  physical  model 
includes  %'  ±«o,  which  is  excluded  from  the  model  of  the  problem 
used  here.  Morgan  states,  however,  that  there  will  be  a  change 


in  the  fields  in  the  TM  case.  Wait  assumes  that  the  fields  on 


the  surface  are  relatively  unchanged.  This  is  of  course  ap¬ 
proximately  true  even  in  the  TM  case  (at  least  as  far  as  the  or¬ 
der  of  magnitude  of  the  fields  is  concerned)  for  the  usual  set 
of  physical  parameters  considered  here. 

In  the  TE  case.  Wait  and  Morgan  again  assume  on  physical 
grounds  that  the  fields  are  unchanged.  Lerner  and  Max  obtain  a 
similar  result  by  solution  of  the  static  problem. 

Another  general  method  used  to  solve  the  "rough  surface" 

12 

problem  is  the  perturbation  theory  (see  Appendix  B).  Winter 
used  the  simplest  form  of  this  technique  to  obtain  the  fields 
on  the  sea  surface.  The  method  used  by  Winter  is  applied  below 
to  the  scalar  functions  to  learn  if  the  results  agree. 

"Since  the  roughness  scales  of  the  sea  surface 
are  extremely  small  compared  with  the  wavelength, 
the  electromagnetic  fields  in  free  space  are 
scarcely  altered  by  the  surface  irregularities. 


winter,  D.  F.,  "Low  Frequency  Radio  Propagation  into  a 
Moderately  Rough  Sea,"  Radio  Propagation-- Section  D,  Journal 
of  Research,  National  Bureau  of  Standards,  Vol.  67/  no.  5, 


Sept. -Oct.  1965/  P*  551* 
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Hence,  the  derivatives  of  <p(z)  can  be  calculated 
from  the  solution  to  the  smooth  sea  problem  with 
little  error  ..." 


This,  as  stated  in  the  appendix,  is  not  necessarily  true. 
The  above  quotation  is  probably  true  in  the  far  field;  however, 
there  is  some  question  as  to  its  correctness  in  the  near  field. 
TM  Case: 

<p(x,o)  =  0  qp(x, y)  =  2i  sin  k  ^y  e""^xix 

<py(x.s)  =  <Py(x,o)  +  q>yy(x,o)S  +i  (pyyy(x,o)52  +  .  .  . 

=  (2i*yi  -  0  -  2ikyi|  (ky.?)2  +  .  .  .)  e-ikxix 
«  2ikyle-ikxix 

<Px(x,5)  =  <Px(x,o)  +  9xy(x,o)?  + -|  9xyy(x,o)?2  +  .  .  . 

-  (0  +  2i(-ikxi)kyiS  +  0  +  .  .  .)  e-ikxix 

»  2 ikyie’ lkxiX (' ikxi? ) 

<Pn(x<§)  2ikyi(1  +  kxi?  Sx )  e~ lkx iX  «  2iky^e~lkxix 
<&  2ikyi(l  +  ksl)  e‘ikxix. 

The  perturbation  method  yre Ids ,  as  as s\imed,  an  "unperturbed 
field;"  unfortunately,  this  is  not  necessarily  correct  in  this 


case. 
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Perturbation  theory:  H  «  2ik  ?  (-ik  ?) 

y  y  x v  x5 ' 

Results  above:  H  &  2ik  § 

y  y*x 

TE  Case: 

cpy(x,o)  =  0  a>(x,y)  =  2  cos  kyiy  e“ikxix 

<p(x,y)  =  <p(x,o)  +  <p  (x#o)  §  +  cp^Cx^o)?2  +  ... 

=  (2  +  0  +  -i(-k2j  )?2  +  . . .  )  e~^xix 

c.  y.4. 

0  -  ik  .  x 

«  2  e  xi 

which  is  the  unperturbed  field  and  which  agrees  with  the 
result  obtained  in  this  investigation. 

2*7  Estimates  of  the  Effect  of  the  Finite  Conductivity  of 
the  Sea  on  the  Fields  in  the  Air. 

As  was  noted  in  the  section  describing  the  pnysical 
model  (1.3)/  the  fields  on  the  sea  surface  have  been  com¬ 
puted  under  the  assumption  that  the  sea  was  a  perfect 
electrical  conductor  (i.e.,  surface  impedance  was  zero). 

In  this  section,  the  effect  of  finite  conductivity  on  the 
solutions  is  considered.  The  impedance  boundary  condition 
obtained  in  Section  2.4  may  be  used  to  estimate  the  effect 
of  the  finite  conductivity  of  the  sea  water  (and  therefore 
non-zero  surface  impedance)  on  the  solutions. 

The  integral  equation  that  the  wave  functions  must 
satisfy  is 

q>(p)  =  2?i(p)  +  2^  {cp(p')  Gn*  (k0|p-p  1  |  ) 

C 

-  V(p')  8(*oIP-P,l)}0f 


(2.7.1) 
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In  the  TM.  case,  cp  represents  a  rectangular  component  of 
the  vector  potential  A  and  Htan  is  directly  proportional 
to  cpn*  The  boundary  condition  $  =  0  on  C  was  applied  and 
the  solution 

<pn(x)  «  2  ikyi(l  +  kg  A  cos  kgx)  (2-7-2) 

was  obtained. 

Applying  the  impedance  boundary  condition  to  equation 
(2-7-2)  implies  that 
2k  . 

q>/x)  «  (1  +  k  A  cos  k  x)  (2 -7 - > ) 

'  '  c 

then 

f  <p(p  ■  )  Gn.(|p-0'|)-H‘~-|  (\A)2  (2-7-4) 

C  c 

where  Gn, (jp-p' |)  =  lnjp-p'  |,  the  static  two-dimensional 

Green's  function  (the  approximation  of  the  dynamic  Green's 
function  by  the  static  Green's  function  is  considered  in 
Section  2.8). 

As 

^  <pGn,<H*  «  j«pi!  (2-7-5) 

c 

M  «  l<PiI  (2-7-6) 

the  original  approximation  of  <p  =  0  in  C  yields  an  accurate 
integral  equation  and  the  original  solution  (2-7-2)  is 
verified.  In  this  case  there  is  little  change  due  to  the 
introduction  of  the  non- zero  surface  impedance. 
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In  the  TE  case,  cp  represents  a  rectangular  component 
of  the  vector  potential  F.  cpn  was  assumed  to  be  zero  on 
C  and  the  solution 


<p(x)  «2cpi(x) 


(2.7.7) 


was  obtained.  The  impedance  boundary  condition  gives 
*c 

<Pn  -  T'  <P 

Again  using  the  static  kernel. 


(2.7.8) 


<Pn  Gd  i  «  (-^)  2{2  |  ik  } 


<p  Gd  t 
Yn 


«  1 


~ik  k 
y  c 


(2.7.9) 

(2.7.10) 


Again,  the  solution  is  basically  unchanged  from  the  zero 
surface  impedance  case.  There  is,  however,  a  small  correc¬ 
tion  term  in  the  TE  case. 


2.8  The  Justification  of  the  Use  of  the  Static  Kernel  in  Some 

Integral  Equations. 

In  general,  electromagnetic  fields  are  generated  by  and 
support  nonstationary  currents,  that  is,  currents  that  oscil¬ 
late  at  such  frequencies  as  to  make  the  interaction  between 
current  elements  in  different  parts  of  space  significantly 
affected  by  the  finiteness  of  the  velocity  of  propagation  of 
electromagnetic  effects.  The  quasi- stationary  state  is  a 
special  case  of  the  general  nonstationary  state  in  which  the 
velocity  of  propagation  may  be  treated  as  being  infinite. 

In  the  quasi- stationary  state  the  currents  oscillate  slowly 
enough  so  that  the  approximation  that  all  significant 
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interactions  between  currents  are  effectively  instantaneous 

is  accurate.  As  instantaneous  interaction  between  two 

separate  elements  is  the  same  as  the  continuing  or  constant 

interaction  in  the  stationary  state  (in  the  stationary  or 

steady  state,  the  electromagnetic  fields  and  currents  do 

not  vary  with  time;  therefore,  the  interactions  are  unaffected 

by  the  finite  velocity  of  electromagnetic  interactions). 

For  monochromatic  fields  of  radian  frequency  a),  the  condition 

l"5) 

for  the  quasi-stationary  state  is  that:  ^ 


(£>l 


max 


«  V 


max 


« 


2ir 


where 


(2.8.1) 


V  —  velocity  of  propagation  of  electromagnetic  radiation 

o 

in  the  medium  (3*10°  meters/sec.  in  free  space). 

I  —  maximum  distance  between  currents  which  signifi- 
cantly  interact. 

For  the  electromagnetic  fields  in  the  air,  it  is  not  clear 
that  an  satisfying  the  requirements  (2.8.1)  can  be  found. 

In  fact,  it  is  clear  for  <fK  and  q>r  the  quasi-stationary  state 
does  not  hold  as ; 


H 


(k  jx-x'l)  e  i^nks  +  kx^X  dx'  = 


n  -ik  x 
2e  x 


^r(nk  +  k  )2  -k2 
S  X '  o 

(2.8.2) 


x ' 


1"5 

^King,  R.  W.  P.,  "Quasi-Stationary  and  Nonstationary 
Currents  in  Electric  Circuits,"  in  Handbuch  der  Physics, 
Band  XVI,  Elektrische  Felder  und  Wellen  (Berlin,  Springer, 

lSWTpT  165“ - 
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®.  and  <p  are  involved  in  the  n  =  0  term  in  which  k 

may  not  be  neglected,  which  clearly  implies  that  for  these 

fields  the  qua si- stationary  approximation  does  not  hold. 

However,  the  scattered  fields  cp  are  quasi- stationary. 

s 

The  approximation  involved  may  be  seen  in  the  integral 
equation  (2.8.2)  (kg  »  kQ  so  that  J (nko  +  k^)2  “  k2  «  nkg) 
for  terms  |n|  >0. 

The  more  accurate  general  result  that  the  solution  of 

2  2  2 
V  <p^  =  o  yields  an  accurate  solution  to  (v  +  k^)q>s  =  0 

may  be  seen  from  the  boundary  value  problem. 

(V2  +  k2)  <Ps(x,y)  =  0 

<p  (x,y)  =  cp  (x  +  L,y)  boundary  conditions  (2.8.3) 

s  s 

2  <ps  (x  #  y )  =  -  (r^~~ )  <PS  (x » y ) 


o 

■^-7  ?<.(*. y)  = 

W>'  -  <X>2}  *,<*•*> 

(a) 

= 

(¥)2  d-  &>2> 

(b) 

(2.8.4) 

a2 

if  L  «  A 

ay 

_  /  2n7r  \ 2  _ 

<P-  «  (  )  <P_  or 

S  L  s 

V2«p  =  0 

s 

(2.8.5) 

2 

then  solutions  of  V  cp  =0  closely  approximate  the  solutions 

s 

o  o  2 

of  (V  +  k^;9s  =  0,  as  long  as  L  «  A.  As  V  cpg  =  0,  by  the 

use  of  Green's  theorem: 

9s(p)  =  J^(cPg(p')  Gn.(p,p')  -  ^gn'Cp')  G(p'p') 

C 


(2.8.6) 
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where 

G(p  /  p  ' )  =  -to  I  p-p  1  I 

is  the  "static"  Green's  function. 

As  shown  above,  the  stationary  state  should  be  used 

only  to  compute  the  "static"  part  of  the  field  (i.e.,  <p  ), 

s 

and  the  and  cp^  should  be  removed  before  using  static 
approximations . 

The  use  of  periodic  boundary  conditions  played  a  major 

role  in  determining  that  cp  could  be  computed  accurately 

s 

by  use  of  quasi- stationary  equations.  It  therefore  seems 
wise  to  investigate  the  correctness  of  the  assumption  of 
a  periodic  surface.  For  L  «  X,  the  periodic  assumption 
would  not  seem  to  effect  the  fields  in  any  major  way  (i.e., 
a  small  change  in  the  surface  some  distance  from  where  the 
fields  are  computed  will  have  little  effect  on  the  fields), 
in  general,  the  above  statement  is  not  true;  particularly 
when  X  «  L,  the  periodic  assumption  may  lead  to  a  great 
change  in  the  fields. 

2.9  The  Effect  of  Motion  of  the  Sea  Surface. 

It  was  previously  assumed  that  the  sea  surface  was 
stationary.  This  approximation  seemed  to  be  reasonable, 
because  the  velocity  of  the  electromagnetic  wave  in  air 
is  many  times  the  velocity  of  the  sea  wave.  However, 
even  with  the.  great  difference  in  velocities,  there  still 
is  some  effect  on  the  electromagnetic  fields  in  the  air 
due  to  the  motion  of  the  sea  surface.  This  effect  is 
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greatly  multiplied  in  the  sea  as  the  velocity  of  the  elec¬ 
tromagnetic  radiation  is  much  less  in  the  sea  than  in  the 
air.  However,  it  is  sufficiently  accurate  to  solve  the 
boundary  value  problem  for  a  stationary  surface  and  then 
assume  the  computed  field  moves  with  a  velocity  related 
to  the  velocity  of  the  surface. 

In  many  cases,  o  a  stationary  observer  in  the  sea, 
it  may  appear  that  the  sea  surface  is  flat,  but  that  the 
source  of  the  electromagnetic  wave  is  approximately  moving 
up  and  down  in  the  sea.  If  the  wave  were  to  propagate 
into  the  sea  as  e  T  ,  where  D  is  the  depth  below  the  sea 
surface,  this  analysis  would  be  correct?  however,  as  will 
be  shown  later,  this  is  not  exactly  true.  The  difficulty 
with  the  general  case  (i.e.,  using  the  actual  solution  to 
the  boundary  value  problem)  is  the  "equivalent  velocity 
of  propagation"  can  be  obtained  only  numerically,  so  no 
general  result  can  be  stated. 

If  it  is  assumed  that  the  fields  propagate  approxi- 
-  vD 

mately  as  e  ,  then  from  the  usual  doppler  theory 


f 

i  -  v  /c 


where 


(2.9.1) 


f  is  the  frequency  of  the  electromagnetic  wave 
f  the  observed  frequency  (to  a  stationary  observer 
in  the  sea) 

v  the  assumed  phase  velocity  of  the  electromagnetic 
field  as  a  whole  (the  equivalent  velocity  that  the 
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sources  of  the  field  would  have)  and  carries  a 
plus  sign  if  the  source  is  approaching  the  observer 
and  negative  sign  if  the  source  is  moving  away  from 
the  observer 

c  phase  velocity  of  electromagnetic  radiation  in  the 
sea . 


v  may  be  obtained  as  follows: 

The  velocity  that  an  equivalent  source  (that  is,  a  source 
that  would  yield  the  same  fields),  moving  up  and  down  in  the 
sea,  is  given  by: 


v  *  It  D  =  A  cos  (ksvst)  = 


-  Ak  v  sin(kvt)  (2.9*2) 
s  s  s  s 


M  s  ksA  VS 

C  =  =  f(2ir6c)  =  £±§22  =  1600  jr  =  5-104  Jf^rc  (2.9.3) 

yf 

For  f  =  1  kc/s 

c  =  5*10^  meters/sec.,  a  decrease  by  a  factor  of  approxi- 

4 

mately  10  from  the  phase  velocity  in  free  space. 

v  ^  J>0  meters/sec.,  even  for  high  sea  states, 
s 

j v |  ^  JO  =  4. J  meters/sec. 

999*9  *  fQ  £  1000.1  cps.  (2.9.4; 

This  amount  of  doppler  shift  may  not  a^  r  to  be  too  large; 
and  this  is  a  "worst  case"  calculation.  The  "largeness"  of 
this  effect  is  due  to  the  shortening  of  the  wave  length  in 


;he  sea. 


3.0  COMPUTATION  OP  THE  FIELDS  IN  THE  SEA. 


3.1  Introduction. 

In  the  two-dimensional  case  the  problem  of  calculating 

the  fields  in  the  sea  reduces  to  the  solution  of  a  Dirichlet 

type  boundary  value  problem. ^  Unfortunately,  again  because  of 

the  rough  sea  surface,  the  classical  method  of  separation  of 

0 

variables  can  not  be  cirectly  applied.  However,  once  the 
fields  are  below  the  lowest  point  of  the  rough  sea  surface, 
separation  of  variables  can  be  used  directly  to  compute  the 
fields.  The  method  of  separation  of  variables  is  used  below 
to  obtain  some  interesting  and  general  results. 

To  compute  the  fields  in  the  sea,  the  integral  representa¬ 
tion  of  the  wave  functions  may  be  used.  The  integrals  would 
then  be  numerically  evaluated  to  obtain  the  fields.  The  major 
reasons  why  such  an  approach  was  not  taken  are: 

1.  As  the  fields  near  the  surface  were  to  be  computed, 
the  kernel  of  the  integral  would  have  to  be  evaluated  for 
small,  intermediate  and  later  large  arguments.  The  evaluation 
of  the  Hankel  function  of  complex  argument  is  a  somewhat  in¬ 
volved  though  straightforward  problem,  involving  large  amounts 
of  computation. 

- - ^ - 

A  Dirichlet  boundary  value  problem  is  to  find  cp  such 

that 

(V2  +  k2)  cp(p)_  =0 

lim  cp(p)  =  f(p)  a  given  "relatively"  arbitrary  function. 
p->C 
p  €  V 

2 

See  Appendix  B. 
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2.  The  amount  of  computation  is  greatly  increased 
by  the  fact  that  the  integral  must  be  evaluated  for  each  point 
at  which  the  field  is  computed,  or  at  least  for  points  in  the 
near  field(i.e,,  near  the  sea  surface). 

Basically,  for  the  above  reasons,  even  though  the  numeri¬ 
cal  evaluation  of  the  integral  representation  has  been  effec- 
tively  used  previously,^  the  method  was  not  used  here. 

There  still  remains  a  wide  selection  of  methods  that  may 
be  used  to  compute  the  fields  in  the  sea.  The  most  general  of 
these  is  the  method  of  finite  differences  as  used  to  obtain  a 
solution  to  certain  boundary  value  problems.  The  finite  dif¬ 
ference  method  is  a  widely  used  technique  for  the  numerical 

4 

solution  of  boundary  value  problems. 

Briefly,  this  method  consists  of  replacing  the  partial 
differential  equation  with  a  partial  difference  equation.  This 
approximation  involves  an  error,  usually  called  discretization 


^Banaugh,  Robert  P.,  "Scattering  of  Acoustic  and  Elastic 
Waves  by  Surfaces  of  Arbitrary  Shape,"  (Ph.  D.  Thesis,  Univer¬ 
sity  of  Wisconsin),  1962. 

Mei,  Kenneth  Kwan-hsiang,  "Scattering  of  Radio  Waves 
by  Rectangular  Cylinders,"  (Ph.  D.  Thesis,  University  of 
Wisconsin),  1963. 

Lerner  and  Max?  loc.  cit. 

4 

Collatz,  L.,  The  Numerical  Treatment  of  Differential 
Equations,  (Berlin:  Springer,  19b0j . 

Kantorovich,  L.  V.,  and  V.  I.  Krylov,  Approximate  Meth¬ 
ods  of  Higher  Analysis,  (New  York:  Intersdience,  . 

Forsythe,  G.  E.,  and  W.  R.  Wasow,  Finite-Difference 
Methods  for  Partial  Differential  Equations,  (New  ’ 
and  Sons,  196O) . 
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error.  The  problem  then  becomes  one  of  solving  a  set  of  dif¬ 
ference  equations.  Because  of  the  accuracy  needed,  the  number 
of  equations  may  become  quite  large,  on  the  order  of  100  or 
more.  The  technique  usually  used  to  solve  this  set  of  equa¬ 
tions  is  by  interation  procedures,  though  it  should  be  noted 
that  100  x  100  matrices  may  be  directly  inverted  on  computers 

5 

available  today.  In  what  follows  the  general  results  obtain¬ 
able  by  the  classical  method  of  separation  of  variables  are 
first  discussed.  This  approach  should  give  some  insight  use¬ 
ful  in  considering  the  computation  of  the  fields  in  the  sea. 

The  technique  used,  the  so-called  “method  of  lines,"  is  a 
modification  of  the  finite  difference  method.  The  partial 
differential  equation  is  approximated  by  a  difference  differ¬ 
ential  equation.  The  finite  difference  approximation  is  used 
in  the  direction  approximately  parallel  to  the  rough  surface  to 
obtain  the  "propagation"  or  separation  constant  to  be  used  in 
the  differential  equation  which  characterizes  the  fields  in  a 
direction  approximately  normal  to  the  sea  surface. 


z  ®—+~x  x  direction  approximately  parallel  to  sea  surface, 
v  y  direction  approximately  normal  to  sea  surface. 
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Figure  5.1.  Grids  for  Computing  the  Fields  in  the  Sea. 
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3.2  Classical  Separation  of  Variables. 

The  boundary  value  problem  for  the  fields  in  the  sea  may 
be  stated  as  to  find  cp  (a  rectangular  component  of  the  field 
vectors)  such  that  it  satisfies  the  Helmholtz  equation 

(v2  +  k2)  cp  =  0  (3.2.1) 

and 

<¥>(x,y)  =  +  L,y)  elkxiL  «  cp(x  +  L,y)  (3.2.2) 

Equation  (3.2.2)  is  a  periodic  boundary  condition  (i.e,,  it 
requires  cp  to  be  periodic  with  a  period  L  within  a.  constant 
factor  eikxiL) .  This  follows  directly  from  the  physical  prob¬ 
lem,  as  discussed  in  section  2.2.  The  factor  e~i3cxiL  may  be 
neglected  as  kQL  <<  1. 

cp  must  satisfy  the  radiation  condition  in  the  ( +y)  direc¬ 
tion  below  the  lowest  point  (trough)  of  the  sea  surface.  Above 
the  trough  of  the  sea  wave,  both  "inward"  and  "outward"  waves 
may  exist. ^  As  the  sea  is  a  "highly"  conducting  medium,  the 
second  radiation  condition  holds  and  the  fields  decay  exponen¬ 
tially. 

cp  must  take  on  the  correct  boundary  value  on  surface  5(x) 

lira  cp(x,y)  =  f(x) 

(x,y)  ->  C 

(x,y)  e  S  (3.2.3) 

-  7-  . 

°See  Appendix  B. 
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Separation  of  variables  may  be  applied  to  the  boundary 

value  problem  in  the  region  below  the  lowest  point  (trough)  of 

the  sea  surface;  this  is  possible  there  because  the  radiation 

condition  applies j  In  the  region  above  the  trough  of  the  sea 

surface  but  below  the  crest  separation  of  variables  applies; 

however,  there  are  too  many  constants  to  determine  by  classi- 
7 

cal  methods.1  It  is  nevertheless  possible  to  obtain  very  gene¬ 
ral  qualitative  results  by  this  method. 

From  the  boundary  condition 


cp(x,y)  tp(x  +  L,y) 


(3.2.4) 


cp(x,y) 


(5.2.5) 


The  qualitative  results  follow  from  a  discussion  of  the  pro¬ 
perties  of  k  .  Since  k  .  and  k  .  are  negligible  compared  to 

y  x  x  y  x 

k  and  k  : 
c  s 

~  (nks)^  n  >  0  (a) 

ky  ^  kc  n  -  0  ( b)  (3.2.6) 

7 

'See  Appendix  B. 
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For  n  =0,  k  «  kc;  the  field  is  approximately  constant  in  the 
x-direction  resulting  in  a  plane  wave  propagating  approximately 
straight  down  into  the  sea  as  in  the  flat  interface  case. 

If  the  field  is  not  approximately  a  constant  in  the  x- 
direction,  the  propagation  constant  becomes 


Yy  =  */  Y  -  + 


=  J  fz  - 


(3.2.7) 


As  the  term  added  to  increases  in  magnitude.  Re  Fy  1 
increases  and  the  particular  mode  characterized  by  this  pro¬ 
pagation  constant  decreases  more  rapidly  with  depth  than  the 

n  =  0  mode.  This  implies  that  as  y  •>  ®  the  elect romagnetic 

o 

fields  become  a  plane  wave  with  propagation  constant  Y. 
Asymptotically,  (y  -»  00 )  the  major  perturbation  in  the  n  =  1 
mode;  this  is  independent  of  the  shape  of  the  sea  surface. 

The  major  perturbation  near  the  sea  surface  will  depend  on 
the  shape  of  the  sea  surface,  and  will  be  attenuated  with 
depth,  the  fields  tending  towards  a  plane  wave  propagating 
downward  (the  n  =  0  mode) . 

Intuitively,  a  reasonable  breakpoint  in  these  types  of 
propagation  is  when  Re  [y  1  for  n  =  1  is  1.85  ReY  =  . 


As  Re|_Y  I  >  Rei^Yj  the  modes  n  >  0  are  attenuated  more 
than  the  n  =  0Jmode  (plane  wave)  .  Asymptotically  (y-*00)  the 


higher  order  modes  become  zero  more  rapidly  than  the  n  =  0 
mode  (i.e.,  the  ratio  of  the  amplitude  of  the  higher  order 
modes  to  the  n  =  0  mode  approaches  zero  asymptotically) . 
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,rr 


— x 


a /a  +  i 


Re 


->  a  =  1  ->  L  ^  »/7  7rn6 

6 


(3.2.8) 


For  L  £  then  n  =  1  mode  is  considerably  more  attenua¬ 

ted  as  it  propagates  down  into  the  sea  than  the  n  =  0  mode  and 
the  higher  order  modes  are  even  more  rapidly  attenuated.  A 
qualitative  idea  of  what  this  means  is  given  in  the  chart  below: 


f  =  3  CPS  6  =  1^5  meters  U  =  65O  meters 


5  kc/s 

4.6 

20 

18.6  kc/s 

1.85 

8 

50  kc/s 

1.45 

6.5 

For  L  =  L  the  added  attenuation  of  the  lower  order  modes 
should  “rapidly"  bring  the  field  into  a  plane  wave.  In  the  VLF 
range,  for  L  =  10L,  the  effect  of  the  added  attenuation  on  the 
mode  solution  should  be  evident,  even  at  the  depths  considered. 
If  L  ^  100L,  at  the  depths  considered  in  this  report  little  de¬ 
cay  in  the  perturbation  due  to  the  rough  sea  surface  would  be 
expected . 


5.5  Method  of  lines 

In  one  method  of  numerical  solution  of  partial  differen¬ 
tial  equations  (the  method  of  lines),  one  of  the  variables, 
say  x,  is  discretized,  while  the  other  variable  y  is  left 
continuous.  After  the  finite  difference  approximations  are 
substituted  for  the  x  derivatives,  the  partial  differential 


equations  become  a  coupled  system  of  ordinary  differential 
equations,  i.e.,  difference  differential  equations. ^ 

Usually  when  the  method  of  lines  is  used  on  an  automatic 

digital  computer,  the  problem  is  discretized  in  the  y-direc- 

tion  also  as  a  finite  difference  method  is  used  to  solve  the 

10 

system  of  ordinary  differential  equations.  In  the  modifica¬ 
tion  of  the  method  of  lines  used  below,  the  problem  is  discre¬ 
tized  in  the  y-direction;  however,  the  difference  differential 
equations  are  solved  by  a  method  closely  akin  to  separation  of 
variables  and  Euler's  method. ^  In  this  problem,  as  is  usual 
in  the  method  of  lines,  the  number  of  subdivisions  in  the  y- 
direction  greatly  exceeds  the  number  of  such  divisions  in  the 
x-direction,  i.e.,  the  discretization  distance  in  the  x-direc- 
tion  is  larger  than  the  discretisation  distance  in  the  y-direc 
tion.  The  problem  is  then  concerned  with  a  rectangular  net 
with  relatively  long  rectangles .  This  is  very  useful  in  sat¬ 
isfying  the  boundary  condition  at  the  sea  surface  (which  is  a 
slightly  rough  surface,  i.e.,  almost  flat  surface  with  rela¬ 
tively  large  variation  in  x,  resulting  in  only  small  varia¬ 
tions  in  y  =  ?(x)  )  . 

The  boundary  value  problem  was  stated  mathematically  in 
in  equations  (3.2.1),  (3.2.2)  and  (3.2.3). 
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As  shown  in  the  preceding  section,  some  general  qualita¬ 
tive  results  may  be  stated,  based  on  the  relative  values  of 
physical  parameters.  Using  the  numerical  results  obtained 
by  the  use  of  the  method  given  in  this  section,  quantitative 
values  will  be  placed  on  the  terms  "relatively  little 
change,"  etc.,  by  comparison  of  the  theory  to  solutions  of 
some  representative  problems . 

The  method  used  discretizes  the  x  variable  by  replacing 
tiie  x-derivatives  by  a  finite- difference  approximation 

Yx  (x'y)  "  q)(x,y)  D-°-  ’>(x'y>  +  E  (3-3-1) 


where  D.O.  is  a  difference  operator  and  E  an  error  term. 

The  D.O.  depends  on  y,  and  will  be  given  explicitly  along 

with  the  error  term  E  later.  The  partial  differential 

equation  then  becomes  a  set  of  ordinary  differential  equa- 

2 

tions;  the  coupling  of  the  equations  is  through  the  y  term. 

A 

Letting 

cp(x,y)  =  X(x)  Y(y)  as  in  the  separation  of  variables 
method 


+  {y2  -  Yx>}  - 0 


dy‘ 

.2 


2  2 
=  Y  -  Y. 


(3-3.2) 


(3-3.3) 


■y  1  'X 

Intuitively,  because  the  sea  surface  is  only  slightly 
rough  and  because  of  the  relative  values  of  the  physical 
parameters  of  the  problem,  the  propagation  should  be  basically 
in  the  y-direction  (approximately  normal  to  the  sea  surface). 
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Y(y)  =  Cj^e'V  +  C2eV  -  (5.?. 6) 

In  the  computer  program,  it  proved  convenient  to  choose  yy 
such  that  Re  <0,  as  the  propagation  is  in  the  +y  direc¬ 

tion;  the  radiation  condition  then  takes  the  form  =  0. 

Below  the  lowest  point  (trough)  of  the  sea  surface, 
the  radiation  condition  must  be  applied,  in  which  case  the 
method  of  separation  of  variables  could  be  used  in  its 
classical  form.  However,  to  obtain  numerical  results,  the 
use  of  the  classical  method  of  separation  of  variables  does 


not  prove  convenient. 
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The  boundary  value  problem  is  best  solved  numerically, 

even  in  this  case,  where  the  separation  of  variables  could 

12 

be  used.  The  problem  is  a  marching  type  problem.  If 
<p(x,y)  is  known  at  the  lattice  points  for  y  =  y  , 
qp  (x ,  y^  +  6y)  can  be  computed  from  the  equations 

Yx  D-°-  <a> 

Yy  =  HI  -  (b)  (5-3.7) 

'?(x,y1  +  6y)  =  <p(x,y,)  eYy(6y)  (a) 

This  process  is  repeated  until  the  depth  y,  the  great¬ 
est  depth  at  which  the  fields  are  to  be  computed,  is  attained. 

Sea  Surface  Crest 


To  compute  the  fields  between  the  crest  of  the  sur¬ 
face  and  the  trough  of  the  surface,  the  problem  is  consid¬ 
erably  more  complicated.  The  radiation  condition  does 


_  - 

A  marching  type  boundary  value  problem  is  a  prob¬ 
lem  in  which  the  knowledge  of  boundary  conditions  may  be 
used  directly  to  compute  numerically  the  field  near  the 
boundary;  this  may,  in  turn,  be  used  to  extend  the  field 
further  in  a  step-by-step  fashion. 
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not  apply  in  this  region  and  the  method  of  separation  of 
variables,  at  least  its  classical  form,  cannot  be  used. 

The  method  used  is,  of  course,  based  on  the  same 
difference-differential  equations  assuming  some  value  is 
given  to  <p(x,y)  for  all  x,y  in  the  region  considered.  To 
compute  a  new  set  of  values  for  <p(x,y),  let 

Y  =  Vy2  -  y2  (a) 

*  X 

Yx=?th7  D-°-  q)(x'y)  (b) 

(3.3.8) 

<p(x,y  -  Sy)  =  C±  +  C2  (c) 

<P(x,y  +  6y)  =  Cie-Vy26y  +  C2eYy26y  (d) 

The  above  two  equations  are  solved  for  and  C2 ; 
then  a  new  value  of  qp(x,y)  is  computed: 

<p(x,y)  =  C^e^y67  +  C2eYy6y  (3*3*9) 

This  process  is  repeated  for  each  x,y  in  the  region  con¬ 
sidered  until  the  process  has  been  repeated  for  the  whole 
region.  The  process  is  stopped  when  there  is  very  little 
difference  between  the  old  values  of  9  and  the  newly  com¬ 
puted  values  of  9 . 

3.4  Error  Estimates. 

When  using  numerical  methods,  the  question  of  error 
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estimates  is  of  primary  importance.  If  no  reasonable 
error  bounds  can  be  found  the  method  used  and  the  results 
obtained  may  be  useless.  In  this  section,  the  error  es¬ 
timates  for  the  numerical  process  described  in  the  pre¬ 
vious  section  are  given.  Because  a  large  digital  computer 
was  used  carrying  many  extra  significant  figures,  round 
off  error  may  be  neglected.  The  major  error  term  then 
arises  fro.,  the  discretization  process. 

As  different  finite  difference  approximations  were 
used  in  the  different  regions  of  the  problem,  two  sets  of 
estimates  are  g.i>'en. 

Region  1.  (Between  the  crest  and  trough  of  the  sea 
wave . ) 


<p(x  +  *vy)  +  <P(x  -  K  •>.  "  2<*(x/y) 

D.O.  <p(x,y)  = - - - 5 - - - 


where  h  is  the  discretion  constant  in  the  x  variable 

X 


(?•*». 1) 


<Pxx(x,y)  =  d.o.  <p(x,y)  +  B 


where  E  is  an  error  term 


E 


‘I’xxxx(x'Y> 


if  h  is  sufficiently  small  and  cp  smootn  enough  ( it  is 
assumed  this  is  true  in  the  remainder  of  this  section? 
this  is  easily  verified) 


iE 


9 


XX 


1 

T 1 


9. 


xxxx 


9. 


XX 


X 


1 

~  12 


(ivxihxr 


"  ^T" 


) 


the  effective  propagation  constant  for  propa- 


where,  yx  - 
gation  in  the  x-direction. 


6  «  \£- 
x  1  V 
7x 


Let 


--  a  relative  error  term 


Yx  --  the  computed  value  of  yx 

2  ~2  2 

Y  =  Y  +  e 
'  x  Tx  Yx 

e  2  z  ±  (V 

\  6  V 

2  2 
ty  —  the  relative  error  in  y. 

x 


X 


From  the  above,  it  is  clear  that  the  less  the  variation  of 

the  field  in  the  x-direction,  the  larger  h,  may  be  made 

and  still  have  sufficient  accuracy  in  the  calculations. 

However,  the  main  question  of  accuracy  involves  y  and, 

2 

of  course,  cp(x,y)  not  y*. 

2  2 1 

As  |  y I  >  ! Y-J  under  most  sets  of  conditions  for  the 
problem  considered  here. 


-J 


o 

Y^  -  Y 


2  =  yVi 

x  '  v  y  J 


,  v  2 

1  -  -  (-£)  +  ... 
2  Y 


where  y  --  the  effective  propagation  constant  for  propa¬ 
gation  in  the  y-direction 
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1  2 
™  Sx 


linger  most  conditions 


For  h  <  46 
x  x 


2  ^  .01  and  e  ,  £  .0001  . 

Yx  Yy 


The  final  error  estimates  involve  the  field  cp(xfy)  which 
of  course  is  the  estimate  desired. 


p  =  |e  y  y 

€cp  16 


Y  h  ,  -y  (l  -  e  )h  -y  h 

xri  Nr'  v  •  \r  1  \r 


i 


'y  y  -  e 


y  y 


}l  -  lvyUYyhy 


e  ^  -r—  0  h 

cp  6  v  y 

y  y 


For  the  condition  above,  if  h  <56 

y  y 

e<p  ^  “01  " 

Under  most  conditions.  €  is  considerably  less. 

CD 

Region  2.  (Below  the  bottom  of  the  sea  wave.) 


D.Oo  <p(x,y)  = 


jcp(x  +  2hx,y)  +  ®(x  -  2h^fy)} 


12  h 


+  l6{cp(x  +  hx,y)  +  <p(x  -  hx,y)}  -  30<p(x,y) 

y. 

(3*4.2 ) 


12h 


q>  =  D*°*  <P(x,y)  +  e 

aA 

with  E  «  ^  ^xxxxxxhx 


For  hx  <  26x 

€  2  ^  .003 

Yx 

€v  *  .0003  • 

Yy 


For  h  <  506 

y  y 

ecp  *  ,015 • 


s ..  —  relative  error  in  the  field 


In  this  region,  the  first  derivative  is  also  calcu¬ 
lated.  This  allows  taking  the  curl  of  the  field  vectors. 


D.O.  cp  = 


cp(x  +  3hx,y)  -  -  3*vy) 


60V 


-  9{cp(x  +  2tVy)  "  «>(x  “  2hx,y^}  +  ^5{cp(x  +  hx,y)  -  cp(x  -  hx,y)} 


cp  =  D.o.  cp  +  E 


(3.4.3) 


E  ***  -^xr  Cp 


y5  ^xxxxx  x 


F  SI  x 

Dx  75 


l^xxxxxi  ^4  ^  1  ,hxN4 

■— - '  hx"  73  '-H 

^X  '  J  X 


For  h  <26 
x  x 


eD  £  .01. 

x 


Similarly,  for  D 


\  £  73 


For  h  <  26 

y  y 


S  .01. 


Based  on  these  rough  error  estimates,  all  calculations  are 
sufficiently  accurate  (within  1%)  .  Nearly  all  calculations 
are  more  accurate,  as  the  conditions  given  above  are  the 
extremes  of  those  encountered.  It  should  be  noted,  how¬ 
ever,  based  on  the  calculations  made,  the  error  in  each 
computation  usually  took  on  the  maximum  value  computed  for 
it  (i.e.,  the  error  was  very  nearly  equal  to  its  bound 
given  above) . 

3.5  Some  Conclusions  and  Verifications  Based  on  the  Com¬ 
puter  Calculation 

In  the  next  chapter  the  numerical  results  of  the  com¬ 
puter  calculations  are  given;  the  errors  and  the  verifi¬ 
cation  of  approximations  made  in  the  computer  program  are 
briefly  discussed  in  light  of  the  experience  of  the  compu¬ 
ter  runs . 

The  first  question  considered  in  numerical  solution 
of  differential  equations  is:  how  fine  must  the  grid  be 
made  to  obtain  accurate  results?  This  is  usually  deter¬ 
mined  in  the  finite  difference  calculations  by  subdividing 
the  finite  difference  interval  until  no  changes  occur  in 
the  solution.  The  "subdividing"  method  was  used  in  the 
computer  calculations .  The  accuracy  required  was  three 
significant  places .  This  method  was  applied  to  both  in¬ 
tervals  in  the  x  and  in  the  y  variables .  The  effect  of 
interval  change  in  the  x  variable  is  somewhat  small  as 
long  as  the  finite  difference  formulas  are  reasonably 
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accurate.  However,  the  changes  in  the  interval  in  the  y 
variable  could  greatly  effect  the  numerical  result.  This 
sensitivity  occurred  only  i:i  problems  where  the  lower  order 
modes  were  decaying  rapidly.  The  reason  for  this  pheno¬ 
menon  is  that  if  a  large  interval  is  chosen  in  the  y  vari¬ 
able  and  if  the  decay  rate  As  rapid  and  not  "corrected"  in 
the  large  interval,  it  causes  the  field  to  "overshoot". 

This  effect  may  be  cumulative,  in  which  case  the  compu¬ 
ted  field  rapidly  becomes  zero,  or  if  the  interval  is  small 
enough  the  computed  field  appears  to  "hunt"  after  reaching 
equilibrium  (i.e.,  plane  wave).  In  either  case,  the  errors 
are  easily  seen  in  the  computed  values. 

An  assumption  that  is  verified  in  the  computed  re¬ 
sults  is  that  |y.  |  <<  |v  |,  even  where  relatively  rapid 

x  y 

decay  of  lower  order  modes  occur  (  j y  I  is  largest  there, 

|yxl  <  10  j  y  j ,  however  it  should  be  noted  that  this  is 

for  k  A  <  1/7,  if  k  A  was  not  restricted  in  any  way,  the 
s  s 

above  conclusion  would  not  hold) . 

It  would  seem  reasonable,  that  the  computer  program 
would  accurately  compute  the  field,  as  the  program  is 
based  directly  on  the  physical  processes  involved  in  the 
propagation  of  the  wave.  However,  it  is  also  clear,  that 
if  the  original  assumptions  and  restrictions  on  the  para¬ 
meters  considered  do  not  hold,  the  program  may  yield  in¬ 


accurate  results . 
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4.0  SELECTED  RESULTS  OF  COMPUTER  CALCULATIONS. 

4.1  Introduction . 

The  purpose  of  this  chapter  is  to  present  numerical  cal¬ 
culations  of  the  electromagnetic  fields  in  the  sea.  From 
these  calculations,  a  correlation  between  what  the  theory  pre¬ 
dicts  for  the  basic  mode  of  propagation  in  the  sea  (see  Section 
3.2)  and  the  computed  electromagnetic  fields  in  the  sea  may  be 
found.  Selected  results  are  presented,  and  these  should  be 
sufficient  to  see  the  basic  propagation  pattern  in  the  sea. 

As  stated  previously,  the  sea  surface  was  assumed  to  be 
described  by  the  equation  ?(x)  =  A  cos  (k  x) .  For  the  calcu- 
lation  given  in  this  chapter,  A  was  chosen  such  that  the  maxi¬ 
mum  slope  bf  the  wave  was  1/7. 1  These  calculations  should  then 
bound  the  electromagnetic  fields  in  the  sea  produced  in  an  ac¬ 
tual  physical  situation  (i.e.,  the  difference  in  propagation 
patterns  caused  by  the  changes  in  the  fields  due  to  the  rough 
surface  is  maximum) . 

First,  to  see  the  correlation  between  the  "propagation 
theory"  and  the  actual  propagation  effects,  a  set  of  curves 
presenting  the  numerical  solution  of  the  scalar  boundary  value 
problem  at  the  lowest  level  (trough)  of  the  sea  wave  is  given. 

Secondly,  a  set  of  results  presenting  the  electromagnetic 
fields  in  the  sea  at  different  depths  is  presented.  These 

See  Appendix  A. 
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fields  are  normalized  to  see  the  deviation  from  a  plane  wave. 
Lastly,  sane  comments  on  the  numerical  results  are  given. 

4.2  Scalar  Fields  at  the  Lowest  Level  (Trough)  of  the  Sea  Wave 

The  symbol  a  is  used  to  connote  attenuation  ( in  db) : 

--  Relative  ratio  of  the  field  at  the  point  considered 
to  the  field  on  the  sea  surface  vertically  above  it  in  the  TE 
case . 

&2  --  Relative  ratio  of  the  field  at  the  point  considered 
to  the  field  on  the  sea  surface  vertically  above  it  in  the  TM 
case;  this  is  basically  to  present  the  propagation  effects. 

--  -The  "corrected"  value  of  a^,  that  is,  the  relative 
ratio  of  the  field  at  the  point  considered  to  a  fixed  normali¬ 
zation  value;  this  is  to  account  for  the  difference  in  the 

field  on  the  sea  surface  in  the  TM  case  (then  a,  and  a_  on  the 

l  5 

same  scale) . 

ac  --  The  computed  relative  ratio  of  the  field  at  the 
point  considered  to  the  field  on  the  sea  surface  vertically 
above  it  assuming  the  wave  propagates  straight  downward  with 
the  propagation  constant  of  sea  water  at  the  frequency  of  the 
electromagnetic  wave  and  that  the  field  on  the  sea  surface  is 
a  constant  (the  same  normalization  constant  used  to  compute 
and  a^)  . 

is  used  to  connote  phase  shift  ( in  degrees) : 
j61  --  Relative  phase  of  the  field  at  the  point  considered 
to  the  field  on  the  sea  surface  vertically  above  it  in  the  TE 


case . 
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^2  --  Relative  phase  of  the  field  at  the  point  considered 
to  the  field  on  the  sea  surface  vertically  abcce  it  in  the  Tte 

case . 


j8  --  The  computed  relative  phase  of  the  field  at  the 
point  considered  to  the  field  on  the  sea  surface  vertically 
above  it,  assuming  the  wave  propagate^  -tr sight  dev  award  with 
the  propagation  constant  of  sea  water  at  the  frequency  cf  the 
radio  wave  and  that  the  field  on  the  sea  surface  has  a  constant 
(zero)  phase. 

The  curves  presented  are  fors 
f  =  18.6  kc/s 

A  =  .3  meters  L  ~  22  meters 

A  =  1.5  L  -  64 

A  =  4.3  L  -  180 


f  =  3  cps 

a  =  15  I,  =  6^0 

For  sea  water  at  18.6  kc/s,  6  =  1.35  meter.-;  therefore 
based  on  the  "propagation  theory,"  the  fields  under  the  surface 
with  L  -  12  should  show  a  marked  decrease  in  the  perturbations 
caused  by  the  rough  sea  surface  and  indeed  this  is  the  -a-e. 

For  L  s  64,  there  should  be  little  decrease  in  the  perturbation 
and  the  calculated  change  is  very  sxa.ll  a\  ;.an  be  -sen  on  the 
graphs,  even  for  the  4.3  meter  waves. 

For  sea  water  at  3  cps,  6  —  145  meters  and  serve  decrease 
in  the  perturbation  for  the  15  meter  wave  noted. 


♦ 


Figure  4.2  s  Scalar  Field  at  the  Trough  Level  of  the  Sea  Wcve 


Figure  4.3.,  Scalar  Field  at  the  Trough  Level  of  the  Sra  wave 
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4,3  The  Electromagnetic  Fields  in  the  Sea. 

Due  to  the  vector  nature  of  the  electromagnetic  fields 

a  set  of  curves  is  necessary  to  describe  the  electromagnetic 

fields  at  a  fixed  depth.  The  description  of  the  fields  is 

presented  in  terms  of  E  or  H  depending  on  the  polarization 

z  z 

and  direction  of  propagation  of  the  incident  field. 

Again,  a  is  used  to  connote  attenuation  ( in  db) . 

ft,  Relative  ratio  of  the  field  (in  this  case  E  or 
i.  z 

R  )  at  the  point  considered  to  a  normalization  value, 
z' 

ac  --  The  computed  relative  ratio  of  the  field  at  the 
point  considered  to  the  field  on  the  sea  surface  vertically 
anove  it  assuming  the  wave  propagates  straight  downward 
with  the  propagation  constant  of  sea  water  at  the  frequency 
of  the  radio  wave  and  that  the  field  on  the  sea  surface  is 
a  constant . 

/3  again  is  used  to  connote  phase  shift  ( in  degrees)  . 

/31  --  Relative  phase  of  the  field  at  the  point  considered 
to  the  field  on  the  sea  surface  vertically  above  it  in  the 
TE  case . 

j3c  --  The  computed  relative  phase  of  the  field  at  the 
point  considered  to  the  field  on  the  sea  surface  vertically 
above  it,  assuming  the  wave  propagates  straight  downward 
with  the  propagation  constant  of  sea  water  at  the  frequency 
of  the  radio  wave  and  that  the  field  on  the  sea  surface  has 
a  constant  ( zero)  phase . 

A  second  set  of  curves  is  presented  to  represent  the 
other  field  components..  These  are  normalized  by  the  facror 
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±1  /  \ 
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c  z  z 


a  =  { 


20  log 


20  log 


hi  h  , 

1  'c  X| 


E 


E 


a 


20  log 


20  log 


E 


x 


ri  H  . 
c  zi 


phase(H  )  +  45°-  phase  E 
X  z 


phase(E  )  -  45°  -  phase  H 

z 

phase(H  )  +  45°  -  phase  E 
Y  ^ 


|  phase(E  )  -  45°  -  phase  H  -  90° 

V  X  z 

where  phase  H  is  the  phase  of  H  . 

a  denotes  the  deviation  of  the  ratio  of  the  horizontal  fields 
from  the  plane  wave  ( flat  interface)  case  ( in  db) .  a  denotes 
the  ratio  of  the  vertical  field  to  the  horizontal  field  (in 


db)  . 

Similarly,  and  ]3  are  measures  of  the  phase  deviation 
from  the  plane  wave  case  for  the  horizontal  and  vertical 
fields,  respectively. 

On  each  set  of  curves  in  this  section  a  complete  des¬ 
cription  of  the  sea  surface  (A  and  l)  is  given,  along  with 
the  frequency,  polarization  and  direction  of  propagation  of 
the  incident  electromagnetic  wave.  The  depth  D  (below  the 
trough)  at  which  the  fields  are  computed  is  also  given  on 


each  set  of  curves. 
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The  curves  presented  in  this  section  are  all  for  the 
direction  of  propagation  of  the  incident  plane  wave  in  the 
x-y  plane.  The  sets  of  curves  differ  with  respect  to 
polarization.  To  obtain  the  fields  due  to  an  incident 
plane  wave  with  propagation  vector  in  the  y-z  plane  the 
"dual"  of  the  solutions  is  taken  (i.e.,  for  a  vertically- 
polarized  wave,  the  curve  for  the  horizontally-polarized 
wave  must  be  used).  Rather  than  belabor  the  discussion, 
a  very  simple  procedure  will  be  given  in  Section  4.4  for 
using  the  curves  to  compute  the  horizontal  fields  in  the 
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Figure  4.9„  Electromagnetic  Fields  in  the  Sea. 


Vertical  Polarization 


Horizontal  distance  fzcro  crest  of  sea  wave  (in  sea  wave  wavelengths) 
Figure  4.8.  Electromagnetic  Fields  in  the  Sea. 
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Figure  f.12,  Electromagnetic  Fields  in  the  Sea 
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Figure  4.14.  Electromagnetic  Fields  in  the  Sea. 
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Figure  4.15.  Electromagnetic  Fields  in  the  Sea. 


Figure  4.16.  Electromagnetic  Fields  in  t'iie  Sea. 
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Figure  4.20.  Electromagnetic  Fields  in  the  Sea. 
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4.4  Interpretation  of  Numerical  Calculations. 

From  the  curves  of  Section  4.2,  there  is  a  clear  indi¬ 
cation  even  at  these  relatively  small  depths  that  the  basic 
theory  of  propagation  (Section  3*2)  is  qualitatively  correct. 
The  field  for  L  -=  12  meters  at  f  =  18.6  kc/s  (L  =  12, 

./2tt6  =  8)  and  L  =  640  meters  at  f  =  3  cps .  (L  =  640, 

//7tt6  =  o5C)  shows  the  decay  of  the  higher  order  modes . 

The  other  curves  (L  =  64,  180,  640  at  f  =  18.6  kc/s, 

_/7tt6  =  8)  show  the  perturbation  caused  by  the  roughness 
of  the  sea  surface  is  "relatively  unchanged."  The  small 
changes  (decaying  of  the  perturbation)  that  do  occur, 
occur  slowly  and  are  only  beginning  to  become  evident. 

The  curves  of  Section  4.3,  which  describe  the  electro¬ 
magnetic  fields  present  the  picture  of  a  plane  wave  propagat¬ 
ing  downward.  As  Y Y,  the  relationships  for  the  horizontal 

u 

fields  are  basically  that  of  a  plane  wave  (i.e.,  =  v^H^) ; 

the  deviation  from  this  condition  is  the  order  of  tenth  of  a 
db  and  one  degree  at  the  depth  considered.  As  !Yvl  «  IYI, 
the  vertical  fields  should  be  relatively  small;  they  usually 
are  at  least  20  db  below  the  horizontal  fields. 

The  practical  interest  is  then  in  computing  the  hori¬ 
zontal  fields,  the  vertical  fields  being  so  small.  As  the 
horizontal  fields  are  basically  related  by  T)  (neglecting 
a  small  factor)  the  electric  field  may  be  computed  from  the 
magnetic  field  and  vice  versa.  A  very  simple  procedure  for 
computing  the  fields  in  the  sea  from  the  curves  of  Section 
4.3  is; 


1.  Calculate  the  tangential  field  assuming  a  flat  sea 
surface  (the  tangential  electric  field,  for  example). 

2.  E  at  depth  D  is  related  to  the  E  computed  on  the 
flat  surface  by  the  and  for  the  horizontal  polarized 
wave  given  in  Section  4.3  (this  is  independent  of  the  polar 
zation  of  the  actual  incident  wave) . 

3.  Similarly,  E  at  depth  D  is  related  to  the  E  on 
the  flat  surface  by  the  a  and  for  the  vertically 
polarized  wave. 
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5.0  THE  BASIC  THEORY  FOR  THE  THREE-DIMENSIONAL  PROBLEM 


5.1  Introduction  to  the  Three-Dimensional  Problem. 

It  would  seem  clear  that  the  three-dimensional  elec¬ 
tromagnetic  boundary  value  problem  will  be  considerably 
more  complex  and  complicated  than  the  two-dimensional 
electromagnetic  boundary  value  problem  treated  previously. 
Comparing  the  two  problems,  some  striking  differences  and 
similarities  appear. 

In  the  two-dimensional  problem,  the  TE  and  TM  modes 
can  be  completely  decoupled.  With  this  convenience,  the 
two-dimensional  electromagnetic  boundary  value  problem  can 
be  reduced  to  a  set  of  considerably  simpler  uncoupled 
scalar  boundary  value  problems.  In  the  three-dimensional 
electromagnetic  problem,  the  TE  and  TM  modes  can  no  longer 
be  decoupled.  Coupling  occurs  through  the  application  of 
boundary  conditions  on  the  rough  surface.  The  three-dimen 
sional  problem  must  then  be  solved  "all  at  once".  However 
the  coupling  between  the  components  of  the  field  is  "weak" 
The  dominant  magnetic  field  component  is  virtually  inde¬ 
pendent  of  the  other  field  components.  Therefore,  the 
dominant  component  of  the  magnetic  field  may  be  computed 
by  assuming  the  other  magnetic  field  components  are  zero. 
The  field  pattern  is  basically  one  nf  a  "dominantly"  TM 
mode  or  at  least  it  can  be  considered  as  such. 
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A  difference  in  the  mathematical  formulation  of  the  two 
problems  is  that  the  use  of  vector  potentials,  which  proved 
quite  useful  in  simplifying  the  two-dimensional  electromagnetic 
problem,  no  longer  provides  a  clear  means  of  simplification  in 
the  three-dimensional  problem.  The  major  reason  for  this  is  the 
effect  of  the  boundary  conditions  at  a  rough  surface „ 

Other  difficulties  in  the  three-dimensional  electromagnetic 
problem  are: 

i.  The  solution  is  now  represented  as  an  integral  over  a 
two-dimensional  space,  so  that  two  integrations  are  necessary, 
not  just  one  as  previously  needed  for  the  two-dimensional  case. 

ii.  The  solution  is  a  vector  and  a  function  of  two  space 
variables  making  the  solution  more  complicated  and  therefore 
somewhat  more  difficult  to  interpret. 

As  a  result  of  the  above  difficulties,  relatively  little 
work  has  been  done  on  the  three-dimensional  problem.1  The  work 
done  usually  involves  the  use  of  a  perturbation  technique  which, 
as  previously  pointed  out,  may  not  be  "accurate"  in  the  near 
field.2 

The  basic  approach  is  based  on  integral  equation  techniques. 
Many  of  the  assumptions  used  in  the  three-dimensional  problem 
were  discussed  in  the  two-dimensional  problem  and  therefore  a 

Lysanov,  Y.  P.,  "Theory  of  the  Scattering  of  Waves  at 
Periodically  Uneven  Surfaces,"  Soviet  Physics  Acoustics,  Vol . 
no.  1  (Jan. -March,  1958),  pp.  1-8 . 

2Hiatt,  loc.  cit . 


Winter,  loc.  cit. 
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somewhat  abbreviated  discussion  is  given  below. 

The  mathematical  description  of  the  physical  model  is  ba¬ 
sically  an  extension  of  the  two-dimensional  model  to  three  di¬ 
mensions.  It  is  assumed  that: 

1.  The  "basically  spherical"  earth  is  replaced  by  a  "basi¬ 
cally  flat"  earth. 

2.  The  sea  surface  ?(x,y)  is  given  by: 

S(x,y)  =  A  cos  ksxx  cos  hsyy 

3-  The  source  of  electromagnetic  energy  is  a  plane  wave 
incident  on  the  rough  surface . 

4.  For  the  computation  of  the  electromagnetic  fields  in 
the  air  but  on  the  sea  surface,  the  sea  is  assumed  to  be  a  per¬ 
fect  electric  conductor  (a  =  «,  z  =  0) . 

s 

5.2  Formulation  of  the  Electromagnetic  Boundary  Value  Problem 
in  Terms  of  Integral  Equations 

A  convenient  starting  point  in  formulating  the  integral 
equations  for  the  three-dimensional  problem  is  the  vector  form 
of  the  Helmholtz's  formula, 

®(r)  =  {t(r‘.)  Gn,(kJ  r-r'|  )  -  fn,(r')  G(kJ  r-r'l  )}ds' 

(5.2.1) 

S(r)  -  J  {t(r>)  Gn,(kJ  r-r'|  )  -  Sn,(J')  G(kJ  r-r«|  )}ds- 

(5.2.2) 

where 

S  is  closed  surface  bounding  volume  V 


r '  e  S 
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n  --  outward  normal  as  previously  defined. 

'  E  and  H  are  given  in  terms  of  their  rectangular  compo¬ 
nents  . 

Again  dividing  the  surface  S  into  two  parts,  s’  and  Sw,  "S’ 

i 

represents  the  rough  sea  surface  and  Sa  "infinite"  hemisphere. 

Following  the  method  used  in  Section  1.4,  equations 
(5.2,1)  and  (5.2.2)  become 

f(r)  =  E^(r)  +  J_{f(r')  Gnl(kj  ?-r'l  )-fn,(r’)  G(kJ  r-r’l  )}ds' 

(5.2.3) 

3(r)  -  31(r)  +  Jjt(r’)  Gn,  (kj  r-r ’|  ) -8n, (r>)  G(kJ  r-r'|  )}ds- 

S  (5.2.4) 

where  r  e  V,  r‘  e  s’. 

Equations  (5.2.5)  and  (5.2.4)  may  be  shown  to  be  equiva¬ 
lent  to  the  Stratton-Chu  type  equations . 


if  =  If.  -  [*  [”z(n'  x  §)G+(n‘  x  if)  x  grad  G+(n'*^)  grad  G^ds 

"S  '(5.2.5) 

§  =  +  r  [y(n'  x  II)  G-(  n 1  x  it)  x  grad  G~(n5-it)  grad  G]ds 

s  (5.2.6) 

The  P  are  "singular"  integrals  (as  were  the  P  in  the  two- 
J  J  U 

dimensional  problem),  the  details  are  available  in  Appendix  C. 

Applying  the  results  of  Appendix  C  to  (5.2.5)  and  (5.2.6) 

f  =  2[f±  -  J_{z(n'  x  t)G+(n'  x  f)  x  grad  G-t-(n’  •ll)gradG]'dsJ 

S  (5.2.7) 

t  =  2[8±  +  J_{Y(n'  x  ®)G-(n‘  x  S)  x  grad  G+(n'  -i)gradG}ds] 

S  (5-2.8) 
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where  r  and  r'  e  s’. 

The  i?  and  if  gi^en  by  (5.2.7)  and  (5.2.8)  satisfy  the  vec¬ 
tor  partial  differential  equations,  radiation  and  source  con¬ 
ditions  of  Maxwell's  equations?  there  remains  only  the  boun¬ 
dary  condition .  12 

The  boundary  condition  at  a  perfect  electric  conductor 

are: 

i)  the  tangential  electric  field  is  zero:  n  x  if4  =  0  (5.2.9) 

ii)  the  normal  magnetic  field  is  zero:  n*lf+  =  0  (5.2.10) 

where  the  +  sign  on  if*  and  if+  den  ote  a  limiting  process  ( see 
Section  1.2). 

Applying  (5.2.9)  and  (5.2.10)  to  (5.2.7)  and  (5-2.8), 
the  following  integral  equations  are  obtained: 

if  =  2[lL  -  J  z(n'  x  if)  G  +  (n  •  if)  grad  G JdsJ  ^  2  n) 


i  =  2[it1  +  j  _  -  (n‘  x  2)  x  grad  G  as 


(5.2.12) 


6.0  THE  SOLUTION  OF  THE  INTEGRAL  EQUATIONS 


6.1  Introduction 

In  this  chapter,  integral  equation  (5.2.12)  will  be 
solved;  this  solution  gives  the  tangential  magnetic  field 
on  the  air-side  of  the  sea  surface.  The  technic- e  used 
is  similar  to  that  used  in  the  two-dimensional  problem. 

The  integral  equation  in  the  sea  again  yields  an  impedance 
type  boundary  condition  on  the  sea  surface. 


6.2  The  Integral  Equations  in  the  Air 

To  compute  the  tangential  magnetic  field  on  the  sea 
surface,  equation  (5.2.12)  is  considered  (note  that  the 
electric  field  is  absent  from  this  equation,  yielding  an 
integral  equation  in  only  the  magnetic  field). 


For  the  scattered  field  kQ  equal  to  zero  (i.e.,  co  =  0  and 
the  Green's  Function  becomes  the  static  or  stationary 
Green's  Function)  is  a  valid  approximation  (see  Section 
2.8).  In  this  case,  the  integrals  that  must  be  evaluated 
(in  5*2.12)  take  the  form: 
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I  = 


00  n  CO 


.  00  U—  oc 


£(x'y) 


[x2  +  y2  +  ?2(x,y)]  /2 


dxdy 


(6.2.1) 


Following  the  method  used  for  the  two-dimensional  problem 
(6.2.1)  is  approximated  by 


I 


p  03  p  03 
CO  V*.  00 


£(x,y) 


y2)V2 


dxdy 


(6.2.2) 


I  * 

This  approximation  is  vali.d  as: 


i  r  r “f(*.y){- 


(x2  +  y2  +  *2Y>/'2 


(x2  +  y2)3/2 


|dxdy | 


^  C 


f!  00 


f  (x,y) 


u-“  (x2  +  y2  +  ?2)5/2 


dxdy  | 


c  *  1  4  3/2  sup(?2,?2) 
sup  (?2.  ?2;  s  1/50 
-97  s  c  s  1.03  • 

To  simplify  the  computations,  the  problem  i.s  subdivided 
into  two  parts  depending  on  the  polarization  of  the  incident 
wave. 

1.  Vertical  polarization.  Without  loss  of  generality 
it  may  be  assumed  that  the  incident  magnetic  field  is  in  the 
x  direction  and  that  the  propagation  vector  is  in  the 
y-z  plane.  The  incident  field  is  a  TM  to  y  electromagnetic 
wave.  As  the  currents  flowing  on  the  sea  surface  produce 
c  magnetic  field  mainly  in  the  x-z  plane,  the  y  component 
of  the  magnetic  field  is  small  and  the  total  wave  remains 
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basically  TM  to  y.  In  this  case,  the  integral  equation 
(5*2.12)  "simplifies"  to 


H 

H. 

X 

l 

H 

>  =  2  < 

o  > 

y 

H 

Zi 

1°, 

-G  (n  H  -n  H  )  +  G  (-n  H  +n  H  ) 

zx  z  x  x  z'  yv  y  x  x  y' 

2  I  {  G  (-n  H  +n  H  )  -  G  (-n  H  +n  H  )  "Ids’ 
JL«  J_«  1  2  z  y  y  z'  xv  y  x  x  y'  J 


-G  (-n  H  +n  H  )  +  G  (n  H  -n  H  ) 

V  '  y  y  y  2  *  \  2  x  5C  Z  ' 


(6.2.5) 


where 


nx,  n^,  nz  are  the  rectangular  components  of  the  outward 
normal. 


n 


t  _ 


«*<*  ) 


1 - 5 - o— “  «,(-=■) 


similar  equations  hold  for  n  and  n 

y  z 


G^,  Gyr  Gz  are  rectan9u^-ar  components  of  grad  G. 


G  = 
x 


(x-X '  ) 


{(x-x'  )2+  (y-y '  )2+(S(x,y)  -I  (x '  ,y ' )  )2}5/2 


similar  equations  hold  for  G  and  G  . 

y  z 

Following  the  method  used  previously,  let 


m,n 


x,m,n 


<  h 


y/m,n 


z,m,n 


-i(k  ,mx  +  k  ny)  -  ik  .y 


sy 


yi- 


(6.2.4) 


123 


fi  -X 


xi,  m,  n 
h  . 

yi,m, n 


21,  m,  n 


-i(k  mx  +  k  ny)  -ik  .y  //■  0  ,\ 
e  '  sx  xy  J '  e  yiJ  (6.2.5/ 


_  /  \  *  -i  ,  A  /  i  (k  x  +  k  y )  , 

?  (x,  y)  =  A  cos  ksXx  cos  kg,ry  =  -^  |e  v  sx  .  syJ  '  + 

i(k  x-k  y)  ,  i(-k  x  +  k  _y)  ,  i(-k  x  -  k  y)\  ,r  0  r\ 
e  '  sx  sy^ ;  +  e  v  sx  sy-2  '  +  e  v  sx  sy 11 J-  (6.2.6) 

Placing  equations  (6-2.4),  (6.2.5)  and  (6.2.6)  into 
(6.2.3)  with  the  use  of  the  formula 


P"  r"  {xrX,)  e-i('-sxx'+kxyi',)cbc.dy' 

J_co  <J_oo  ( (x-x  '  )  +  (y-y 1  )  )>/2 


2i7r  k 


sx  -i(k  x+k  v) 
- e  '  sx  sy2 1 


(6.2.7)a 


k"  +kA 
1  sx  xy 


■^To  obtain  equation  (6.2.7)  consider 


CO  n  CO 


0  p 

CD  cJ_  03 


[X-X  ' 


((x-x')Z+(y-y')Z)3/2 


e-i(ksxx'+ksyy'>dx'dy' 


-i  ^ksxx+ksy^^  f“  f  *  x  '  gj  (-VX  — %£-  dx  •  dy  ■ 
J-.  J-«  ((x')Z+(v')2)3/2 


*  .  ik  x' 

x  e  sx 


p«  x '  s  iii  k  x  * 
■  =  2i  9  sx 

^o  ( (x  '  )  +  (v  * 


((x')2+(y')2)3/2 


((x')2+(y')2)3/2 


-  2i[  | 


V  [I  V  v 

r(3/2)  sx  % 


Ocsyy')]  Re[y']  >  0 


(Bateman  Manuscript  Project,  Vol.  I,  McGraw-Hill, 
New  York,  1954;  p.  69/  equation  No.  11.) 
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for  =  < 


r ik2i?(x,y) 

o 

o 


v  -  ik  .  y 
>  e  yi^ 


the  followina  results  are  obtained. 


Hx  -  2{l  + 


k  2A 
sx 


hi  +  ksy 


cos  k  x  cos  k 


sx  sy 


,y}e  lkyiY  (a) 


H 

y 


,  k  k  A  ... 

2{~~^=  *sx*  sin  ksyy}  *  Yi' 

hi  +ksy 


H  «  (-2k  A  sin  k  x  cos  k  y)  e~^kyiY 
z  v  sx  sx  xy*  ’  1 


(6.2.8) 

(*>) 

(c) 


p  co  x  1  sin  kgxx 1 

Noting  the  - * -  -  -a— v —  dx  1  is  an  even  func- 

Jo  ((*’)2  +  (y')2}3/2 


tion  of  y'  for  real  y' 


n m  |  ik  x 1 
x 1  e  sx 


((*')*  +  (Y')2)3/2 


dx'  =  2iksx  Ko(ksxly'i) 


y'  -  real  y'  ^  0 


where  K  is  a  modified  Bessel  function, 
o 


~J  w  f  ^  (0 

Ko^ksxly'  I  )e*Y{  dY'  =  2l  K0(’Ksxy' )  cos  ksyY'  ay’ 


=  2[ 


t r 
1 


^rjrroT2) 


] 


(Bateman  Manuscript  Project,  Vol.  I,  McGraw-Hill, 
New  York,  1954;  p.  49/  equation  No.  40. ) 


f"  f“ - - __  e-ikSxx'e-iksyy'  dx'dy' 

■J.»  J.»  ((x-x')2  +  (y-y')2)3/2 

2i7rk 


sx 


J(,\x)2  +  (*=j2 


-i(k  x  +  k  y) 
e  v  sx  syJ ' 


sy ' 
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2 .  Horizontal  polarization . 

Again  without  loss  of  generality,  it  may  be  assumed 
that  the  incident  magnetic  field  is  in  the  x-z  plane  and 
the  propagation  is  in  the  x-z  plane.  In  this  case  the  inci¬ 
dent  field  is  TM  to  y  and  the  total  field  is  basically 
TM.  to  y.  The  integral  equation  becomes 


f  \ 

H 

✓ 

H  . 

X 

XI 

<  H 

>  =  2  < 

0 

Y 

H 

H_  • 

^  Z  / 

k  Z  V 

'  Ji±  e-ikzi?(x'y)  ) 

o 


H.  =  < 


0 

e-ikzi?(x'y> 

o 


-  ik  .x 
e  xi 


(6.2.10) 


Placing  (6.2.5)/  (6. 2,6),  (6.2.7)  ana  (6.2,11)  into  (6.2.10- 
the  following  results  are  obtained: 
k_  .•  sx 

■I  1  +  - — -  cos  k  x  cos  h 

syJ 


zi 


H  «  1  + 


°  */k  ^  +  k  ^ 

sx  sy 


cos  k  x  cos  k  y\e_1^yi^  (a) 

sx  sy-^J  J  v  ' 


k  .  rk  k  A  •  . 

H  «  2-^— — f -  ern  k  x  sm  k y\e  yi/ 

y  ol^T77i  sx  sy  J 

sx  sy 

H  w  -2— ^i{k  A  cos  k  y  sin  k  x )}e*" ^yi^  (c) 

z  k  l  sx  xyJ  sx  J  '  ' 


(*) 

(6.2.11) 
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The  tangential  magnetic  field  is  basically  the  same  for 

horizontal  and  vertical  polarization  (neglecting  the 
k  . 

factor  the  basic  difference  is  in  the  direction  of 

k 

o 

the  fields  (relative  to  the  direction  of  propagation). 
Vertical  polarization: 


tangential  magnetic  field  « 
2{l  +  - 


k  2A 
sx 


k  ^  +  k  2 
sx  sy 


cos  k  x  cos  k  y 
sx  sy 


'l  -  ik  .  y 

J  e  y3. 


(6.2.12) 


Horizontal  polarization: 


tangential  magnetic  field  ** 
k 


2-fi  /l 

k  L 
o 


2a 

SX 


2  • 2 
z  +  k  L 
sx  sy 


cos  k  x  cos  k  yTe~ '*'kxix  (6.2,13) 
sx  sy^J  s 


At  this  point  there  is  an  obvious  need  to  check  these 

results  with  the  results  of  the  two-dimensional  problem. 

There  are  two  cases : 

1.  k  =  0,  k  ±0 
sy  sx 

H  perpendicular  to  the  direction  of  surface  variation 


x 


H 


tan 


=  0  (two-dimensional  case) 
n  v  ' 


-  ik  .  y 

T/1  «* 


pa  2(1  +  k  A  cos  »k  x)e"  yi 
'  SX  sx  2 


which  agrees  with  the  results  for  the  limiting  case  (kgy=0) 

in  equation  (6.2.12). 

2 .  k  =  0 ,  k  5*  0 
sx  sy 

parallel  to  the  direction  of  surface  variation 

X 

H.  Ra  2<P. 
tan 


which  again  agrees  with  the  three-dimensional  result. 


127 


A  major  difference  caused  by  polarization  is  in  the 
kzi 

factor  — present  in  the  horizontally  polarized  case, 
o 

This  factor  is  present  in  the  flat  surface  or  unperturbed 
case. 

As  the  basic  modes  of  propagation  in  the  sea  are 
the  same  for  the  two  or  three-dimensional  problem,  no 
numerical  computations  are  given  in  the  three-dimensional 
case. 


6.5  The  Integral  Equations  in  the  Sea. 

The  same  integral  equations  (5.2.3)  and  (5.2.4)  hold 
for  the  electromagnetic  fields  in  the  sea  as  well  as  in 
the  air;  in  the  sea,  E.  and  I?,  are  equal  to  zero  and  k  =  k  . 

JL  1  C 

The  development  of  this  section  parallels  that  of  Section  2.4, 
only  is  somewhat  more  complicated  because  of  the  vector 
nature  of  the  fields  and  the  fact  that  the  integrals  are 
now  over  a  two-dimensional  space.  However,  the  scalar  prob¬ 
lem  in  the  sea  may  be  treated  first. 

The  integral  equation  for  the  scalar  wave  function 
cp(r)  in  terms  of  its  boundary  values  is 


<P(r)  =  2  {cp(r')  Gn,  (kjr-r1  |  )  -  cpn ,  (r  ’  )g  (kc  |  r-r 1  j  )}  ds  ’ 

§  (6.3.I) 

for  r  and  r 1 eS . 

Using  the  same  assumptions  as  in  Section  J>.4 


S 


n  2tt  p®  -ik  p 
?„,(*')  G(kjr-?'|)ds'  «  „n(i)  1 


n 


d$ 


4irp 


PdP 


9„(5>  .  . 

"  71k-  reS 


(6.3.2) 
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J  <P(r‘  )  Gn,  (3cc  |  r-r  '  |  )ds  '  « 
S 


pdp 

(6.3.3)1 

(6.J.4) 


Pn  _  -(x-x')gx  -  (y-y1  )gy  +  g(x,y)  -  g(x',y') 
p  (x-x')  +  (y-y')  +(§(x,y)  -  g(x‘,y')) 

For  x-x'  and  y-y'  sufficiently  small 

5(x,y)  =  5(x',y')  +  ?x(x-x')  +  gy(y-y'  ) 

+  ?xy(X'X')(y-y,)  +  ?yy(y-y')2} 

+  ... 


P 


n 


P 


*  T 


(r 


X 


1  1 
2  R 


where 


R  -  radius  of  curvature  of  the  section  of  the  surface 

X 

(i.e.  curve)  in  the  x  plane. 


radius  of  curvature  of  the  section  of  the  sur¬ 


face  (i.e.  curve)  in  the  x=y  plane. 
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For  |kc|R  «  1/  the  "impedance"  type  boundary  condition  is 
again  encountered: 


,-<o  (r)>  cp  (r) 

*Cr)  «  ^ 

c  c 

reS 

(6.3.5) 

S(r) 

t(*)  »  i  £ 

reS 

(6.3.6) 

c 


Equation  (6 «3*5)  appears  to  be  a  direct  "carry  over"  from 
the  twc- dimensional  problem  where  the  same  equation  was 
derived.  In  the  three-dimensional  problem,  however,  equation 
(6.3.6)  does  not  directly  imply  Etan  =  T)  Htan  as 

curl  8  /  ^  Stan  * 

Then,  to  obtain  the  impedance  type  relationship  between  the 
electromagnetic  quantities  directly,  the  Stratton-Chu  equa¬ 
tions  are  used.  Only  one  of  the  two  vector  Stratton-Chu 
integral  equations  is  necessary. 

§(r)  =  2 J  {Y(n'xf(r'))  G(k,|r-?'|)  - 

s 

(n'xH(r'))  x  grad  G(kc|r-r’|)  -  (n‘  *H(r '  )  )grad  Gj-ds ' 

(6.3.7) 

(nxl?)  x  grad  G  =  (n*grad  g)e?  -  (grad  G*E?)n  .  (6.3-8) 

The  relation  obtained  directly  involves  only  the  tan¬ 
gential  fields,  so  the  second  term  on  the  right-hand  side 
of  equation  (6.3*8)  may  be  neglected.  Also,  using  the 
boundary  condition  at  a  perfect  electric  conductor  n-H  =  0 
so  the  third  term  on  the  right-hand  side  of  equation  (6.3*7) 
is  taken  to  be  zero. 


^This  follows  from  (curl  -  ~)eT  ^  0* 
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£(r)  =  2  r{Y(Sxf(r'))  G(k  |r-r'|)  - 


(n' *grad  G(kc  |r-r  '  |  )£?(r '  )}•  ds '  (6.3*9) 

(n'xE(r'))  G(kc|r-r'|)  ds  '  «  nxE?  (r )  Jg(1cc  jr-r '  |  )  ds  ' 


fa 


“c2^)  Tilt 


(6.3*10) 


f (n‘ *grad  G(k  | r-r '  |  )  H(r 1  )ds  '  M  H(r)  G  , (k  |r-r'|)ds' 
U_  c  n  c  * 


Fa 


S 

_  -Mm 

2ik  R 
c 


(6.3*11) 


Then  from  placing  equations  (6.3. 10)  and  (6.3. 11 )  into 

(6-3*9) 

Htan  I1  -  ^  Etan  (6’5*12> 

c 

for  |kjR  »  1  (see  discussion  Section  2.4)  . 


h  ~  Y  f 
tar.  Ik  Etan  w  r\ 

r*  >r 


'tan 


(6.5.13) 


as 

ko  »  jjyz  nc  =  jm.  . 

Equation  (6. 3*13)  yields  an  estimate  on  the  tangential  elec- 
trie  field  at  a  good  electric  conductor  in  terms  of  the 
tangential  magnetic  field  there. 


7.0  CONCLUSION 


Fran  the  solutions  given  in  Chapters  2  and  6  for  the 
tangential  magnetic  field  on  the  rough  sea  surface,  the 
perturbation  of  this  field  due  to  the  roughness  of  the  sea 
surface  is  less  than  3  db  for  most  sea  conditions.  The 
perturbation  depends  both  on  the  polarization  of  the  incident 
field  and  the  direction  of  propagation  of  the  incident  field. 

From  both  the  numerical  solutions  ( Chapter  4)  and  the 
heuristic  theory  of  propagation  in  the  sea  (Chapter  3),  it  is 
seen  that  the  perturbation  of  the  fields  in  the  sea  caused 
by  the  roughness  of  the  sea  surface  decays  rapidly  with  depth 
if  the  sea  wave  wavelength  is  less  than  or  the  order  of  mag¬ 
nitude  of  the  skin  depth  of  the  sea  at  the  frequency  considered; 
if  the  sea  wave  wavelength  is  many  orders  of  magnitude  larger 
than  the  skin  depth,  there  is  little  decay  (at  the  depths 
considered)  of  the  perturbation,  so  that  the  phase  and 
amplitude  of  the  fields  in  the  sea  vary  directly  with  the 
height  of  the  sea  vertically  above  them.  These  results  are 
in  agreement  with  the  work  of  Lerner  and  Max.'*' 


- y - 

Lerner,  R.  M.,  and  J.  Max,  "Very  Low  Frequency  and  Low 
Frequency  Fields  Propagating  near  and  into  a  Rough  Sea,"  a 
paper  presented  to  the  URSI  Spring  1963  Meeting. 
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1J>2 

The  above  results  are  also  in  agreement  with  the  work  of 

2 

Winter.  Winter  used  a  statistical  description  of  the  electro¬ 
magnetic  fields  in  the  sea,  obtaining  the  basic  conclusion 
that  the  fields  in  the  sea  are  "on  the  average"  greatly  per¬ 
turbed  by  the  rough  sea  surface  (from  the  curves  of  Chapter  -4, 
the  perturbation  may  be  40  db  or  greater) „  The  statistical 
description  of  the  electromagnetic  fields  in  the  sea,  in 
this  case,  do  not  convey  very  much  information.  This  is 
particularly  true  in  light  of  the  fact  that  the  sea  surface 
varies  so  slowly  (the  sea  wave  velocity  is  very  small  compared 
to  the  velocity  of  the  electromagnetic  wave  in  the  sea) , 
permitting  the  observer  to  follow  the  variation  of  the 
sea  wave  by  observing  the  changes  in  the  electromagnetic 
field . 

For  VLF  signals,  the  perturbation  of  the  fields  in  the 
sea  due  to  large  sea  waves  is  relatively  unchanged  (within 
a  few  db)  at  depths  of  ten  meters  or  so.  The  large  perturba¬ 
tion  in  the  signal  due  to  the  rough  sea  surface  could  then  be 
"corrected"  in  part  by  monitoring  the  sea  surface  height  above 
the  receiving  antenna. 

No  "correction"  is  necessary  for  the  perturbation  due  to 
small  sea  waves  as  these  perturbations  are  small  ( a  few  db) . 

For  ELF  signals,  even  the  perturbation  due  to  large  sea  waves 
is  small,  so  no  "correction"  is  needed. 

^Winter,  D.  F . ,  "Low-Frequency  Radio  Propagation  into  a 
Moderately  Rough  Sea,"  Radio  Propanation,  Vol.  670,  No. 

Sept. -Oct.  1965,  PP.  551-562. 


APPENDIX  A 


GENERAL  NATURE  OF  THE  SEA  SURFACE  WITH  A  DISCUSSION 
OF  THE  MATHEMATICAL  DESCRIPTION  OF  THE  SEA  SURFACE 

In  this  appendix,  the  nature  of  the  sea  surface  is  discussed. 
The  present  state  of  knowledge  about  the  sd?a  surface  is  given? 
particular  attention  is  given  the  statistical  nature  of  the 
mathematical  description  of  the  sea  surface.  Also,  at  the  end 
of  the  appendix  is  presented  a  brief  list  of  useful  formulas 
pertaining  to  the  mathematical  description  of  the  sea  surface 
and  some  tables  giving  pertinent  data. 

Even  casual  observation  shows  the  great  irregularity  of 
the  sea  surface?  no  single  wave  retains  its  identity  long?  the 
period,  form,  etc.,  vary  greatly  even  for  consecutive  waves. 
Indeed,  the  sea  surface,  in  a  rough  sea,  seems  to  vary  almost 
randomly  in  both  time  and  space. 

The  study  of  the  sea  surface  can,  for  purposes  of  the  dis¬ 
cussion  below,  be  placed  into  three  general  categories: 

i.  The  study  of  the  sea  surface  by  classical  hydrodyna¬ 
mics  . 

ii.  The  study  of  the  sea  surface  by  probabilistic  methods. 

In  this  category  only  the  linearized  problem  will  be  consi¬ 
dered  as  this  allows  us  to  obtain  general  results.  By  assum¬ 
ing  a  linearized  free-surface  boundary  condition,  the  problem 
becomes  linear  and  the  sea  surface  can  then  be  described  by  a 
known  random  process  (a  stationary  Gaussian  process) . 
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iii.  Other  theories  which  consider  the  non-linear  effects 
of  the  sea  surface,  methods  of  generating  the  sea  surface,  etc. 

In  this  case  the  results  obtained  are  far  less  useful  for  a 
description  of  the  sea  surface  than  the  results  of  Category  ii. 
However,  such  results  place  the  limitations  of  the  theory  of 
Category  ii  in  the  correct  perspective . 

Before  going  into  a  more  detailed  study  of  the  sea  surface 
a  few  general  comments  are  appropriate.  First,  the  present 
theories  assume  a  fixed  meteorological  condition.  Then  the 
description  of  the  sea  surface  given  by  the  theory  holds  as 
long  as  this  assumption  is  approximately  true.  A  theory  based 
on  the  correlation  of  the  changes  in  meteorological  conditions 
with  changes  in  the  sea  surface  could  be  used  but  would  be 
somewhat  complex.  This  may,  however,  be  necessary  in  some 
cases.  Second,  the  presently  known  theories  consider  only  the 
gross  meteorological  conditions  (e.g.,  average  wind  velocity, 
average  fetch,  etc.)  .  In  practice,  this  is  all  that  can  be 
assumed  without  making  the  problem  inordinately  complex.  Third, 
no  theory  at  present  is  complete  in  the  description  of  the  sea 
surfaces.  Under  certain  meteorological  conditions  one  theory 
may  be  approximately  correct?  however,  it  fails  when  the  meteoro¬ 
logical  conditions  change.  It  should  be  noted  that  no  theory 
now  available  gives  the  complete  description  of  the  detailed 
properties  of  the  sea  surface?  only  the  gross  features  can  be 
mathematically  described.  A  brief  discussion  of  the  three 
general  categories  follows: 
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i.  Classical  Hydrodynamics,  When  the  depth  of  the  sea  is 
large,  a  solution  of  the  hydrodynamic  equations  which  represent 
the  sea  surface  is  a  trochcidal  wave.^  The  parametric  equation 
of  the  trochoid  is: 

x  =  r0  -  a  sin  0 

z(x)  =  r  -  a  cos  0  (A.l) 

where 

r,  a  --  fixed  parameters  for  a  given  trochoid. 

0  —  parameter  that  varies  (i.e.,  generates  the  curve). 

The  trochoid  is  a  two-dimensional  wave  which  "could  exist" 
in  a  swell . 

A  more  general  deterministic  model  ( as  opposed  to  a  sta¬ 
tistical  model)  of  the  sea  state  would  include  the  irregularity 
of  the  sea;  however,  such  a  formulation  is  too  complex  for  prac¬ 
tical  ,w3.  The  deterministic  model  is  considered  below.  Leav¬ 
ing  out  details,  the  following  equation  of  a  simple  harmonic 
progressive  wave  is  a  solution  to  the  linearized  hydrodynamic 
equations: 

?(x,y,t)  =  A  cos  (x  cos  0  +  y  sin  0)  -  ~  t  +  ej  (A. 2) 
where 

§  --  sea  surface 

A  —  amplitude  of  the  simple  harmonic  progressive  wave 


xLamb,  H.,  Hydrodynamics,  Dover  Publications,  Inc.,  New 
York,  1945,  P.  4 7JT" 

Kerr,  D.  E.,  Propagation  of  Short  Radio  Waves,  McGraw- 
Hill,  New  York,  1951,  p.  48f5 


J 


T  --  period  of  the  simple  harmonic  progressive  wave 
L  --  wavelength  of  the  simple  harmonic  progressive  wave 
c  --  velocity  of  the  simple  harmonic  progressive  wave 


9  --  direction  of  propagation  measured  with  respect  ot 
the  tx  axis 

e  —  phase  atx=y=t=0  (arbitrary) 


If 


p  =  xx  +  yy  —  position  vector  in  the  horizontal  plane 
it  =  k  x  +  k  y  —  propagation  vector  in  the  horizontal  plane 

S  X  V 


k  =  lit 

s  s 


^  2  2 
—  =  —  wave  number  of  the  sea  surface^ 

g  L 


g  —  acceleration  due  to  gravity 


co 


2ir 

=  -rfr  —  radian  frequency  of  the  wave 


then 


S(P,t)  =  A  cos  (it- p1  -  cot  +  s) 


(A- 5) 


P  =  (x,y)  in  the  ordered  pair  notation  used  previously. 

If  we  assume  the  waves  are  progressing  in  the  +x  direction 

(i.e.,  k  ^  0)  which  is  reasonable  if  the  wind  is  in  the  +x 

direction,  the  general  solution  is-^ 
tt/2 


§(p,t)  =  j1  J*  [a(co,  9)  cos(it-o’-cot)  +  b(co,  9)  sin(it-^p-cot)  J 
-tt/2  0 


dco  d9 
(A. 4) 


■  •  > 

^Longuet-Higgins,  M.  S.,  "The  Directional  Spectrum  of  Ocean 
Waves,  and  Processes  of  Wave  Generation,"  Proceedings  of  the 
Royal  Society,  Vol  265,  no.  1J22,  Jan.  20,  1962,  p.  2bb. 

■3 

^Pierson,  W.  J.,  Jr.,  "Wind  Generated  Gravity  Waves," 
Advances  in  Geophysics,  Vol.  II,  1955,  Academic  Press,  New 
York,  p.  107' 
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where 

a(a),8)  and  b(co,  9)  are  the  spectra  of  S(p,t) 

If  ?(0,y,t)  is  known,  a(<n,t)  and  b(oo,t)  can  be  obtained 
and  from  them  ?(s,y,t).  However,  it  is  clear  that  the  deter¬ 
ministic  model  of  the  sea  surface  as  given  above  is  not  prac¬ 
tical.  As  has  been  the  case  with  many  problems  whose  complexity 
defies  deterministic  solution,  one  next  attempts  to  formulate 
the  problem  in  probabilistic  terms .  We  will  now  consider  such 
an  attempt . 

ii.  Probabilistic  Description  of  the  Sea  Surface.  In  the 
past  few  years  there  has  been  an  increased  tendency  to  treat 
many  natural  phenomena  as  random  processes .  The  main  feature 
of  such  a  process  is  an  indeterminacy  in  the  expected  behavior 
of  a  single  occurrence  coupled  with  strong  statistical  proper¬ 
ties  for  a  large  number  of  occurrences .  The  " indeterminacy"  in 
the  sea  surface  comes  from  the  complexity  of  its  mathematical 
description  Ce.g.,  finding  5(0, y,t)]. 

Chronologically,  experimental  data  first  led  to  the  asser¬ 
tion  that  the  sea  surface  could  be  represented  approximately  as 
a  stationary  multivariate  Gaussian  process.  From  a  theore¬ 
tical  view,  it  can  be  shown  by  using  equation  (A. 2)  with  the 

assumption  that  the  random  vs  -iable  e(o,9)  has  a  uniform  dis- 

4 

tribution,  that  the  general  solution  for  _.ie  sea  surface  is: 

00  ir 

?(p,t)  =  P  J*  cos  (^‘p--  oot  +  e(a>0)v/A^((jo,  9)  dco  d9  (A.  5) 

0  -7T 

Ibid.,  p.  122. 
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This  expression  is  a  multivariate  Gaussian  process  where 


2 

A ‘(to,  9)  is  the  energy  density  spectrum.  This  [as  does  (A. 3)1 
assumes  linearized  equations  and  boundary  conditions;  if  these 
assumptions  hold,  the  surface  is  given  by  (A. 3) .  However,  as 
stated  above,  the  experimental  data  implies  that  the  surface  is 
only  approximately  Gaussian,  the  error  being  due  in  large  meas¬ 
ure  to  the  nonlinear  effects  neglected  by  this  theory. 

The  same  results  can  be  given  in  the  form: 


5(p,t)  = 


-7T 


a(o).  0)  cos  (it-'p  -  ovt  +  e(o>,  0)  djo  d0  (A. 6) 


0 


where 


1  2 

y  A  (to,  9)  =  H(to,0)  is  the  energy  density  spectrum. 


Using  this  Gaussian  mode1  of  the  sea  surface,  many  of  the 

5 

general  properties  of  the  "sea  state"  have  been  calculated. 

The  basic  question  to  be  resolved  is:  Can  the  multivariate 
Gaussian  process  describe  accurately  the  real  sea  surface?  If 
not,  from  what  standpoint  is  it  deficient? 

To  answer  this  question,  we  consider  a  specific  model  of 
the  sea  surface  (i,e.,  a  given  spectrum  E(to,0))  and  then  compare 


9 Ibid . ,  p„  93. 

Longue t -Higgins,  M.  S.,  "The  Statistical  Analysis  of  a 
Random  Moving  Surface,"  Trans.  Royal  Society  of  London,  Series 
A,  Vol.  249,  1936-57,  P.  3^1- 

Longue ' -Higgins ,  M.  S.,  "Statistical  Properties  of  an 
Isotropic  Random  Surface,"  Trans.  Royal  Society  of  London, 
Series  A,  Vol.  250,  1957-58,  p,  157* 
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this  random  process  with  the  known  properties  (experimental 
data)  of  the  sea  surface.  Naturally,  for  a  different  spectrum 
the  properties  of  the  sea  predicted  by  the  random  process  will 
be  different.  The  best  known  spectrum  of  the  "sea  surface"  is 
a  semi-empirical  expression  given  by  Neumann  for  a  fully  devel¬ 
oped  sea:^ 


E(co)  da)  = 


7 r 


7 


0  2/  2  2 
-2g  /oo  v 


dco 


03 


(A. 7) 


where: 

v  is  the  velocity  of  the  wind  "generating"  the  "sea," 
and  c  is  a  constant . 

The  total  energy  E  for  this  spectrum  then  becomes: 


E  = 


CO 

P 


-  ^/2 

E(m)  d03  =  ^73-5  3v 


u  a' 

0  c  g 


(A.8) 


This  seems  to  be  in  agreement  with  some  experimental  data  if 
2  S 

c  =  3*05  m  /sec  .  To  study  in  detail  what  this  model  predicts, 
we  will  consider  the  shape  of  the  spectrum  and  some  of  its  re¬ 
sults  . 


Neumann's  spectrum  rises  rapidly  at  co  =  ^—7—  and  has  a 

-6 

g/v.  ror  xctrge  co,  oh  go;  a 

5  ' 


maximum  at  a  =  J~ 7  g/v.  For  large  co,  co(co)  »  co  and  for  small 


9  0  2 

00  e(co)  ~  e-2'  ~/<>Sv  (i.e.,  the  exponential  predominates),  and 

therefore,  there  is  very  little  energy  in  the  low  frequency 
(long  wavelength)  part  of  the  spectrum. 

Neumann’s  spectrum  is  for  a  fully  developed  sea  only.  In 
actuality,  of  course,  the  sea  may  not  be  fully  developed.  The 


Pierson,  op.  cit 


p.  148 . 
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growth  of  the  sea  waves  depends  upon  the  fetch  (distance  over 
which  the  wind  blows),  the  duration  (length  of  time  the  wind 
has  been  blowing)  and  naturally,  the  wind  velocity.  There  are 
methods  that  consider  the  problem  of  duration-limited  and 
fetch-limited  seas  (but  only  approximately) . 

From  empirical  data  an  approximate  directed  spectrum 
(which  depends  on  the  direction  relative  to  the  wind  as  well 
as  on  the  duration  and  fetch  of  the  wave)  is:*^ 


E(co,0) 


,2,2  2 

-2g  /a)  v 


CD 

0 


cos 


0  for  -  X  <  9<4rcoJ  <  co  <  “ 


t r 

7  "i 


otherwise 


(a.9) 


where 


0  --  polar  angle  with  reference  to  the  x  axis  (the  wind 
is  blowing  in  the  +x  direction) . 
cd^  -  intersection  radian  frequency  (a  function  of  duration 
and  fetch) . 

The  above  model  is  at  best  a  good  approximation  to  the 

sea  surface  under  certain  conditions .  There  are  better 
o 

models;  however,  they  are  more  complex  concerning  the  direc¬ 
tional  part  of  the  formula  and  other  errors  have  been  noted. 


^Ibid . ,  p.  155. 

Q 

°Longuet -Higgins,  M.  S.,  "The  Directional  Spectrum  of 
Ocean  Waves,  and  Processes  of  Wave  Generation,"  Proceedings 
of  the  Royal  Society,  Vol.  265,  no.  1^22,  Jan.  50,  1962,  p.  286. 

Pierson,  W.  J.,  Jr.  (Ed.),  "The  Directional  Spectrum  of  a 
Wind  Generated  Sea  as  Determined  from  Data  Obtained  by  the 
Stereo  Wave  Observation  Project,"  New  York  University  Meteoro¬ 
logical  Papers,  Vol.  2,  no.  6,  June,  1950  (unpublished) . 
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This  is  particularly  true  for  high  sea  states  where  the  sea  may 
not  be  fully  developed.  Other  empirical  formulas  are  available 
to  represent  this  case it  is  still  not  known  conclusively 
how  well  the  statistical  results  of  this  model  hold  for  actual 
sea  waves.  This  model  of  the  sea  surface  contains  some  energy 
in  the  high  frequency  part  of  the  spectrum  and  predicts  the  sea 
to  be  "completely  covered  by  short  wavelength  ripples."  Physic¬ 
ally  this  seems  reasonable  even  though  the  above  spectrum  is 
known  to  be  deficient  in  the  high  frequency  range.  The  error  in 
the  high  frequency  part  of  the  spectrum  is  due  to  the  fact  that 
non-linear  effects  are  prevalent  at  these  frequencies?  we  con¬ 
sider  this  effect  in  the  next  section. 

iii)  The  Generation  of  the  Sea  Surface  and  Non-Linear 
Effects . 

Phillips^  has  described  the  generation  of  the  sea  surface 
by  considering  turbulent  fluctuation  of  the  air  (wind)  above 
the  air-sea  interface.  At  present  this  theory  seems  to  be  ac¬ 
cepted  as  correct  for  the  " original"  formation  of  the  sea  sur¬ 
face.  However,  once  the  waves  have  been  formed,  their  growth 
occurs  through  other  mechanisms . 

_ 

Pierson,  W.  J.,  Jr.,  "A  Study  of  Wave  Forecasting  Methods 
and  of  the  Height  of  a  Fully  Developed  Sea,"  Deutsche 
Hydrograph is che  Zeitschrift,  Vol ,  12,  no,  6,  1959,  P.  244. 

10 

Phillips,  0.  M.,  "On  the  Generation  of  Waves  by  Turbulent 
Wind,"  Journal  of  Fluid  Mechanics,  Vol.  2,  1957,  P.  417. 
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At  present  there  are  two  principal  theories  on  the  growth 
of  waves: 

(a)  Phillips'  resonance  theory.'1''1' 

12 

(b)  Miles'  sheer-flow  instability  theory. 

These  two  theories  do  not  agree  in  many  respects  (e.g„, 
Phillips'  theory  gives  a  "growth  rate"  proportional  to  time  and 
Miles'  gives  a  "growth  rate"  exponential  with  time).  While 
both  mechanisms  play  a  role  in  the  growth  of  waves,  Phillips  has 
recently  given  a  description  of  the  domains  of  dominance  of 
each  mechanism.  These  theories,  while  they  do  give  some  under¬ 
standing  of  the  sea  surface,  are  unfortunately  of  no  qualita¬ 
tive  help  at  present  in  describing  the  sea  surface . 

We  will  consider  only  one  other  mechanism  of  energy  trans¬ 
fer,  the  breaking  of  the  waves.  Under  sufficiently  high  winds, 
the  energy  transfer  from  the  breaking  of  waves  reaches  an  equi¬ 
librium;  then  for  that  range  of  frequencies  where  the  nonlinear 
(breaking)  effects  are  important  the  energy  spectrum  is  given 
by: 

E(co)  =  ag^<D-^  where  a  =7.4  10”^  (empirical  constant)  . 

The  same  basic  conclusions  (under  basically  the  same  restric¬ 
tions)  ware  obtained  by  Mikhailov,^  who  used  the  theory  of 

- XT - ~ 

Phillips,  0.  M.,  "Resonance  Phenomena  in  Gravity  Waves," 
Proc.  of  Symposia  in  Applies  Math.,  Vol .  XIII,  Amer.  Math. 
Society,  McGraw-Hill,  New  York,  1962,  p,  91  • 

12 

Miles,  J.  W.,  "Generation  of  Surface  Waves  by  Shear 
Flows,"  Proc.  of  Symposia  in  Applied  Math.,  Vol  XIII,  Amer. 

Math.  Society,  McGraw-Hill,  New  York,  1962,  p.  79* 

'!' -'’Mikhailov,  V.  I.,  "On  the  Theory  of  Scattering  of  Elec¬ 
tromagnetic  Waves  on  the  Sea-Surface,"  Bulletin,  Academy  of 
Science,  U.S.S.R.,  Geophysics  Series,  i960,  p.  8l8. 


turbulence .  There  is  close  agreement  between  these  theories 

[Phillips  has  E(a))  w  Mikhailov  cites  experimental  work 

which  gives  E(oo)  »  to  both  theories  are  within  these 

bounds.  Note  Neumann  in  this  range  has  E(cd) 

14 

Tick  has  considered  the  problem  of  non-linear  effects 
from  a  perturbation  point  of  view.  He  adds  to  the  "linear 
(-Gaussian)  spectrum"  a  correction  spectrum  due  to  the  non¬ 
linear  effects.  However,  the  statistics  of  the  "non-linear" 
spectrum  are  not  unknown. 

All  theories  (as  those  above)  which  do  not  directly  con¬ 
sider  capillary  waves  (waves  of  very  short  wavelength)  do  not 
hold  for  frequencies  in  the  “capillary"  range.  While  much  work 
on  capillaries  has  been  done,  at  present  there  is  nothing  avail¬ 
able  on  correlating  this  theory  with  the  overall  sea  surface. 

And  there  is  no  correct  information  available  on  the  "capillary" 
range  of  the  energy  density  spectrum.  Work  on  radar  return  and 
light  reflection  from  the  sea  surface  has  shown  that  the  slopes 
of  the  sea  surface  are  nearly  Gaussian  and  that  the  sea  surface 
curvature  is  highly  non-Gaussian .  The  usual  experimental  data  * 
taken  for  description  of  the  sea  surface  does  not  include  fre¬ 
quencies  in  the  "capillary"  range. 


- — - 

“  Tick,  L.  J.,  "A  Non-Linear  Randan  Model  of  Gravity  Waves 
I,"  Journal  Math,  and  Mech.,  Vol .  8,  1959,  P.  645 . 

Pierson,  W.  J.,  Jr.,  "A  Note  on  the  Growth  of  the  Spec¬ 
trum  of  Wind-Generated  Gravity  Waves  as  Determined  by  Non- 
Linear  Considerations,"  Journ.  Geophysical  Research,  Vol.  64, 
no.  8,  August  1959,  p.  IOO7 . 
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The  review  of  the  past  work  done  on  the  mathematical  des¬ 
cription  of  the  sea  surface  indicates  the  point  of  departure 
for  the  study  of  propagation  of  electromagnetic  waves  near  the 
air-sea  interface.  The  complexity  of  the  air-sea  interface 
alone  is  enough  to  make  one  resort  to  statistical  analysis? 
however,  the  fact  that  both  the  linear  model  and  the  experi¬ 
mental  data  give  a  Gaussian  process  for  the  gross  feature  of 
the  sea  surface  implies  that  a  useful  model  would  be  a  satis- 
tical  one . 

If  one  is  willing  to  neglect  the  non-linear  effects  of 
the  sea  surface  and  use  the  Gaussian  model,  the  Ergodic  theo¬ 
rem1^  implies  that  the  ensemble  analysis  of  statistical  pro¬ 
perties  is  the  same  as  space  (or  time)  analysis  of  the  statis¬ 
tical  properties  (note  that  this  requires  a  stationary  process, 
i.e.,  equilibrium  state).  There  are  questions  as  to  the  vali¬ 
dity  of  the  Gaussian  model  as  regards  the  statistical  proper¬ 
ties  of  the  sea  surface.  However,  if  one  is  only  interested 
in  the  gross  statistical  features,  then  the  Gaussian  descrip¬ 
tion  may  be  adequate.  However,  such  description  will  give  only 
the  "smooth  average"  shape  of  the  sea  without  any  details  such 
as  ripples  or'  minor  variations.  The  experimental  data  on  sea 
surface  variation  is  meagre?  the  processing  of  such  data  for 
meaningful  results  is  long  and  costly?  therefore,  it  is  doubt¬ 
ful  if  an  accurate  description  of  the  sea  surface  can  be 

- rr - 

^Lee,  Y.  W.,  Statistical  Theory  of  Communication,  John 
Wiley,  New  York,  I960 ,  p.  26J. 
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readily  obtained.  The  shape  of  the  surface  of  the  sea  depends 
upon  past  as  well  as  upon  present  conditions;  however,  to  a  de¬ 
gree  of  approximation,  past  effects  though  still  significant 
may  j  discounted.  Also,  by  describing  the  sea  condition  by 
means  of  only  a  few  average  parameters,  the  sea  surface  will 
not  be  represented  accurately  in  every  detail,  but  from  the 
practical  approach,  this  may  be  all  that  can  be  done. 


For  an  excellent  review  of  much  of  the  spectral  ana¬ 
lyses  of  the  sea  surface,  see  Proceedings  of  a  Conference  on 
Ocean  Wave  Spectra  (Englewood,  N.  J.,  Prentice-Hall,  1965) . 


Some  Useful  Formulae 


For  reference,  some  useful  relations  pertaining  to  the  sea 
surface  are  given  below. 

If  §(r)  =  A  cos  (it*r  -  ayt)  is  a  solution  to  the  linearized 
hydrodynamic  equations,  then: 

2 

a)  =  gk  tanh  k_h 
s  s 


for  kgh  >>  1  (large  depth) 

co2  =  gkg  (1) 

where 


We  then  have  for  a  single  harmonic  wave  in 
Vp  =  ^  ^  =  *2^  T  --  phase  velocity  of  the 


deep  water: 
sea  wave. 


v  = 


da) 

3k 


=  TT  V'  = 


-  T  --  group  velocity  of  the  sea  wave. 


L 


For  relatively  "regular"  seas  (after  they  have  consolida¬ 
ted  into  a  " regular"  series  of  connected  troughs  and  crests)  . 


H/height\  1  .  H  .  1 

tlerigtV  '•  TZ  <  L  <  35 

1  H  1  i 

(some  references  give  y  <  —  <  where  y  is  the  theoretical 
limit  for  stability  (greater  slopes  will  "break")). 


For  a  swell: 


1  <  H  <  1 
J5  K  L  <  "2DTT* 


Using  Neumann’s  spectrum  for  a  fully  developed  sea: 


E  2  .242^) 5 


v  --  knots 


Tave  -  -285  v 


T  --  sec. 


Lave  =  |(5.X2 )  T2, 


L  --  ft. 


'Pierson,  W.  J.,  G.  Neumann  and  R.  W.  James,  Observing 
and  Forecasting  Ocean  Waves,  Hydrographic  Office  Publication 


APPENDIX  B 


REVIEW  OF  THE  THEORY  OF  ROUGH  SURFACE  PROBLEMS 

In  this  appendix,  a  review  of  the  present  theory  of 
boundary  value  problems  involving  a  rough  surface  is  pre¬ 
sented.  Particular  attention  is  given  to  the  accuracy  and 
applicability  of  the  methods  presented  to  the  electromag¬ 
netic  wave- rough  surface  interaction  problem  considered 
in  the  main  sections  of  this  report. 

"In  the  mathematical  sense  the  problem  at  hand  is 
extremely  complex,  since  it  is  impossible  to  use  the  method 
of  separation  of  variables  to  obtain  a  solution  of  the  wave 
equation  which  satisfies  the  boundary  conditions  specified 
on  the  uneven  surface.  The  methods  for  obtaining  an 
accurate  solution  of  this  problem  in  general  form  has 
(sic)  not  yet  been  found.  However,  there  are  a  number  of 
theoretical  papers  (which  have  appeared  essentially  during 
the  last  five  to  six  years)  in  which  a  number  of  approxi¬ 
mate  methods  for  computing  the  field  have  been  developed." 
The  above  is  a  quotation  from  a  survey  paper,  "Theory  of 
the  Scattering  of  Waves  at  Periodically  Uneven  Surfaces," 
by  Iu.  P.  Lysanov."'"  This  paper  treats  for  the  most  part 
the  problem  of  a  scalar  field  in  the  presence  of  a  periodi¬ 
cally  uneven  (rough)  surface  with  "natural"  boundary 

■'"Lysanov,  Iu.  P.,  'Theory  of  the  Scattering  of  Waves 
at  Periodically  Uneven  Surfaces,"  Soviet  Physics  Acoustics, 
Vol.  4,  No.  1,  Jan. --March,  1958,  p.  1. 
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conditions.  When  the  boundary  conditions  of  the  actual 
problem  are  not  the  "natural"  ones,  the  "natural"  boundary 
conditions  are  used  as  an  approximation.  The  "interface" 
boundary  conditions  are  not  considered.  The  problem 
treated  is  then  the  simplest  problem  considering  a  rough 
surface;  even  then  this  problem  is  not  solvable  usr'ng 
known  techniques.  Although  work  on  this  problem  has 
continued  since  1958,  little  advancement  has  been  made 
and  no  basically  new  methods  or  theories  have  resulted. 

The  investigation  has  generally  been  confined  to  refine¬ 
ment  of  techniques  given  in  Lysanov's  paper  and  experi¬ 
mental  work. 

Basically,  Lysanov's  approach  will  be  followed.  He 
considers  six  general  classifications  or  techniques: 

1.  Rayleigh  Method.  This  is  the  oldest  approach 

and  due  to  Rayleigh.  Rayleigh  assumed  that  the  reflected 

(scattered \  field  could  be  represented  as  a  sum  of  outward 

(from  the  surface)  directed  plane  waves.  However,  this  is 

only  true  for  the  field  above  the  highest  point  of  the 

surface.  Lippmann  slowed  that  the  assumption  of  Rayleigh 

was  not  quite  correct.  Unfortunately,  he  could  not  obtain 

an  "accurate"  solution  to  the  problem  by  the  use  of  his 

2 

variational  technique.  Barantsev  has  also  obtained  the 
- 2 - 

Barantsev,  R.  G.,  "Plane  Wave  Scattering  by  a  Double 
Periodic  Surface  of  Arbitrary  Shape, "  Soviet  Physics 
Acoustics ,  Vol.  7/  No.  2,  Oct. -Dec.,  1961,  pT  125* 
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same  result  by  use  of  Laplace  transforms.  In  Barantsev's 
approach  to  the  problem,  as  in  Rayleigh's,  an  infinite 
number  of  algebraic  equations  must  be  solved.  In  the  case 
where  the  irregularities  are  small  (the  ma  imum  height  of 
the  surface  wave  is  much  smaller  _han  a  'ave  lergth  of  the 
radiation  field),  Rayleigh  obtainea  a  solution  for  the 
first  mode  (  ecular  reflection).  It  agrees  quite  well 
with  experimental  data  in  the  "far  field;"  however,  near 
the  interface  many  evanescent  (surface)  waves  are  present 
and  an  accurate  solution  is  not  practicable  (even  if 

■X 

Rayleigh's  assumption  was  considered  correct).-7 

2.  The  Method  of  Small  Perturbations.  The  boundary 
conditions  are  specified  on  the  uneven  surfa  e,  z  =  ( x ) , 

and  then  are  transferred  to  the  plane  z  =  0  by  means  of 
an  expansion  into  power  series  with  respect  to  § .  It  has 
been  shown  by  Lysanov  that  this  method  leads  to  solutions 
identical  to  Rayleigh's  results.  This  method  has  been 
used  by  Feinberg  for  solving  the  problem  of  propagation 

_ 

^La  Casce,  E.  0.,  "Some  Notes  on  the  Reflection  of 
Sound  from  a  Rigid  Corrugated  Surface,  "  Jourr^j.  >'£ 
Acoustical  Society  of  America,  Vol.  53/  No.  12,  Dec., 

1961,  p.  1772. 

La  Casce,  E.  0.,  B.  D.  McCombe,  R.  L.  Thomas, 
"Measurements  of  Sound  Reflection  from  a  Rigid  Corrugated 
Surface, "  Journal  of  Accoustical  Society  of  America, 

Vol.  53/  NoT  T2~,  Dec . ,  19613  p3  1766. 

Senior,  T.  B.  A.,  "The  Scattering  of  Electromagnetic 
Waves  by  a  Corrugated  Sheet,"  Canadian  Journal  of  Physics, 

voi.  57;  1959  p-  787  (ad  258  moyr 


of  electromagnetic  waves  along  the  earth's  surface  and 

h 

more  recently  by  Senior  who  considered  the  effect  of 

surface  roughness  on  the  reflection  of  electromagnetic 

waves.  Senior  is  able  to  replace  surface  roughness  by  a 

change  in  the  impedance  at  the  boundary.  The  surface 

5 

impedance  is  then  a  function  of  position.  Hessel  and 
others  have  also  considered  the  problem  of  varying  surface 
impedance?  Hessel  seems  to  account  for  Wood's  anomaly  this 
way.  However,  this  method  is  useful  only  in  the  far  field. 
Bass^  also  uses  the  method  of  perturbations  to  obtain 
results  parallel  to  Senior,  and  they  again  are  valid  only 
in  the  far  field.  It  should  be  stressed  here  that  the 
method  of  small  perturbations  may  be  in  error  in  the  near 


^Senior,  T.  B.  A.,  "Impedance  Boundary  Conditions  for 
Imperfectly  Conducting  Surfaces,"  Applied  Science  Research, 
Section  B.,  Vol.  8,  i960,  p.  4l8. 

Senior,  T.  B.  A.,  "Impedance  Boundary  Conditions  for 
Statistically  Rough  Surfaces,"  Applied  Science  Research, 
Section  B.,  Vol.  8,  i960,  p.  437* 

Hiatt,  R.  E„,  T.  B.  A.  Senior,  and  V'.  H.  Weston, 
"Studies  in  Radar  Cross  Sections  XL, "  --  "Surface  Roughness 
and  Impedance  Boundary  Conditions,"  University  of  Michigan 
Research  Institute,  Ann  Arbor,  Michigan,  July  i960  (un¬ 
published)  <, 

^Hessel,  A.,  A.  A.  Oliver,  "On  the  Theory  of  Wood's 
Anomalies,  in  Progress  Report  No.  19/"  R*  452.1961, 
Polytechnic  Institute  of  Brooklyn  (unpublished), 

AI)  256  809. 


Bass,  F.  G.,  V*  G.  Bocharov,  "On  the  Theory  of  Scatter¬ 
ing  of  Electromagnetic  Waves  from  a  Statistically  Uneven 
Surface,"  Radiotekhnika  Elektronika,  Vol.  J>,  No.  2,  1956/ 
p.  18 G. 


Therefore, 


7 

field  and  is  valid  only  in  the  far  field. 

3 

Winter,  who  applied  the  method  of  small  perturbation  to 

compute  the  fields  on  the  rough  sea  surface,  could  be  in 

error  in  his  calculations.  The  results  obtained  by 

Winter's  method  were  compared  with  the  results  obtained 

in  this  report  and  there  was  a  difference  (see  Section  2.6). 
q 

Lysanov^  has  extended  the  theory  of  small  perturba¬ 
tions  to  include  point  sources. 

3.  The  Method  of  L.  M.  Brekhovskikh .  This  could 
better  be  called  the  method  of  geometric  optics.  Basically, 
it  uses  Kirchhoff's  principle  (or  approximations).  For 
the  frequency  range  we  are  interested  in,  the  use  of 
Kirchhoff's  principle  could  well  be  in  error.  In  any 
case,  for  the  solution  to  the  problem  on  the  interface 
this  method  would  assume  away  the  problem.  The  problem 
becomes  one  of  evaluation  of  an  integral  which  cannot  be 
done  in  closed  form  near  the  surface. 


^Bass,  F.  G.,  "On  the  Theory  of  Combination  Scattering 
of  Waves  on  a  Rough  Surface, "  Izvestia  VUZ,  Radiofizika, 

Vol.  4,  No.  1,  1961,  AD  262  417. 

Lysanov,  op.  cit.,  p.  3- 

^Winter,  D.  F.,  ''Low  Frequency  Radio  Propagation  into 
a  Moderately  Rough  Sea, "  Radio  Propagation  --  Section  P, 
Journal  of  Research,  National  Bureau  of  Standards,  Vol.  6r[D, 

No.  5,  Sept.-0<-t.  1963,  p-  551- 

^Lysanov,  Iu.  P.,  "On  the  Field  of  a  Point  Radiator 
in  a  Laminar- Inhomogeneous  Medium  Bounded  by  an  Uneven 
Surface, "  Soviet  Physics  Acoustics,  Vol.  J ,  No.  3/ 

Jan. -March,  1962 /  p.  235 • 


4.  The  Integral  Equation  Method.  The  problem  of  the 
scalar  boundary  value  problem  with  "natural"  boundary 
conditions  can,  by  use  of  Green's  theorem,  be  formulated 
conveniently  as  an  integral  equation.  This  equation  is 
exact;  however,  it  cannot  be  solved  without  making  some 
approximations.  If  the  rough  surface  satisfies  certain 
conditions  (is  not  too  rough),  an  approximate  integral 
equation  is  obtained  which  can  be  solved.  There  is  a 
difference  in  the  field  assumed  on  the  surface  by  geometric 
optics  (Kirchhoff's  approximation)  and  the  field  on  the 
surface  obtained  by  the  integral  equation  method.  This 
has  led  Meecham^  to  question  the  use  of  Kirchhoff's 
approximation.  This  criticism  could  also  be  extended 
to  the  method  of  small  perturbations  in  the  near  field. 

A  variation  of  the  integral  equation  method  was  used  in 
the  main  part  of  this  paper. 

Recently,  a  paper  solving  the  exact  integral  equation 
for  a  sinusoidal  surface  was  presented; ^  the  numerical 
results  involved  approximate  solutions  to  an  infinite  set 
of  algebraic  equations. 


10 

Meecham,  W.,  "On  the  Use  of  the  Kirchhoff  Approxima¬ 
tion  for  the  Solution  of  Reflection  Problems,"  J.  Rational 
Mech.  Analysis,  Vol.  5,  1956,  p.  525'  . 


^Uretsky,  Jack  L. ,  "Reflection  of  a  Plane  Sound  Wave 
from  a  Sinusoidal  Surface, "  submitted  as  a  Letter  to  the 
Editor,  J.  Acoust.  Soc.  Am.,  March  1965- 
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5.  The  Method  of  Images.  The  method  of  images  can 
be  used  to  investigate  the  fields  in  the  presence  of  an 
'  uneven  surface  with  a  sufficiently  simple  shape;  this  has 
been  done  by  Twer sky  for  a  perfectly  reflecting  plane 
covered  with  ha If- cylinders  or  hemispheres  (with  little 
interaction  between  scatterers).  Biot  considered  the 
perfect  conduction  plane  covered  by  hemispheres  with 
strong  interaction  (as  there  would  be  in  the  case  of  the 
sea  surface).  Twersky  also  has  considered  multiple 
scattering  in  a  very  general  way.  This  method  could  be 
used  to  find  the  fields  on  the  surface  assumed  to  be  a 
plane  covered  with  hemispheres.  However,  to  obtain  a 
more  realistic  model  of  the  sea  surface  would  involve 
higher  multi-pole  expansions  and  does  not  seem  practical 
(although  Biot  thought  this  method  could  be  used  for  the 
air-sea  interface  problem).  It  does  give  a  good  idea  of 
what  is  happening  to  the  electromagnetic  field  in  the  air 
near  the  air-sea  interface  (Biot  uses  a  static  (co  =  0) 
solution) . 


T2  ! 

Biot,  M,  A.,  "Some  New  Aspects  of  the  Reflection 
of  Electromagnetic  Waves  on  a  Rough  Surface, "  Journal 
Applied  Physics,  Vol.  28,  December  1957/  P*  1455* 

Biot,  M.  A.,  "On  the  Reflection  of  Electromagnetic 
Waves  on  a  Rough  Surface,"  Journal  Applied  Physics, 

Vol.  29/  June  1958/  p.  998- 
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6.  Method  of  Matching  Fields-  This  method  can  be 
used  only  in  the  case  where  the  rough  surface  is  such  that 
the  space  may  be  separated  into  regions  in  which  the  wave 
equation  allows  solution  by  the  method  of  separation  of 
variables  (when  the  wave  equation  is  written  in  an  appro¬ 
priate  coordinate  system).  As  this  is  not  possible  for  the 
air-sea  interface,  we  will  not  consider  this  method  further. 

Lysanov  also  references  some  experimental  papers. 

Since  1958  other  experimental  papers  have  been  published 
which  seem  to  imply  that  Rayleigh's  theory  (including  the 
second  mode)  is  correct  for  the  low  frequency  problem 
except  near  Ar  =  A^  (Ar#  wavelength  of  radiations; 

Aw,  wavelength  of  the  surface).1-^  Since  1958  many  papers 
on  the  statistical  analysis  of  the  reflection  of  sound 
from  the  rough  sea  surface  have  been  published;  these  use 
either  Rayleigh's  assumption  or  Kirchhoff's  approximation 

for  the  solution  of  the  boundary  value  problem  and  as  such 

14 

are  only  useful  in  the  far  field. 


15 

La  Casce,  loc.  crt. 

14 

Proud,  J.  M.,  Jr.,  R.  T.  Beyer,  and  Paul  Tamarkin, 
"Reflection  of  Sound  from  Randomly  Rough  Surfaces, "  Journal 
of  Applied  Physics,  Vol.  51/  No.  5/  March  i960,  p.  545* 

Marsh,  H.  W.,  "Exact  Solution  of  Wave  Scattering  by 
Irregular  Surfaces,"  Journal  of  the  Acoustical  Society  of 
America,  Vol.  55/  No.  5/  March  1981,  p.  550* 

Marsh,  H.  W.,  M.  Schulkin  and  S.  G.  Kneale,  "Scatter¬ 
ing  of  Underwater  Sound  by  the  Sea  Surface, "  Journal  of  the 
Acoust.  Soc.  of  Amer.,  Vol.  55/  No.  5/  March  1951,  p.  55^. 

Beckmann,  Petr,  and  Andre  Spizzichino,  The  Scatter¬ 
ing  of  Electromagnetic  Waves  from  Rough  Surfaces,  Pergamon 
Press  Ltd.,  London,  1965- 
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One  of  the  few  papers  considering  "interface"  boundary 
conditions  on  an  irregular  surface  was  published  by  Wait  J 
in  1959*  this  paper  has  been  used  to  calculate  the  effect 
of  the  rough  surface  on  the  electromagnetic  fields  in  the 
sea.  Wait's  basic  approach  is  the  work  of  Leontovich  on 
approximate  boundary  conditions  for  a  good  conductor. 

Wait  obtains  the  field  on  a  plane  in  the  sea,  and  solves 
the  wave  equation  in  the  sea  with  this  as  a  boundary  condi¬ 
tion.  A  short  discussion  of  this  work  is  given  in  a  paper 

1 /T 

by  Whalen.  Some  comment  on  the  accuracy  of  Wait's 

assumption  is  given  in  the  main  part  of  this  paper.  Other 

work  directly  considering  the  rough  sea  surface-VLF  radio 

17  3  8 

wave  interaction  are:  Lerner  and  Max  '  and  Winter.  '  Winter 
used  the  method  of  small  perturbation  to  compute  the  fields 
on  the  sea  surface  and  the  stochastic  Stratton-Chu  integral 
equations  to  obtain  statistical  results  for  the  fields  in 
the  sea.  The  use  of  the  method  of  small  perturbation  is 


^Wait,  J.  R.,  "The  Calculation  of  the  Field  in  a  Homo¬ 
geneous  Conductor  with  a  Wavy  Interface, "  Proc.  IRE,  Vol.  47, 
No.  6,  June  i960,  p.  1155- 

16 

uWhalen,  J.  L.,  "Measured  Effects  of  Ocean  Waves  on 
the  Phase  and  Amplitude  of  VLF  Electromagnetic  Radiation 
Received  Below  the  Waves,"  USL  Tech.  Memorandum  No.  9^'1~ 
67-61,  5  August  1961  (unpublished). 

17 

'Lerner,  R.  M.  and  J.  Max,  "Very  Low  Frequency  and  Low 
Frequency  Fields  Propagating  Near  and  Into  a  Rough  Sea, " 
a  paper  presented  to  the  URSI  Spring  1965  Meeting. 

^Winter,  loc.  cit. 


discussed  in  Section  2  of  this  appendix  and  Winter's  statis¬ 
tical  results  are  considered  in  the  conclusion  chapter  of  the 
main  part  of  this  report.  The  basic  method  of  Lerner  and 
Max  was  outlined  in  Section  2.6  of  this  report,  and  their 
results  are  considered  in  the  conclusion  part  of  this  report. 


APPENDIX  C 


SINGULAR  INTEGRALS 


Although  the  results  given  in  this  appendix  are  avail¬ 
able  in  the  literature,  it  seems  useful  to  make  readily 
available  to  the  reader  the  basic  theory  and  results  per¬ 
taining  to  "singular"  integral  equations.  These  results 
have  been  used  in  the  body  of  the  paper  to  obtain  the  inte¬ 
gral  equations  of  Sections  1.4  and  5-2. 

The  starting  point  in  considering  "singular"  integrals 
is  the  concept  of  improper  integral.1  By  "singular"  inte¬ 
gral  is  meant  an  integral  whose  value  changes  discontinu¬ 
ous  ly  (i.e.,  a  "discontinuous  integral"). 

The  field  of  a  single  layer  is: 

F  (r)  ^  f  cp  (r  '  )  G  (k  |  r-r '  |  )dS  '  (C.l) 

s  J 

S 

r  e  r  '  e  S 

V.  -  interior  volume  bounded  by  closed  surface  S. 

1 

Vg  -  exterior  volume  bounded  by  closed  surface  S. 

<p(r)  -  strength  of  the  single  layer  sources. 

F  (r)  -  field  of  the  single  layer  sources, 
s 


"Slikhlin,  S.  G.,  Integral  Equations,  Pergamon  Press, 
London,  195?/  P-  113- 
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lim  F  (r )  =  lim  f  [*  <p  (r  ’  )  G  (k  |  r-  r '  |  )dS  1  + 

r-S  s  r-S  L  °S-N 

N  -*0  e 

e 

cp(r'  )  G  (k  | r-r '  |  )  dS  '  J  (C.2) 

e 

N  -  neighborhood  of  point  r  in  S  (in  this  neighbor- 
hood  G(k|r-r'|)  becomes  unbounded  and  the  integral  is  improper). 
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aS  ’  =  lim  cp(r)  j  G(kr’ )  dS  '  = 

N  ,  -O  4  . 
e  e 


lim  cp  (r }  G  (r '  )  dS  ’  =  0 

N  ,  -0  4  , 


as  G  is  an  even  function. 

N  t  -  neighborhood  of  point  r ‘  in  S . 

6 

G(r)  -  static  Green's  function  (the  approximation 

G(r')  fa  G(kr')  is  valid  for  kr '  «  l),  which 

holds  in  N  , . 
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<p(r '  )  G(k  j  r-r 1  j )  dS ' 
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(c-5) 
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The  field  of  a  double  layer  is  represented  by  a  singu¬ 
lar  integral. 


Fd(5) 


MS 


<p(r')  G(k|?-r'  I)  dS’ 


(c.4) 


<p(r)  -  strength  of  double  layer  sources, 

n  -  outward  normal  to  S. 

F,  (r)  -  field  of  double  layer. 
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n  -  "solid  angle"  of  surface  S  —  (for  smooth  surface 
0  =  2  tt) 
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Courant,  R.  and  D.  Hilbert,  Methods  of  Mathematical 
Physics,  Vol.  2,  Interscience,  New  York,  1962,  p.  253- 

-'Wikhlin,  S.  G.,  Linear  Integral  Equations,  Delhi- 
Hindustan,  i960,  p.  I56T 
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sJy  equations  (C  *3/  Gild  (C.6). 

Equation  (c.8)  is  the  equation  used  in  Sections  1.4 
and  5-2. 

The  extension  of  these  results  to  the  coupled  electro¬ 
magnetic  equations  (stratton-Chu)  is  straightforward.^ 


- n - 

Honl,  H.,  A.  W.  Maue,  and  K.  Westpfahl,  "Theorie  of 
Beugung,"  in  Handbuch  Per  Physik,  Band  XXV/1.  (Berlin: 
Springer,  196 1 ) ,  p .  2 l8T 
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